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Abstract. We prove that the charge-scalar field equations are globally stable 
on (3 + 1) dimensional Minkowski space for small initial data in certain gauge 
covariant weighted Sobolev spaces. These spaces can be chosen to be almost 
scale invariant with respect to the homogeneity of the equations, and our 
result is valid for initial data with non-zero charge that is also non-stationary 
at space-like infinity. The method of proof is a tensor-geometric approach 
which is based on a sharp family of weighted bilinear L 2 space-time estimates. 



1. Introduction 

In this paper, we study the global in time behavior of small amplitude clas- 
sical solutions to the massless Charge-Scalar-Field equations (CSF) equations on 
Minkowski space. These are also sometimes referred to in the literature as Maxwell- 
Higgs or the Maxwell Klein Gordon equations. They are constructed as follows: 
We let M. denote the (3 + 1) dimensional Minkowski space with metric g = 
(—1,1,1,1) and compatible Levi-Civita connection which we refer to here as V. 
Then let: 

V = MxC, 

denote a complex line bundle over M. with hermitian inner product (v)y an d 
compatible connection D. That is, one has the formula: 

(1.1) X(ip,if>)v = (Dxr/>,4>)v + (il>,Dx<l>)v , 

for vector-fields X on M.. The connection D can be extended in a natural way to 
sections of the tensor bundle V ® TfM. of complex valued tensors of type (k,l) on 
Minkowski space such that one has the Leibnitz rule: 

D x {<f)®T) = Dx(c/))®T + <j)®Vx<j) ■ 

In this setup, we call a section cj> of V a complex scalar field if it satisfies the 
equation: 

(1.2) n c </> = D a D a (/> = , 

where D 2 <f) is the covariant Hessian of <fi, which is computed regarding D<p as a 
complex valued one- form on M. Since V is globally trivial, we can let ly denote 
a unit normalized global section. That is one has (ly, ly) = 1. Such a choice is 
referred to as a gauge. In terms of this, any section 4> to V can be identified with a 
complex valued function on M. where we abusively write <p = 4>lv- This allows us 
to define complex conjugation <f> as one would for ordinary complex numbers. We 
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will use this operation to define a symmetric and symplectic inner product on V as 
follows: 

R<&$>v = -(# + #), %4>,^)v = 

Also, with respect to this frame for V, and for any frame {e a } on M., we define the 
connection one-form {A Q } via the relations: 

(1.3) D a l v = s/=lA a . 

Notice that the compatibility condition immediately implies that the {A a } 

are real. Using this notation, the covariant derivative of any complex valued tensor 
T with respect to a vector-field X in M. can be written as: 

(1.4) D X T = V X (T) + V^1A(X)-T . 

Of course, the choice of ly and hence representation of D x is somewhat arbitrary 
because one can always perform a local (or global) unitary transformation of it. If 
we let ly ly denote the transformation given by ly = e~ Ix ly, where x is some 
real valued function on A4, then 11. 4|) shows that the potentials {A a } transform 
as: 

A a = A a — Sa{x) ■ 

Now it turns out that this ambiguity which is inherent to the equation (|1.4f> does 
not need to be resolved in order to proceed with a detailed analysis of the scalar 
field equation l|1.2f> . This is because the connection (|1.1|) has a basic geometric 
invariant on which it is possible to base literally all analysis of l|1.2[) so long as one 
is content to define all analytic objects of interest in terms of geometric quantities. 
This is the approach we shall take here. The basic geometric invariant is, of course, 
the curvature of the connection D which arises from the operation of commuting 
covariant differentiation. Specifically, there exists a (real) two-form F on A4, such 
that for any two vector-fields X, Y we have the relation: 

(1.5) D x D Y <t>-D Y D X (t>-D [Xj Y](j> = V^1F(X, Y) ■ <j> . 
In terms of the frame {e a } we can simply write: 

(1.6) D a Dp<f) — DpD a (j) = 

Here D a Dp(j) denotes the (a, (3) component of the Hessian D 2 <f> and should not be 
confused with the repeated directional covariant differentiation D ea D e/3 (j). To see 
the difference, a short computation of l|1.5|) in the frame {e a } shows that one has 
the identity: 

De a D e0 <j>- D e0 D ea <f>- D[ eate ^(p = v /z T(e Q (A l a) - ep(A a ) - c^A T ) • <f> , 

where [e m e^] = c^e 7 are the structure "constants" of the frame {e Q }. In terms 
of the indexed version l|l.f)|) . using (|1.4|l this last line reads: 

D a Dp4>-DpD a <l> = y/^l(y a A fi - VpA a ) ■ <f> , 

Either way, this gives the well known "Bianchi identity" : 

(1.7) F = dA. 



To define the full CSF system, we now couple the complex scalar field <f> to the 
curvature of the connection D in such a way that one ends up with a Lagrangian 
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field theory. For an intuitive approach to how this is done, we consider F as a 
solution to Maxwell's equations. Because F satisfies the Bianchi identity (|1.7|l . we 
have that: 

(1.8a) V p F a0 = J a , 

(1.8b) V p *F a p = , 

for some one-form J which obeys the rule V Q J a = 0. In the above formulae *F 
denotes the Hodge dual of F which is given by the expression *F a p — \ G Q/3 7<5 F 1 g. 
Here Ga/3-yS denotes the volume form on Minkowski space. As usual, we use this 
duality operation and contraction along the time-like vector-field dt to define the 
electric and magnetic field strengths: 

(1-9) Ei — Foi , Hi = *F 0i . 

Using these quantities, the system <|1.8[) takes the following familiar form which will 
be of use in the sequel: 

(1.10a) d t E - V x x H = J , (div)J5 = J Q , 

(1.10b) d t H + V x x E = , (div)ff = . 

Here V x x =s/ dj is the usual curl operation, and J denotes the spatial part of 
the current vector J. We note here that this latter form of the equations will be 
particularly useful when dealing with questions concerning the initial data of Lie 
derivatives of F a p . In general it is the more natural language with which to discuss 
the initial value problem for the system l|1.8|) . 

For the Maxwell field F, there is an energy- momentum tensor Q[F] which is 
given by: 

(1.11) Q afi [F] = I(F Q7 F^ + *F Q7 % 7 ) = F ai F^ - ^g aP F yS F^ 5 . 

A brief calculation using the identities (|1.8a|) - (|1.8bp shows that one has the diver- 
gence identity: 

(1.12) V a Q al3 [F] = - F^F . 

Furthermore, the complex scalar field i|1.2fl also has an energy-momentum tensor 
which is analogous to the usual energy-momentum tensor of the (non-covariant) 
D'Lambertian. This is given by: 

(1.13) QaM = ^{D a ^4>)-^g a pD^<fD~^ . 

A direct application of the field equation (|1.2fl . the compatibility condition 
and the commutator identity (|1.5f) shows that one has the divergence law: 

(1.14) v a Q af }[4>] = F^{^W4>) . 

Notice that this formula does not contain any explicit reference to the connection 
coefficients <|1.3H • This is because it is a purely tensorial identity, and in particular 
does not depend on a choice of coordinates or frame. 

The natural choice of coupling F to (j) now comes if we simply stipulate that the 
current vector on the right hand side of (jl.8a|l be given by J a = ^((f)D a <f)), With 
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this extra condition satisfied, the total energy momentum tensor Q = Q[F] + Q[<j>] 
becomes divergence free: 

(1.15) V Q Q Q/3 = . 

Notice that this coupling is legitimate because as an immediate consequence of the 
equation (|1.2|l one has that this choice of current vector J a satisfies the continuity 
equation: 

(1.16) V a J a = , 

which is the only prerequisite for a vector J to show up as the right hand side of 
the Maxwell equation l|1.8af> . We write the resulting system of equations, which we 
henceforth refer to as the CSF equations, together as: 

(1.17a) V^F af3 = 3(0A^) , 

(1.17b) n c = 0. 

We note here that implicit in the system l|1.17fl is the Bianchi identity (|1.7fl or 
(|1.8b|) . This is because we are assuming a-priori that F is the curvature of the 
connection which gives rise to O c . 

In this work we study the global in time Cauchy problem for the system H1.17|) . 
Since we have stated the equations in such a way as to deemphasize a choice of 
gauge, this deserves some explanation. Because we are interested in an evolution 
problem, we are reduced to a discussing the notion of gauge covariant initial data. 
To specify these, we first let D_ denote a connection on the initial time slice bundle 
{0} x R 3 x C. Since this is embedded as a hyper-surface in the original bundle 
M. x C = V, we will also specify the initial normal derivative D t . The initial data 
for the system (|1.17|l can then be written in the form: 

(1.18a) F Oi {0) = E, , *F 0i (Q) = H t , 

(1.18b) 0(0) = O , A0(O) = O • 

In the above notation, we have used the labels (E, H) to denote quantities which 
only depend on a; £ K 3 . These should not be confused with the space-time (E, H) 
used in the formulas (|1.1U|) . Now, from the form of the system i|1.17|l it is easy 
to see that this initial data cannot be specified freely. It must also satisfy the 
compatibility conditions: 

(1.19) = 3(0o ^o") , V*fli = 0. 

We will call a data set (E, H, 0o, 0o) which satisfies <|1.19f) admissible. The question 
we will be concerned with here is to describe in as detailed a way as possible the 
global in time behavior of solutions to the CSF equations whose initial data satisfies 
certain natural smallness assumptions. As usual, these will be stated in terms of 
regularity. In general, we define the covariant and gauge covariant weighted Sobolev 
spaces: 



(1.20) 
(1.21) 
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Here / denotes the usual multiindex notation, while V x denotes the restriction of 
the Levi-Civita connection V to the time slice {0} xll T denotes an arbitrary 
tensor. We will consider the Cauchy problem for admissible initial data sets which 
are in these function spaces for various values of So and k. 

Now, a basic regularity theorem for the system (|1.17(l says that if the norms 
(|1.20|> - (|1.21fl applied to (admissible) [E, H, </>o, <^o) are finite for certain values of 
so and fc, then a global solution to this system exists. In fact, there is no need to 
impose any of the weights (1 + r 2 ), and almost no smoothness is needed to make the 
argument work. This is the content of the fundamental regularity result of Eardly- 
Moncrief and its later significant refinement due to Klainerman-Machedon 

[§]. While these results in some sense give the strongest possible global existence 
result one could hope to ask for, they contain surprising little information as to the 
nature of the global solution they obtain. In fact, the only information they provide 
is that the unweighted Sobolev H s norm of the solution remains bounded for all 
time. No practical bounds for this quantity are obtained for 1 < s. Furthermore, 
these results provide no information on the profile of the solution, even assuming 
that the initial data is localized at time t = in a way that is consistent with the 
norms H1.20|l - I|1.21(l . Finally, these results do not lend themselves to a discussion 
of the asymptotic behavior of the scalar field <f> in the sense of L 2 scattering. For 
example, there is no doubt that this involves long range corrections 1 due to the 
asymptotic behavior of the electro-magnetic field F a p, but the large data results 
just mentioned provide no information as to what form these corrections should 
take. 

The result we present here is a first attempt to fill in some of the gap in un- 
derstanding the asymptotic behavior of solutions to the system (I1.17f) . at least for 
initial data with small (|1.20() - (|1.21|l norms where the specific values of sq and k are 
to be specified. The method we employ is a new variant of the tensorial-geometric 
approach of Christodoulou-Klainerman 0-01 which also uses in a crucial way the 
space-time energy estimate approach appearing in the recent work of Lindblad- 
Rodnianski |12|. This will be applied to both the electro-magnetic field (|1.8|l and 
to the scalar field l|1.2|l through a novel geometric adaptation of a set of fundamen- 
tal estimates which goes back to work of Morawetz The approach we follow 
also borrows ideas from previous works on the low regularity behavior of systems 
of semilinear wave equations, notably 9 , in that we make crucial use of a certain 
family of sharp weighted bilinear L 2 space-time estimates for tensorial contractions 
to control error terms which come up in our analysis. Some of these estimates are 
proved through the use of a weighted L 2 (L°°) "Strichartz" type estimate which 
holds for certain null-components in a tensorial setting. 



We will not discuss these corrections here, although they are implicit in the gauge-covariant 
approach we use. We believe the analysis we present goes a long way towards constructing and 
estimating the effect of long range corrections in the system 11.171 . which should take the from of 
phase corrections which depend only in the charge. Also, we note here that the existence of long 
range corrections to the CSF equations is consistent with what is known for the Maxwell-Dirac 
system (see |7|) which is similar. 



6 HANS LINDBLAD AND JACOB STERBENZ UNIVERSITY OF CALIFORNIA AT SAN DIEGO 



There has been previous work in general 2 on the global asymptotic behavior of 
the equations (|1.17(l . This is the paper of Shu ^Sj- That program, while similar 
in spirit to ours, was not carried out in detail and was only aimed at resolving the 
case of solutions which are trivial outside of a fixed forward cone. 3 It should be 
noted that we have found no way of closing our argument without making crucial 
use of space-time L 2 (L°°) estimates which serve to eliminate certain logarithmic 
divergences one would encounter when trying to pursue an L°° approach. As far 
as we can see, these types of estimates are necessary even of one assumes compact 
support for the initial data (modulo charge) . There is no discussion of space-time 
norms in the work |15| . or any other mechanism which would allow one to circum- 
vent divergences arising in the purely L°° and energy estimates approach taken 
there. However, we would like to call the readers attention to the fact that the 
work |15| represents in many ways the approach we will take here. In particular, 
some of the more interesting features of the CSF nonlinearity were first discovered 
there, notably the beautiful cancelation which we use in (a slightly different form 
on) line (|9.11|l below. 

Finally, we also mention the work of Psarrelli ^3 011 the asymptotic behavior 
for the equations i|1.17[) with a non-zero mass. This work again assumes compact 
support of the initial data, and is ultimately very different form what we do here 
because the version of l|1.2(l with non-zero mass enjoys a quite different decay es- 
timate than its massless counterpart. It would be an interesting and non-trivial 
problem to extend the work of ^11 to the case of fields which are non-stationary 
at space-like infinity, and it is not clear if our method is a step in this direction 
because of our persistent use of estimates derived from conformal invariance. 

Before moving on the statement of our main result, it is worth mentioning at 
a heuristic level two of the most difficult features of attempting a purely physical- 
space analysis of the equations (|1 . 1 T|l . We hope this gives the reader a bit of insight 
into why certain things are necessary. The first problem has to do with attempting 
to make use of the so called "null-condition" inherent in the system <|1.17ll . Unlike 
other semilinear geometric wave equations, for example the wave- maps |17|. there 
is not a direct correspondence between the semilinear model equations studied in 
[Hj and which satisfy the null condition directly, and the equations of gauge field 
theory. 

There are basically two known ways the null condition for the equations (|1.17(l 
can be uncovered and utilized. The first is through the use of an elliptic gauge, as 

2 "General" meaning the system 11.171 with or without charge. There is also the previous 
work of Christodoulou-Choquet-Bruhat (|2]) which uses conformal compactification. However, 
this type of procedure requires more decay of the initial data and does not allow for either charge 
or dipole moments. 

3 While it is usually the case with scalar wave equations that once one can prove global existence 
using weighted energies for compactly supported initial data, one can move to non-compactly sup- 
ported initial data by simply adding enough weights (at space-like infinity), for tensorial equations 
of the type 11.171 it is not possible to directly pass from the assumption of compactly supported 
initial data to the case of initial data with infinite support. This problem persists even if one 
assumes that the initial data decays at a very fast polynomial rate. The reason is not entirely 
obvious at first, and is due to subtle interaction between the scaling and tensorial properties of 
the equations 11.171 . This will be explained in a moment. 
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in the works on the low regularity properties of this and other related systems Uj- 
[TT)| . This type of procedure leads to non-local versions of the null- forms studied in 
[5] and PP which also have entirely different scaling properties than their non-local 
versions. To analyze the system of equations which arise in this way by using the in- 
homogeneous algebra (|2.28() is quite awkward, and it is not even clear that this can 
even be done correctly. For example, the commutator of the weighted derivatives in 
(|2.28|l with Riesz potentials seems to cause enough of a problem that it effectively 
negates the savings one gets from the standard null-forms. Furthermore, any such 
analysis would have to somehow take into account the long range corrections the 
system (|1.17|) must undergo due to the lack of decay of the electro-magnetic field 
(|1.8(l . These corrections are clearly a guage dependent phenomena, and it does not 
seem that the usual elliptic gauge is most convenient for writing them down. 

The second way the null condition for the CSF system can uncovered is directly 
through its tensorial structure. That is, the null condition makes itself evident 
through the contraction structure of error terms which arise after commuting the 
equations with geometrically defined operations of differentiation. To be useful, the 
geometric differentiation must take the form of Lie differentiation of the Maxwell 
field and gauge covariant differentiation of the scalar field. Because the resulting er- 
ror terms then involve contractions over all possible indices, they can be expanded 
in a frame which takes into account the null geometry of Minkowski space. In 
this way, one immediately sees that two "bad components" can never interact with 
each other. However, since the system H1.17|) involves the interaction of two differ- 
ent types of quantities, namely a two form and a scalar field, it is not entirely clear 
at first that this type of tensorial structure will cause enough "good interactions" to 
take place, or even what the "components" of the scalar field should be. Fortunately 
for us, the CSF equations do contain a deep underlying structure which allows one 
to treat both quantities, F a p and (j), as a single instance of a master "tensorial" 
object. What's more, one can produce a single family of bilinear estimates for this 
"tensor" which covers all possible error terms that can generated by differentiating 
the system (I1.17|) . This type of underlying unity goes well beyond the scalar null 
structure of |S] and p^, and is akin to what one sees for the pure Yang-Mills and 
the Einstein equations (see |16|-P] ). 

The second main difficulty we need to overcome in our analysis is a consequence 
of an interesting interplay between the scaling and decay properties of solutions 
to the system (|1.17ll . While in (3 + 1) dimensions these equation are subcritical 
with respect to the translation invariant conserved energy 4 , which ultimately is 
responsible for large data regularity as mentioned previously, these equations are in 
general critical with respect to decay at space-like infinity. To better understand 
this, consider a linear wave equation of the form: 



4 The conserved energy is a simple consequence of the tensorial conservation law 11.151 and the 
fact that dt is a Killing field on Minkowski space. This conserved quantity is at the level of the 
H 1 norm of <j>. 

5 This should really be called "charge critical" . The conservation of charge for these equations, 
which can be written as q = J r3 ^s(<t>Dt4>) = const., takes place at exactly the critical regularity 
. i 

which is the H 2 norm of rp. 



(1.22) 
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where □ = d a d a denotes the usual D'Lambertian on Minkowski space. This can be 
seen as a model for the equation ljl.2|) where we have fixed a gauge (which gauge 
does not matter for the sake of this discussion). Now, assume that the initial data 
for i|1.22fl decays like r~ m as r — > oo. The question we then ask is how much decay 
should one then expect from the potentials A a in order to guarantee that the solu- 
tion keeps this decay in the far exterior region 2t < r (which should be the simplest 
to control!). By homogeneity, one expects that if A a decays like r~ l for 2t < r, 
then the solution will decay like r ~ m ~ l+1 (integrate the second order equation 
twice). 6 

One immediately sees from this simple analysis that to control things A a must 
decay at least as well as r _1 in the far exterior, and that for this critical rate of 
decay some work is required to avoid logarithmic divergences in the "error" term 
on the right hand side of (|1.22[1 . This decay rate r -1 is not a coincidence, but 
rather a consequence of the fact that the equation (|1.22|) scales like H 2 at the level 
of initial data, for which one sees that r _1 is precisely the critical rate of decay. 
What's more, the entire system (|1.17|l scales like H 2 at the level of the potentials 
A a , and because of the tensorial nature of the CSF equations at least some of the 
potentials must in general decay no better that r . This is a simple effect of the 
elliptic constraint equation: 

(1.23) V l E t = 3(0 A0) . 

Introducing a potential function if for the curl- free part E c f , we see that unless one 
has: 

f 3(0 A0) = , 

JR 3 

it will be the case that (p ~ i. Since this potential function must show up some- 
where in field potentials A a regardless of the gauge, we see that in general an 
equation of the form 11.22fl , and hence 1)1. 2f> , will be critical for decay at space-like 
infinity. 

To combat this critical rate of decay, one must prove sharp space-time estimates 
in order to have an effective strategy. There will in general be no extra "convergence 
factors" with which one can use to integrate over large time intervals. In practice, 
this means that all of our space-time estimates need to be proved in a-priori (diver- 
gence) form. This kind of difficulty should be understood in contrast to other field 
equations, such as the Einstein equations. In (3 + 1) dimensions the long range 
effect of the mass for this latter system leaves plenty of room with respect to the 
scaling properties of the equations (see for example J2]). Therefore, from the point 
of view of decay, the Einstein equations are somewhat more forgiving and do not 
require estimates which are as precise as what we need to make things work here 
(see for example [12] ). 

We are now ready to state our main result: 



For smooth A a this heuristic can be made precise through an easy use of weighted exterior 
energy estimates. 
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Theorem 1.1 (Global Stability of CSF Equations). Let 2 ^ k with k E N, and let 
sq = s + 7 be given such that so < | , < j, and | < s. Let (E,H,(f>o,(j)o) be an 
admissible initial data set, and define the charge to be the value: 



(1-24) q = / 3(0o fa) ■ 

Then there exists a universal constant £k,s,-y 7 which depends only on the parameters 
k, s, 7, such that if (E, H, cj>o, cfto) is an admissible initial data set which satisfies the 
smallness condition: 
(1-25) 

II E<li ll/f fc ' s 0(R3) + || H ||_H- fc . = 0(E3) + II D.4> ||ff">.«0(R3) + || 4>0 ll_ff fc -0(R3) < £fc,S,7 , 

where E = E d $ + E c $ is the Hodge decomposition of E into its divergence free 
and curl free components ( resp.), then there exists a ( unique ) global solution to the 
system of equations (|1.17|) with this initial data set such that if {L,L,ca} denotes 
a standard spherical null frame (see 1(2. 18|) and 1(2.19(1 below), then the following 
point-wise properties of this solution holds: 

(1.26a) 

_ s _3 _1 — s— i —A 

(1.26b) 

-2 -1- - 1 - 1 

\P\ ^ 9 ' 71 Xl<t<r + £k,s n ■ T + S T_ 2 • (w)j 2 , 

(1.26c) 

and: 
(1.27a) 

\D L <t>\ < £k, s ,<y ■ r+ s_i • (w)^ , \Dl4>\ < £ Mi7 • r+Vl'"* • , 
(1.27b) 

IWI < ^,a, 7 • r+^vi* • , 101 < ^ )S , 7 • T+vr s+2 " • m; 1 . 

_ffere we have set: 
1 



I^W - |-£>iW)rXKt<2r+|r>L^rXr<^ , IWI = 6 AB D A cj>D , 

and D denotes the spatial part of the connection D. Also, (a, a, p,o~) denotes the 

components of the null decomposition (12.23(1 of F a p . Finally, the weight functions 
t± and w 7 are defined via the formulas: 

rl = 1 + (f + r) 2 , t 2 _ = 1 + (t - rf , w 7 = t 2 J X t<r • 



Remark 1.2. In the course of our proof of Theorem 11.11 we will also show de- 
cay (peeling) properties similar to ((1.26(1 - 111.27(1 for the higher derivatives of (F, <fi) 
assuming that 2 < k. We have not stated this formally for the sake of brevity. 



Remark 1.3. Note the first term in the decay asymptotic for p on line (jl .26b|) above. 
As r — > 00 on any given time slice t = const, this portion of the electro-magnetic 
field decays only like r -2 . This is the long range effect of the electro-static equation 



10HANS LINDBLAD AND JACOB STERBENZ UNIVERSITY OF CALIFORNIA AT SAN DIEGO 

(|1.23|) . What is somewhat surprising is that this long range effect only propagates in 
the exterior of the forward cone t = r. That is, after rescaling at infinity, the effect 
of the charge is discontinuous across this cone. This fact goes a long way towards 
explaining why it is so much easier to deal with the case of compactly supported 
initial data. Also, notice that in the case where the value (|1.24|) vanishes, this long 
range effect is eliminated. This, of course, is the main reason why the assumption 
of zero charge is amenable to other techniques, e.g. conformal compactification. 

Remark 1.4. Note that we can allow any value so < I of the weight factor in the 
norms (|1.2U|) - (|1.21[1 applied to the initial data so long as it is bigger than the scale 
invariant factor of 4. At this value, the norms become homogeneous (except for the 
extra term of 1 in the physical space weight) with respect to the scaling properties 
of the equations l|1.17fl ■ It is also at this point that essentially all of our estimates 
break down due to logarithmic divergences. Also, for values of sq < 1, note that 
our condition on the initial data ljl.25jl requires much less decay than previous 
works based on weighted energies, e.g. the classical results in jg] for semilinear 
equations. 7 Thus, our main theorem falls just short of a scale invariant result, and 
can be considered a "low regularity" theorem with respect to physical space decay. 
Finally, we should mention that below the level of the scaling (i.e. sq < ^), there 
is no reason for the type of analysis we do here to make sense. In fact, even the 
simple notion of the charge (|1.24(l cannot be defined in this case because this level 
of decay is consistent with the asymptotics 0o ~ r ~ 1+s and 0o ~ r ~ 2+s for some 
< 8, which lead to a divergence in ljl.24|l . 

1.1. A brief outline of the work. We give here a very quick overview of the 
paper. In the next section, we make a list of some more or less standard geometric 
formulas relating to tensors and complex line bundles, as well as to the geometric 
structure of Minkowski space. This section is probably best left as a reference, and 
can be avoided by the reader with a passing acquaintance of this material. 

In the third section, we build a series of weighted L 2 estimates for electro- 
magnetic fields. This material is for the most part a straight-forward adaptation 
of the standard material found in with some new technical devices which are 
crucial for the approach we take in this paper. These are: A fractionally weighted 
version of the usual Morawetz type estimate for null-decomposed electro-magnetic 
fields, the decomposition of an electro-magnetic field into its pure charge and charge 
free components, and fixed time, characteristic, and space-time estimates which re- 
spect this decomposition. These estimates are proved with the help of a certain 
weighted elliptic estimate which we demonstrate in the appendix. 

In the fourth section, we set about proving analogs of the estimates of the third 
section for complex scalar fields. This involves an interesting new proof of the clas- 
sical Morawetz energy decay estimates for scalar fields which is directly based on 
the conformal geometry of Minkowski space. With some work, this leads to a set of 
estimates which are virtually identical with what is available for electro-magnetic 
fields. This device is important for what we do here because it allows us to treat 



7 For instance, we take the gradient D_<j> first before applying the weighted L norm. 
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the curvature and the (conjugated gradient) scalar field on the same footing. 

In the fifth and sixth sections, we prove L°° type estimates for the electro- 
magnetic and complex scalar fields. Interestingly enough, these involve in a crucial 
way the space-time, fixed time, as well as characteristic energy estimates developed 
in the preceding two sections. 

In the seventh section, we recast all of the estimates we prove in this paper in 
a uniform form. First we introduce an auxiliary "tensor" which has the behavior 
of the electro-magnetic-complex-scalar field which we are studying. We then prove 
a series of weighted L 2 bilinear estimates for interactions of various components of 
this "tensor" . It turns out that all of the estimates we need to control the error 
terms of the sequel can be put in a single form, which we call the abstract par- 
ity estimate. This is proved at the end of this section through a simple case analysis. 

In the final two sections of the paper, we prove error estimates for the com- 
mutators of the field equations l|1.17(l with certain Lie derivatives. With the work 
developed in the previous sections, this turns out to be quite easy because every- 
thing is just a special case of the abstract parity estimate of section seven. 

Finally, in an appendix we prove certain Sobolev type estimates which are used 
in the main work. This material is included for the sake of completeness and can 
simply be referred to by the reader who wants more detailed account of certain 
(standard) calculations we use. 

Acknowledgments: The authors would like to thank Sergiu Klainerman and Igor 
Rodnianski who played a large role in initiating and supporting this project. This 
investigation originated from a conversation with those two as to whether one could 
extend the low regularity approach of 9 to discuss the scattering behavior of the 
CSF system <|1.17[1 . It was S. Klainerman's insistence that one should first "do it 
classically" in order to get a feel for the bulk behavior of the field quantities, and to 
see if anything interesting turns up, that lead to our study. The possibility that one 
could use space-time estimates to eliminate logarithmic divergences coming from 
the charge was first pointed out to us by I. Rodnianski. This plays a crucial role 
in our work. We would also like to thank him for mentioning the connection with 
conformal geometry that a previous proof of the divergence identities (|4.9|) - l|4.10|) 
had. We have used this elegant description in Section 0] 
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Study, Princeton. H.L. was supported by the NSF grant DMS-0111298 to the 
Institute and partially supported by the the NSF Grant DMS-0200226. J.S. was 
supported by a Veblin Research Instructorship and an NSF postoctoral fellowship. 



2. Some Geometric Preliminaries 

In this section, we introduce the basic geometric concepts and identities which 
we will use in our analysis of the field equations (|1.17(1 . All of this material is com- 
pletely standard, and we review it here in detail solely for the convenience of the 
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reader. Experienced readers will find it worthwhile to skip this section and simply 
refer to it at places in the sequel where a specific calculation requires one of the 
identities we list here. 



2.1. Lie Derivatives. We begin with some well known formulas involving Lie 
derivatives and divergences of vector-fields and differential forms. For any sections 
X, Y to TM we have the usual formula CxY = [X,Y]. We also denote by Cx^ 
and CxF the Lie derivative of a one- form and two- form respectively. In the frame 
{e a } these can be computed as: 

(2.1a) (Cxoj) a = X(uj a )-u([X,e a }) , 

(2.1b) (L x F) a() = X{F aP ) - F([X, e a ],e ) - F{e a , [X, ep]) . 

For vector-fields X, we form the Lorentzian divergence: 

(2.2) (div)X = V a X a = (trace)VX . 

Also, for each vector-field X on X, we measure the effect of its flow on the 
Minkowski metric by forming its Lorentzian deformation tensor: 

(2.3) C x g = {x) tt . 

By Lie differentiating the contraction g a p = g lS g^ a gs l3 we have that: 

(2.4) (W) a/S = - m ^ , 

where we have set g^ = g a @ . Using this formula, one can easily compute the 
commutator of the Lie derivative Cx and the duality operator * applied to a two- 
form: 

[Cx,*]F a(i = \c x ^ a p~ t sg^g Sy )F^ , 

(2.5) = - e a0y s (x) 7T^F aS + ((div)X) *F aP . 
Now, a well known calculation shows that we have: 

(2.6) (x) tt = 2(symm)Vl , 

where (symm) denotes the symmetric part of the tensor VX. Since the trace of an 
antisymmetric tensor is always zero, using (I2.2|l this gives the following formula: 

(2.7) (div)X = - (trace) {x \ . 
A direct use of this gives the calculation: 

A((div)r)-y((div)X) = \c x {g a ^Z Y g) a p)-\c Y {g a \C x g) a0 ) , 

= \g al3 (^ix,Y]g)ai3 , 

(2.8) = (div)(£ x Y) . 

In particular, for any vector-field X such that (div)X = const., we have the follow- 
ing useful formula: 

X((div)Y) = (div)(£jrY) . 
This discussion goes through almost verbatim for one- forms w. In this case we 



STABILITY FOR CSF ON MINKOWSKI SPACE 



13 



define: 

(div)w = V a u a . 

Then using the identification u a — g a °ujp, we have from i|2.8[l and 12.4[1 that in any 
frame {e a } that: 

(2.9) X((div)w) -uj a e a ({div)X) = - V + (div)(£ x ^) , 

Thus, again assuming that (div)X — const, one can drop the second term on the 
left hand side of (|2.9(l above. 

To extend the formulas (|2.8J) and (|2.9|l to higher order tensors, it is useful to 
define the Lie derivative of the Levi-Civita connection V via the formula: 

(2.10) £ X (V) = CxV-VCx ■ 

It is easy to see that (|2.10ll in fact defines the tensor, which in our zero curvature 
setting is given by: 

C X (V) = VVI . 

This object acts in a multilinear fashion on sections to TfM.. For example, with 
respect to a contravariant two-tensor the action reads: 

{C x (V)T) api = {V a W X 5 )T &1 + (V a V 7 X 5 )T^ . 

The above formula immediately applies to two-forms. In this case, if ( x ^tt is a 
constant multiple of the metric g a p, which also implies that (div)X = const, by 
formula (|2.7|l , then antisymmetry and some simple manipulations yields the follow- 
ing important commutator formula: 

(2.11) C x V p F aP = -^tt^V^F^ + V^CxF)^ . 

This will be used in to differentiate the Maxwell equation (|1.8J) with respect to 
various vector-fields which are Killing and conformal Killing, and will be the basis 
of our analysis for that portion of the system l|1.17fl • 

We would now like to set up analogs of some of the above formulas for vector- 
fields and one-forms which are naturally associated with the complex line bundle 
V. We consider complexified versions of the tangent and cotangent bundles: 

(2.12) T C M = V®TM, T* C M = V®T*M. 

In local coordinates, sections to these bundles can simply be identified with complex 
valued vector-fields and one-forms respectively. Of special significance is the full 
covariant derivative of a section to V, which we denote by D(f>. In accordance with 
the above notation, this can also be seen as the complex exterior derivative of the 
scalar <f> which is defined via the natural formula: 

Dcj>(X) = D x <t> ■ 

This should be understood in analogy with the natural formula for the exterior 
derivative of a real valued scalar / on M. which is given by df(X) = X(f). For 
real scalars, one has the Lie derivative formula: Cxdf = dX(f), and we seek an 
analog of this relation in the complex setting. Therefore, we define a covariant Lie 
derivative, called C c on the bundles l|2.12() via the following formulas which hold 
for arbitrary tensors (in analogy with (|1.4|) 1: 

(2.13) C X (T) = CxT+V^l A(X)T , 
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where CxT is somewhat abusive notation in that we consider T as a complex valued 
tensor, whose coefficients are computed in the trivialization of {A Q }. It should be 
noted that by the Leibnitz rule, the formula l|2.13[l is completely gauge independent. 
In the sequel, we will not make explicit use of (|2.13(l which is instead provided as 
a concrete way of visualizing C c . We note here that alternatively, one could define 
C c abstractly via the Leibnitz rule which would be sufficient for our purposes. 

Now, using the commutator formula (|1.5f) . it is seen that the complex Lie deriv- 
ative satisfies the following natural formula with respect to D: 

(2.14) L X D4> = D(D X 4>) + V^lixF-cj). 
Here (ixF) a = X^F^a is the interior product of X and F. 

Next, using formula (|2.13() we have the following useful expression involving 
complex scalars (f> and complex valued one-forms rj: 

(2.15) = 9 {Dx4> ■ rj) + $* (<£ • • 

Furthermore, notice that from (|2.13l) the operation of complex Lie differentiation is 
well behaved with respect to contractions. For example, for a real valued two-from 
and real and complex vector-fields X, Y and Z (resp.) we have the formula: 

(2.16) D Y F(X,Z) = C Y F(X,Z) + F([Y,X],Z) + F(X,£$Z) . 

We will use this formula instead of the corresponding formula involving the covari- 
ant derivative Vy. This is necessary because when we differentiate F it will always 
be with respect to C and not V. 

Finally, using a more abstract form of the formula (|2.13|) . we compute the co- 
variant divergence of a complex valued one-form via the implicit formula: 

(div) c (0®w) = uj a D a( j) + (div)w • <f> , 

for any complex scalar field <j> and real one-form uj. Using this last line in conjunction 
with the identity l|2.9|l above, we have the following commutator formula for the 
complex divergence and complex Lie derivative for complex valued one-forms r\: 

D x {{^fn)~n a e a (i^)X) = - Dj- X \^ m + (divf(C x r,) + V=lX a F^ . 

In particular, if ^ x ^tt is a constant multiple of the metric g a /3, we have the follow- 
ing formula for the commutator of the covariant wave equation and the covariant 
derivative (or Lie derivative if you like) Dx' 

(2.17) D a D a D x 4> ! 

= (div) C D(D x cl>) , 

= (divf (cScDc/)- V^lixF-4>) , 

= D x D a D a <f> + ( x \ a PD a Dp(t> + v^T {2X a F aP D^ - V Q (^F Q/3 ) • <j>) . 

The important thing to notice here is that the above formula depends only on the 
curvature F and not the choice of gauge potentials {A Q }. 
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2.2. The Null Frame, Lorentz Algebra, and Associated Covariant Deriva- 
tives. Of central importance to the analysis we do here is the freedom to perform 
calculations on the equations i|1.17[l in an arbitrary frame. In particular, the ability 
to use frames which do not arise from any system of coordinates. Since we are 
assuming that our initial data is both smooth and well localized around the origin 
in R 3 , there is a canonical choice of frame for our problem. This is the so called 
standard spherical null frame. The first two members of this are the null generators 
of forward and backward light cones which we define respectively as: 

(2.18) L = d t +d r , L = d t - d r . 

To complete things, we need only define derivatives in the angular directions. This 
can be done in an identical fashion on each time slice {t = const.} so we only need 
to define things on K 3 . If we let {e^}x=i,2 denote a local orthonormal frame for 
the unit sphere in R 3 , then for each value of the radial variable we can by extension 
define: 

(2.19) e A = -e° A . 

r 

Thus the collection {e A } a=i.2 forms an orthonormal basis on each sphere {r = 
const.} on each fixed time slice. We can write the usual translation invariant frame 
{di} in terms of d r and this basis as follows: 

di = uj l d r + ujfe A , 
where uif = e A (x l ), which follows at once from the formula: 

(2.20) cuf = (d h e A ) = dx\e A ) = e A {x l ) . 

In line above we have used in order orthonormality, duality, and then simply the 
definition of exterior derivative. Also, we note here that the {uf} are part of a 
rotation matrix at each point, which implies the following useful formulas: 

(2.21) u^t = , u ' A u>? = 5 A B , = 8/ - w iU i . 



With respect to the full frame {L, L, e A }, the Minkowski metric reads: 

g = -26 L ® 9 L - 26 L ®6 L + 6 AB 9 A ® 9 B , 

where {9—, 9 L ,9 A } is the corresponding dual frame. Associated with this null frame 
are the standard optical functions: 

(2.22) u = t + r , u = t-r . 

Note that one has the identities: 

L(u) = 2 , L{u) = , 

L{u) = , L(u) = 2 . 

In particular both u and u solve the eikonal equation V Q /iV ' a h — 0. In terms of the 
frame (|2.18|I - H2.19(I we have the standard null decomposition of the electro-magnetic 
field tensor F, as well as the gradient D<j>. As usual we define: 

(2.23a) a A = F LA , a A = F LA , 

(2.23b) p = ^F L l , o = i e AB F AB ■ 
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Using the duality operator * , we have the useful formulas for the null decomposition 
of *F, where €ab= \ ^llab is the volume element on the spheres r = const.: 

(2.24a) *c<a = - a B , *QLa = ^A B SLb s 

(2.24b) *p = a , V = -p . 

One can expand out the energy-momentum tensors {EE} and in the null 

directions l|2.18|l . Using the identifications l|2.23|l . for the electro- magnetic field 
these read: 

(2.25a) Q(L,L)[F] = 5 AB a A a B = \a\ 2 , 

(2.25b) Q(L,L)[F] = 5 AB a A a B = \a\ 2 , 

(2.25c) Q(L,L)[F] = p 2 + a 2 . 

Furthermore, for the energy-momentum tensor of the scalar field (|1.13|) . one has 
that: 

(2.26a) Q{L,L)[4>] = \D L ^\ 2 , 

(2.26b) Q{L,L)[<t>] = \DiA\ 2 , 



(2.26c) Q{L,L)[(j)] = 5 AB D A (j)D B (j) = \p<f>\ 



Next, we record here the following standard formulas for the Levi-Civita con- 
nection V with respect to this frame. This will be used many times in the sequel: 



(2.27a) 


V L i 


= o, 






V L L = 


, 


(2.27b) 




= o , 






VlL = 


o, 


(2.27c) 




= o , 






^L_eA = 


, 


(2.27d) 


X7 eA L 


1 

= -e-A , 
r 








1 

— eA 
r 


(2.27e) 




= V eA e B -) 


-^ab(L- 


L) , 






(2.27f) 




= 1 AB e D - 


Itt5ab{L- 
Zr 


■L) . 







Here, V denotes the intrinsic covariant differentiation on spheres, and T its associ- 
ated frame- Christoffel symbols. The key property of the T we will use in the sequel 
is that they are homogeneous functions of degree —1 with respect to the radial 
variable r. 

Our result is based in part on the usual process of obtaining weighted energy 
inequalities, which in turn give L°° estimates of the type (|1.26JI - dl.27|) . The weights 
in these kind of estimates are closely related to the Killing and conformal-Killing 
structure of Minkowski space, as was first realized to full effect in the seminal work 
[Hj. Accordingly, we introduce the inhomogeneous Lorentz algebra: 



(2.28) 



L — {d a , fl a p, S} , 
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where we use x° — t and xq = —t to express: 

(2.29) Q Qj g = X a df3 - Xftd a , 

(2.30) S = x a d a . 

This collection of vector-fields satisfy the well known algebraic relations: 
(2.31a) [d a ,VL M ] = 5 a{ pd^ , 

(2.31b) [n aPl tt ia ] = -8 {ai np a) , 

(2.31c) [d a ,S] = d a , 

(2.31d) p a p,S\ = 0. 

Expanding the formulas H2.29fKI|2.3l)l) out m tne frame l|2.18|l - H2.19|) . and using the 
symbol ujf — 9 A (di) we have the null decompositions: 



(2.32a) 


d = 


\{L + L) , 


(2.32b) 


d t = 


^(L-L)+u;fe A , 


(2.32c) 


n« = 


flfjeA = {xiLof ~ XjUjf 


(2.32d) 




y (uL - uL) + iwf e A , 


(2.32e) 


5 = 


- (mL + uL) . 



It will also be necessary for us to have a precise account of the algebraic relations 
between the frame <|2.18(l - i|2.19|l and the fields l|2.28|l . These can easily be computed 
using the formulas l|2.32|l . We record this computation here as: 

[L.,do] = [e A ,d } = , 



(2.33a) 


[L,do] 


(2.33b) 


[LA] 


(2.33c) 




(2.33d) 


[L,%] 


(2.33e) 




(2.33f) 


[L,Slm] 


(2.33g) 


\L, fiio] 


(2.33h) 


[eB,£ko] 


(2.33i) 


[L,S] 


(2.33j) 


[L,S] 


(2.33k) 


[e A ,S] 



u}f[L,e A ] = ~ — e A 
r 

, A 

[L,e A ] = -*-e A , 
r 



A 



= -\L,SUj] = , 



u 



ujiL u\ e A , 

r 

a 

—WiL + —Mi e A , 
r 

— {uL - uL) - tuj % T AB e D , 
2r 

L , 

L, 

e A . 

In formula (|2.33e|) above we have set c AB = [e^es]- , which is again a smooth 
function homogeneous of degree —1 in the radial variable r. Also, on line (|2.33h[) 
we have used the identities: 

e B (uf) = V eB 9 A {di) = -Si^-TBD"? > Cba = T BA -T AB . 
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Finally, to end this subsection, we list here the various formula for the action 
of the Levi-Civita connection V on the algebra (j2.28(l . All of these formulas are 
computed in a straightforward way using the identities (|2~2Tjl and J232): 



(2.34a) 




= Vl9 q = V eA <9 Q = , 


(2.34b) 




= — Vlflij = -Oy , 


(2.34c) 






(2.34d) 




= Lu t L + LufeA , 


(2.34e) 




= -UiL + ujfe A , 


(2.34f) 




= ^f(L + L), 


(2.34g) 


v L s 


= L , 


(2.34h) 




= L, 


(2.34i) 


V eA S 


= e A ■ 



3. Fixed Time and Space-time Energy Estimates for the Curvature 

In this section, we begin to build the estimates which lie at the center of our 
approach. All of these will be produced in some way through the tensorial conser- 
vation laws (|1.12(l and (|1.14() . In this section we will deal with weighted L 2 type 
estimates which involve the curvature tensor F a p . We assume that this satisfies the 
Maxwell equation (|1.8|l for an unspecified current density J a . 

On Minkowski space, the two most basic energy estimates (really identities) 
for F a p are based respectively on time translation invariance and the conformal 
structure of the Minkowski metric. This is utilized by contracting the energy- 
momentum tensor (|1 . 1 l|l with the following two vector fields (respectively): 

(3.1) T = d = 1 -{L + L) , 

(3.2) Ko = (t 2 + \x\ 2 )d + 2tx l d l = ^(u 2 L + u 2 L). 
One readily computes their deformation tensors to be: 

(t) tt = , {Ko} Tr = Ug . 

Due to the trace- free nature of the energy-momentum tensor i|l.llf) . the one- form 
resulting from such a contraction is seen to be divergence free. Therefore, defining 
the weights: 

(3.3) t\ = 1 + u 2 , t 2 _ = 1 + u 2 
and the hybrid vector-field: 

Ko = T + K Q , 
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we arrive at the following fundamental estimate through the application of the 
geometric Stokes theorem 8 over domains of the type {u ^ uo} H {0 ^ t ^ to} and 
expanding out the resulting quantities using the identities <|2.25[) : 

(3.4) f T 2 + \a\ 2 + T 2 \a\ 2 + (r 2 + +r 2 _)(p 2 +a 2 )dx 
jR3n{t=t } 

+ sup / r 2 + \a\ 2 + t 2 _{ P 2 + a 2 ) dV c(u) 

« JC(«)n{o<t<t } 

< r I \{KoYF a0 J p \dxdt + f (l + r 2 )(\E\ 2 + \H\ 2 ) dx . 

Jo Jm. 3 JR3n{t=o} 

Here dVc( u ) is the Euclidean volume element on the cone {t — r = u}. 

As it stands, estimate (13 .4|) is not terribly useful except in certain restricted 
situations. This is because in general it is not possible to assume that the right 
hand side of l|3.4|l above is finite, even if one assumes that the "dynamic" portion 
of the initial data is compactly supported. This is the effect of the constraint 
equation on line (|l.lUa|l . Indeed, taking a Hodge decomposition E = E d ^ + E c ? of 
E into its divergence- free and curl- free components (respectively), and introducing 
the potential function: 

E cf = Vip , 

we see that this constraint is equivalent to the elliptic equation: 

(3.5) Aip = Jo . 

In particular, it is not possible for E to decay better that E ~ -K unless the 
following quantity known as the charge vanishes: 

(3.6) q{F)(t) = / J (t) dx . 

Notice that this quantity is a constant of motion thanks to the continuity equation 
(|1.16|) . Therefore, we shall henceforth refer to it as q(F) and drop the dependence 
on t. Because of the weights involved, there is no hope of obtaining a finite value 
for the energy (|3.4|l without the extra assumption that q(F) = 0. 

In order to circumvent this problem it is necessary to either modify the energy 
(|3.4|) so that less decay of the initial data is required, or to modify the field equa- 
tions p.8|l themselves so that the behavior resulting from the charge is eliminated. 
A naive approach to the first tactic would be to simply decrease the amount of 
decay required on the right hand side of l|3.4|l . For example, this can be done by 
eliminating the use of the conformal field (|3.2|l . As we shall see in a moment, it 
does not seem possible to do this and still retain the distinct weights (peeling) on 
the different components on the left hand side of (|3.4|l . Furthermore, recall that 
the simple scaling argument used in the introduction implies that it should not be 
possible to close a global existence proof for the system (|1.17|) under the assump- 
tion of a total decay rate for the curvature less than F ~ ^ at space-like infinity. 
Therefore, we use a different approach here. 



8 We shall discuss this type of procedure in more detail in the sequel. See also for a more 
thorough discussion and a derivation of this particular estimate. 
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A more sophisticated approach to the first tactic, which has been used by other 
authors in similar contexts 0], would be to employ some Lie derivatives to the 
field strength F a p before putting it in weighted L 2 in order to kill off the effect 
of the charge. For instance, the Lie derivatives with respect to the rotation fields 
Qij work well to accomplish this because the leading order behavior of the charge 
is spherically symmetric. However, this strategy runs into technical complications 
when it comes to dealing with the boost vector fields Qio, which we will make 
essential use of in our proof of the sharp weighted L 2 (L°°). Therefore, we shall 
follow a tactic which is both conceptually and technically much simpler. This is as 
follows: Solving the charge equation l|3.5() explicitly, we see that: 



If 1 

ip(t,x) =-— -J (t,y)dy. 

4tt J r3 \x - y\ 

Therefore, at fixed time, we have the following asymptotic behavior for E c ? as 
r — > oo: 

Since this behavior is sufficiently simple we shall subtract it off, proving estimates 
for the remaining field strength. However, to do this correctly it is necessary to 
take into account the behavior of the charge in the wave zone t *~ r. Although it 
is not obvious at first, there is a jump type behavior of the charged component of 
E across the cone u = as long as the current vector J a satisfies certain weighted 
estimates which are consistent with the right hand side of H3.4[l above. Keeping 
(|3.7(l in mind, the approximation we make to the charged component E° is the 
following. We first define the charge two-form, denoted by -Fa/3, via the exterior 
derivative (for ^ <): 

(3.8) F = - d ([J^ l x +( a - 1 - 2) da] dt 

In the above formula, y + denotes a smoothed out version of the Heaviside function; 
i.e. a non-decreasing C°° function equal to on (— oo,0] and 1 on [1, oo). Notice 
that by design, the tensor F a p satisfies the Bianchi identity (|I.8bJ) because it is 
exact. Also, one has the explicit formulas for the electric-magnetic decomposition 
(Ol) of F: 

(3.9) E % = q .^ x +( r -t-2) , Hi = 0. 

This last expression shows that F indeed represents a charge with asymptotic (j3.7(l 
which propagates in the exterior of the light-cone u — 0. In terms of the null 
decomposition l|2.23|l . the formula l|3.9() is even simpler and becomes: 

(3.10) p = p(F) = oj^Ei = q- -^x + (r-t~2) , a = a = <j = 0, 
and the associated current vector 9 is easily computed to be: 

Jk = Q-^(x + )'(r-t-2) , J L = 7.4 = 0, 



9 That is, plugging F into the field equation 11.81 . 
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In particular, we have that: 

(3.11) 7 (t,r) = q ■ ( X + Y(r - t - 2) . 

Notice that the current vector J a is only supported in the region u ~ and has 
only a L component. This will allow us to obtain good estimates when it appears 
as an error on the right hand side of expressions such as (|3.4J) . 

The main estimates we need to prove now are for the remaining field strength: 

(3.12) F = F-F . 

This satisfies the field equations (|1.8|) with current vector: 

(3.13) J a = J a J a. • 

A simple calculation using the formula (|3.11(1 shows that the charge for the modified 
field strength F vanishes: 

q(F) = J (Jo -Jo) dx = 0. 

This formula, used in conjunction with certain weighted elliptic estimates for gra- 
dients (see appendix 1 10| l. allows control of the right hand side of l|3.4|l above with 
F replaced by F. We shall return to this in a moment, after we have first discussed 
the previously mentioned idea of decreasing the powers occurring in the weights in 
estimate (13.41) . 

In the sequel we shall consider Maxwell fields F a p such that the remainder F a p 
satisfies various rates of decay; not necessarily the amount which comes from the 
usual Morawetz estimate (I3.4|l . This can be understood as making different as- 
sumptions as to how "asymptotically flat" the initial data for F a p is after one has 
subtracted off the contribution of the charge 1)3. 9f) at time t = 0. The least amount 
of decay which is consistent with a finite value for the charge (13. 61) is essentially 
E ~ r ^( 2+ T) for some < 7. Large values for 7 are available depending on the 
weighted decay of Jo as well as the vanishing of higher moments for this quantity. 
In order to be able to take these various rates of decay into account for the range 
< 7 ^ 1, we would like a version of estimate (|3.4Jl with weights of the form 
(1 + r 2 ) s on the right hand side for i ^ s < 1. It will be essential for us to be able 
to do this in such a way that distinct decay rates, also known as peeling, are still 
available for the various components in the null decomposition (|2.23|) . This will be 
accomplished below through a suitable modification of the Morawetz field l|3.2fl . 



Finally, before stating and proving the main result of this section, let us mention 
in words one further adjustment to the classical Morawetz estimate (|3.4[1 which will 
be of great use in the sequel. This involves adding some space-time energies to 
the left hand side of l|3.4fl . These additional estimates turn out to be more than 
just a mere technical convenience, and are in fact essential for us to avoid certain 
logarithmic divergences which would enter when trying to use the usual procedure 
of matching fixed time energy and pure L°° bounds. 
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We now state the basic general energy estimate we will use for the curvature F a $ 
which is: 



Proposition 3.1 (Generalized Morawetz estimate for the electro-magnetic field 
with non- vanishing charge). Let F a $ be a two- form which satisfies the system (|1.8(l 
with current vector J a , and let q(F) and F be the associated charge and remainder 
field strength defined by formulas (|3.6() and (|3.12l) respectively. Furthermore, let 
< 7, e and ^ < s $C 1 be given parameters such that s + j< | as well as e s — h, 
and define the weights: 



(3.14) 


TO 


= T-/T+ , 






(3.15) 


u? 7 (t,r) 


= X + {r-t)-T 2 J 4 


(1- 


-X + {r-t)) , 


(3.16) 


w 7te (t,r) 


= X + (r-t)-T 2 J + 


(1- 


-X + (r-t))-r 2 J , 


(3.17) 


w!y, e (t,r) 


= x + ( r ~ t) ■ T - 7 1 


+ (1 


-X + (r-t))-T-^- 1 


Define the 


remainder 


null decomposition: 






(3.18a) 










OLA 


= F LA = 


Fla , 


QLa 


= Fla = Fla , 


(3.18b) p 


= ^{F LL 


— 1 ~ 

- Fll) — t^Fll , 


a 


- 2 E FAB 2 



-AB p 
- £ AB 



Now define the charge modified generalized Morawetz type energy content of F in 
the time slab {0 < t s$ t } n R 3 to be: 

(3.19) E^(0,t )[F] = \ q (F)\ 2 

+ sup / (Tl s \a\ 2 +T 2s \a\ 2 + Tl s (p 2 + a 2 )) w 7 dx 

+ sup f (rl s \a\ 2 + r 2s (p 2 + a 2 )) w, dV c(u) 

« J {C(tt)}n{0<i<t o } 

+ r I (*+ M 2 + T l +2t (^l«| 2 + T + P 2 + T +° 2 )) <e dxdt . 
Jo Jft 3 

Then one has the following general weighted energy type estimate: 

(3.20) E^' e \o,t Q )[F] < 

f° I {t 2 + s+1+2 *tZ 2 Vl\ 2 + rl +2 ^ 2s rt s ~ 2 Vif + 7? T-\ /| 2 ) tfl 7)t dxdt 
Jo JR 3 

(l + r 2 y+~< (\E d f\ 2 + \H\ 2 ) dx +||(l + r) s +TJ (0)|| 2 f 



3 n{t=o} 

Here, C(u) denotes the forward facing light-cone t — r — u for fixed values of u. 
Furthermore, \fi\ 2 — S ab JaJb denotes the angular portion of \J\ 2 . 



Proof of estimate (|3.20|l . The proof closely follows the general strategy for proving 
estimates of the type l|3.4|l with a few additions. We first introduce a warped 
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generalization of the conformal killing field (|3.2|) . This is: 

(3.21) K° = \u 2s L+ l -\u\^L. 

Here we will allow s to range as J ^ s ^ 1. In order to gain some intuition as to 
the nature of Kq , notice that this vector-field essentially interpolates between the 
conformal Killing fields K and S, which are its endpoint values in the interior re- 
gion u for s — 1 and s — | respectively. In order for this analogy to have value, 
it must be seen that the deformation tensor of Kq enjoys some positivity property 
when contracted with a trace-free 2-tensor which satisfies the positive energy con- 
dition. This is indeed seen to be the case through the following calculations. We 
let: 

v(u) = u 2s , v(u) = \u\ 2s , 

be short hand for the fractional weights. Then, with respect to the Minkowski 
metric g a p, the deformation tensor of Kq is easily calculated using the identities 
lHUl and $FH\ : 

(3.22) (K "\ = ^—^ g + 2(^—^-(v + v))(e^<g)e L + e L (g>9 L ). 



In the above expression, the dot notation is used to indicate the derivative of v and 
v as functions of a single variable. For example, we have that v = 2s|u| 2s ~ 1 -sgn(u). 
The positivity claim now follows once we have shown that the factor on the second 
term in H3.22[l is non-negative: 

Lemma 3.2. Setting let \ ^ s ^ 1 and v = u 2s and v — u 2s . Then one has that: 
(3.23) ^ (^— --(v + v) 



Proof of (|3.23|l . This follows immediately from freshman calculus. For fixed t, we 
define the function of r: 

f(r) = v(t + r) — v(t — r) . 

Notice that we have /(0) = 0. We also have that f = v + v. Therefore, from the 
mean value theorem there exists some rg G [0,r] such that: 

^ = /'M = v{t + r ) + v(t-r ) . 
r 

The claim (|3.23|) then follows if we know that /' is a non-increasing function of r. 
This function is computed to be: 

(3.24) f'(r) = 2s{r + t) 2s - 1 + 2s\t - r\ 2s ^ ■ sgn(t - r) . 

For s = \ , this function is easily seen to be decreasing by direct inspection. For the 
range | < s ^ 1, (|3.24l) is continuous so we only need to show that /" ^ when 
< r. Calculating the second derivative we see that: 



f"(r) = (2s-l)2s(u 



2-2 _ i ,2a-2 



The desired result now follows from the fact that < (2s — l)2s while (2s — 2) ^ 
and \u\ < u. □ 
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We now form a momentum density associated with the charge free portion F of 
F, as according to equations lj3.12|l and QH.8JI . and the extra 7, e weights defined in 
the statement of Proposition 13 . II above: 

(3-25) Pi s '^[F} = Q af3 [F]{K s f-w^, 

where K 8 = T + Kq and the weight function {y 7i< = is defined as: 

(3.26) fi; 7ie = ({l + (2-u)^) X +(~u) + (l + (2 + u)-^)(l- X +(-u))) + 

(l + u)-^-{(2~u) 2 ^ X +(-u) + l- X + (-u)) • 

Here \ + can De taken to be the same smoothed out Heaviside function used in line 
(13.81) above. In particular, (x + )' is positive and supported on the interval [0,1]. 
This last assumption assures that the term (2 + u)~ 2e (l — X + {~ u )) on the right 
hand side above is never singular. The weight function uL e and its derivatives 
satisfy simple bounds with respect to the weight functions u> 7 and w!y e defined on 
lines 13.15(I - H3.17|I above. First, notice that one has the bounds: 

for a suitable positive constant C . We will restate this estimate as follows: 

(3.27) w 7 ^ e ~ w 7 . 
Furthermore, a brief calculation shows that: 

= 2e(l + u)-^- 1 - {[2-u) 2 ^x + {-u) + l-X + (-u)) , 

(3.28) ~ C-l r 1+2 ^; >e . 
A similar calculation shows that: 

--^L{wj, € ) = -d u (w lt e) , 

= ((2-ur-(2 + u)-^)(x + y(~u) + 
{l + u)-^({2-uf-<+^-l)(x + )'{-u) + 
2( 7 + e)(l + u)- 2t {2 - u^+^x+i-u) + 
2 7 (2-u)^- 1 X + (-u) + 2e(2 + u)- 2f - 1 (l-x + (-u)) , 

(3.29) ~ C~l w'^ e . 

The bound on the last line follows because one has: 

(2 + u) s > < (2 - u) 5 * , 

whenever 8\ ^ 82 and u G [—1,0]. Note again that {x + Y is a positive function. 
Also, it is clear that the constant implicit in the ~ notation on lines H3.27|) - l|3.29|) 
above can be taken to be the same value uniform in e and 7 because we kept the 
effect of these constants in our bounds with the G 1:<L notation. We will continue to 
use the ~ notation in the remainder of this proof with the understanding that the 
same implicit constant is being used throughout. 
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Now, calculating the space-time divergence of the quantity (|3.25|) we have that: 

(3.30) V a P^^[F] = F aP J a {K s f .«J 7 , e + ±Q af3 [F] ^n"? ■ w 1<£ 

\Ql P [F]{kIY ■ L{w 1>€ ) - \Q L p{F](KlY-L{w 1>e ). 

Integrating both sides of this last line over time slabs of the form: 

K(t ,u ) = {0 ^ t < t } n {u < u } , 

and applying the geometric version of the Stokes theorem (Gauss theorem) we arrive 
at the following general integral identity: 

(3.31) f P {sn ' e) [F] dx = [ P^' e) [F] dx 

+ / p^ 7 ' c) [f] dv c(ua) + [ [ (R.H.S.)G30 • 

Jc(u )n{o^t^t } J Jii(t ,u ) 

Notice that the above identity reduces to the usual energy type estimate on the 
time slab {0 ^ t ^ to} H M 3 when to ^ uq- In order to proceed, we now calculate 
each of the terms in (|3.31|) individually. 

When t = 0, we compute that: 

(3.32) P (s ' 7 ' e) [P] ~ (l + r 2 ) s+ ->{\E\ 2 + \H\ 2 ) , 

where again the ~ notation here means that the ratio of the two terms above is 
bounded by on the left and right by C _1 and C respectively for a suitable positive 
constant C . In particular, if the right hand side of such an expression is positive, 
then so is the left. Also, E denotes the electric part of the tensor F while H is the 
magnetic part of F. Note that: 

E = E-E , 

where E is defined as on line H3.9[) above. Likewise, with the help of the identities 
(|2.25|l . we compute that at time t = to'. 

(3.33) P^[F] ~ (^|a| a + ^|a| 2 +^V + ^))-«' 7 . 
We also compute the characteristic energy term: 

(3.34) P^[F] ~ (rl*\a\ 2 +T 2 J{? + o- 2 )) -w 7 . 

It remains to calculate the terms on the right hand side of (|3.3UI) above. Since the 
fist such term does not have a sign, we simply expand it using the null decomposition 
(|3.18a|l -l )3.18bj) and the definition of J given by H3.13|l above as well as the estimate 
(|3~27|l to bound: 

(3.35) \F a pJ a {K S of -w, t A < 

(t^(H.|/| + |p|.|J l |) + r! s (|a|.|/| + |p|.|J^|)).u; 7 + \p\ ■ \q\ ■ t^tZ 1 <V 7 . 

Notice that the extra t_ weight on the last term of the right hand side above comes 
because J is supported in the region where t_ ~ 1. 
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We now calculate the second term on the right hand side of 1(3.30(1 above. Because 
Qap[F] is trace- free and T is Killing with respect to gap, the contraction need only 
be taken with the second term on the right hand side of line ((3.22|) . This yields: 

(3.36) Q aP [F] (^0)^/3 . ^ „ f^ZJL _ + y^j . (p3 + ct2) . ^ _ 

In particular, using Lemma (3.21 and the range restriction i ^ s 1 we have that 
this last expression is non-negative. 

Moving on, we calculate the third term on the right hand side of ((3.30(1 . Using 
the bound 1)3.28(1 and the decomposition 1(2.25(1 we conclude that: 

(3.37) - \Q^[F](K S ,f • L(w lie ) ~ C~l r 1+2c (r 2 J\a\ 2 + ^{p 2 + a 2 )) ■ < e . 
Similarly, using the bound 1(3.29(1 we estimate: 

(3.38) - \q L!} [F]{kIY ■ £K, £ ) ~ C-l (r 2s \a\ 2 + r 2 V + a 2 )) • < £ . 



Therefore, collecting the positive terms 1(3.33(1 - 1(3.34(1 and 1(3.36(1 - 1(3.38(1 onto the 
right hand side of ((3.31|) above, and collecting 1(3.35(1 with ((3.32|) on the left hand 
side and using the shorthand notation ((3.46(1 , which will be introduced in a moment, 
we can estimate: 



(3.39) 



{t=0}nR 3 
to 



(l+r 2 y+^(\E\ 2 + \H\ 2 ) dx + 



JR 3 



\p\ ' M ' T + 2 T_ w w-y dxdt 



(rl s (\a\ ■ | /| + \p\ ■ \J L \) + r 2 J(\a\ ■ | /| + \p\ ■ \J L \)) ■ w, 



dxdt 



> C. 



7>« 



E^^(0,t o )[F] + 



JR 3 



v — V 



— (£> + «))• (p 2 + a 2 )-w 1 dxdt 



We now make a preliminary reduction on estimate 1(3. 39(1 above by first discarding 
the second positive term on the right hand side, and then using a Cauchy-Schwartz 
to peel off a factor of (E^y s,ry ' e '(0,to)[F] from the latter two terms on the left hand 
side of the above expression. Setting: 



(3.40) |||J 



lL 2 [0,i ](L 2 )( S , 7 ,e) - 

V_|J £ | a 



2s-l-2e_2s+l I 



t1 ~\Jl\ + T+T-lftf) w 7 , e dxdt , 

by a simple use of Holders inequality, the condition e ^ s — | , and the weight bound 
w 2 ^ T—w' e u> 7i£ we may replace ((3.39(1 with the estimate: 

E^(0,t )[F] < C 7 4(^^) (s ' 7 ' e) (0,to)[i ? ]-(|||J|||L 2 [o,to](L 2 )(,, 7 , e ) + |g|) 

(l + r 2 y+~<(\E\ 2 + \H\ 2 ) dx 

{t=o}nR 3 



STABILITY FOR CSF ON MINKOWSKI SPACE 



27 



from which easily follows the bound: 
(3.41) 



2 I 1 2 

L 2 [0,t ](L 2 )(s,7,e) + M 



{*=0}nR 3 



(l+r 2 ) s+ ^(|£:| 2 + |iJ| 2 )da; 



In order to deduce l|3.2U[) from i|3.41fl . we need to prove the following bounds for 
the range ~ < s ^ 1 and < 7, where s + 7 < I: 



(3.42) 

and: 

(3.43) 



\q\ 2 < C T ||(l + r)^Jo(0)|| 



(l + r 2 ) s+ ^\E\ 2 dx 



{t=o}nR 3 



(l+r 2 ) s +^| 2 + ||(l + r) fl+ TJ o (0)||% 
{t=o}nR 3 L l 



The first of the above estimates, Q3.42J1 . follows from a simple application of 
Holders inequality with the weight (1 + r)~ s ~ 7 in the integral (|3.6|l at time t = 0. 
Note that | < s + 7 because of our range restrictions so the resulting factor integral 
converges with a bound depending only on 7. 

The second of the above estimates results from a Hodge decomposition and the 
weighted elliptic estimate (|10.2fl from Appendix 1 1UI Expanding out E on the left 
hand side we can bound: 

{l + r 2 ) s+ ->\E\ 2 dx ^ 2 f (l + r 2 ) 8 ^ (\E df \ 2 + \E C} -E\ 2 ) dx . 

{t=o}nR 3 J{t=o}nR 3 

Introducing a potential function ip for E°* as in line (|3.5(l above and recalling the 
definition 1)3. 9fl of E, and using line (|3.42l) to estimate: 



(1+r 2 ) 



{*=o}nR 3 r 
we are reduced to proving the estimate: 



V xX + (r-2)\ 2 dx < C 7 |g(F)| 2 s=: C 7 || (l+r) s +V o (0) || 



2 



(3.44) 



1 



{t=o}nR 3 



(1 + r 2 ) .s+7| V (_ j o(0) + ^ X + (r 2))\ A dx 



< C 7 ||(l + r ) s +-rj (0)||^ f 



By using the usual L$ c — > L 2 Sobolev embedding, as well as the bounds: 



{t=o}nR 3 



l V (^7X + (r-2))| 2 dx < \q(F)\ 2 < C 7 ||(l + r) s +V o (0)|| 



r 2( S +7) | V (^(l- Y +(r-2)))|^ < C 7 |g(F)| 2 sC C 7 || (1 + r)^ J o (0) || 



'{*= 



47rr 
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the second of which follows from the condition h < s + 7, we are reduced to showing 
that: 

(3.45) / r 2( S +7)| V ( 1 Jo(0) + ^£))|2 dx < c \\r^M0)\\\ . 

J{t=o}nm & ^ hi 

Using the definition (|3.6() of q(F), this is precisely the statement of (|10.2|l in Ap- 
pendix Notice that the condition i<s + 7<|is enforced by assumption. 
This ends the proof of estimate (|3.2U|) . □ 



We conclude this section with an important generalization of Proposition 13. II 
This involves estimates for derivatives of the field strength F a p. As it turns out, 
the simple decomposition (|3.12|) is remarkably robust with respect to the operation 
of Lie differentiation. One can show that Z^F^ satisfies bounds similar to (|3.2U|) 
with the appropriate right hand side. This will be used in the sequel to show that 
the various components of the null decomposition of F Q/ g satisfy the expected L°° 
estimates. Together with the fact that F a p is given explicitly, these bounds will 
fully determine the point-wise behavior of the original field strength F a p. We will 
prove: 

Proposition 3.3 (Generalized Morawetz estimate for derivatives of the electro-magnetic 
field ). Let F a p be a two- form which satisfies the system (|1.8|) with current vector 
J a , and let q(F) and F be the associated charge and remainder field strength de- 
fined by formulas (|3.6[) and i|3.12|l respectively. Furthermore, let < 7, e and 
h < s ^ 1 be given parameters such that s + 7 < | and e ^ s — ^ . Define the 
weights tq, uj 7 , u> 7i£ , u> 7 e as on lines (|3.14|) - (|3.17|l above. Next, define the gener- 
alized (not charge modified) Morawetz type energy content of a two-from F in the 
time slab {0^t^« }nt 3 to be: 

(3.46) 

E^> e \0,t )[F] = sup / (r 2 + s \a\ 2 +T 2s \g\ 2 +T 2 + S {p 2 +a 2 )) Wl dx 




dxdt . 



Now, define the k th weighted charge modified Morawetz energy: 



(3.47) Ei s ^\0,t o )[F] = q 2 (F) + £ (S|7,e) ((Mo)[^F] . 

|/|<fe 
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Then, recalling the vector norm ||| • |||L 2 [o,to](L 2 )(s,7,e) from line H3.40|) above, 
have the following general energy estimate: 



(3.48) Ei s ^\0,t o )[F] < Cl e [ £ ( W\£xJ 



i 2 [0,t ](L 2 )( S , 7 , e ) 

I\^k 



{1+r 2y+i+\i\ (\vlE<Vf + \ViH\ 2 ) dx + ||(l+r) s+ ^l J lv£Jo(0)|| 2 f ) 
1 n{t=o} L S ' 

+ £ ll(i + ^ +7+|/|+ M,^(o)lli 2 " • 

Here V x and Vt,x denote the Lie derivatives of the covector J and the two-form 
F with respect to the spatial translation invariant fields {di} and the full set of 
translation invariant fields {d a } respectively. 



Proof of estimate l|3.48|l . The proof is essentially that of estimate (|3.2U|) applied 
to the field C x F a p with some additional calculations at the end to wrap things 
up. First notice that by (|2.11() and the fact that for any Lie derivative one has 
[d, Cx] = 0, we have the following formula for C X F whenever IeL \ {S}: 

^(C x F) aP = (C x J) a - (C x J) a , 
W P *(C X F) = . 

The effect of the Lie derivative £$ is equally easy to account for. In this case 
formula (|2.11(l gives: 

V {C S F) a p = {J-J) a + {C S J) a - {C S J)a , 

VP*(£ S F) = . 

Using the above formulas and following the proof of i|3.2U[l until line l|3.39|l yields 
the inequality: 

Ei^' e \0,t o )[F] < C 7 , e J2 [l^)| 2 + lll^(J- J)|||i 2[ o^](L 2 )(s , 7 , e ) 

xeh 

+ I (l + ry^(\E(£ x F)\ 2 + \H(C x F)\ 2 ) dx' . 

Using the same steps as at the end of the proof of H3.20(l to bound the quantity 
| (7(^)1, we are done once we have shown that: 

(3-49) III^J|||i 2[0;to](i2)(s , 7>e) < \q(F)\ 2 , 
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as well as the bound: 



(3.50) I (l + r 2 y^(\E(C x F)\ 2 + \H(/: x F)\ 2 ) dx 

£ [ / (l + r 2 r+^l(|V^| 2 + |V^| 2 ) dx 



< 



{t=0} 

(i + ry+'y+WviMo) 

+ E II(i + O s+7+|/|+1 vLj(o) 



\i\<k-i 



We begin with l|3.49|l . This will follow from direct computation using point- wise 
bounds on the quantity C X J for X Eh. These point-wise bounds will be provided 
through induction on the value |/|. Notice that when |/| = we may write: 

(3.51) 

j l = q - — xw , Jl = q- r4 ■ xm , J a = q- —3 — xw , 

where x is a C°° an d 0(1) bump function adapted to the origin u = and the 
Q a (io) are smooth function of the angular variable only. In particular, these satisfy 
inductive identities: 

L(Q a ) = L(fi a ) = 0, e A (Q a ) = -Cl a . 

r 

where the £l a have the same properties of the £l a . Also, notice that by direct compu- 
tation one can substitute the decay rates (|3.51|l into the norm ||| • 1 1 ll,2[o t ](i 3 )(s -y e) 
to achieve a bound in terms of We now show that the schematics on line 

(|3.51fl are preserved after each round of Lie differentiating by inductively establish- 
ing the identities: 



(3.52a) (C x J) L = q-Y\ - J £ k ~ • (Xfc)i(«) , 



it=o r 



(3.52b) (£ X J) L = q-J2 ' (Xfc)l («) , 



(3.52c) (/£ J) A = %P ' (x*)a(«) 



fe=o r 



where the f2^ k and (Xfc)a depend on the specific combination and number of vector- 
fields, but have the same properties as the Sl a and x m nne <|3.51[) above. Also, 
the product notation is symbolic and is used to denote a finite sum of products of 
functions with these properties. Assuming that l|3.52|l is true for |/| = Z — 1 it can be 
shown that l|3.52|l holds for |/| = I through direct use of the identities H2.32|) - H2.33|) 
as well as the formula: 

(CxJ) a = X(J a )-J{[X,e a ]) . 
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We leave this as a straightforward, although rather tedious, exercise for the reader. 
Note that everything can be put in terms of the variables u, r, lo so the only thing 
to verify is the homogeneity claim with respect to the power of r. 

It remains to prove the second bound (|3.50() above. This is done in two separate 
steps. First of all, a direct computation involving the formula: 

{C x F) al3 = X(F a0 )-F([X,e a ],e )-F{e a ,[X,e ]) , 

as well as the bracket identities Ij2.31a|l and (|2.31c|) shows that at time t = one 
has the point-wise bound: 

£ {{E^n 2 + \H^ X F)\") < £ (l + r 2 f\(\Vl x E\ 2 + \Vl x Hf 

\I\<k \I\^k 
X6L 

Using the field equations l|1.10fl . this last line can be further reduced to the estimate: 



£ (|£(^F)| 2 + |H(^F)| 2 ) < (l+r- 2 ) m (|V^| 2 + |V^| 2 + |V^ 2 

\I\<k 

+ E (l + r 2 ) m+1 (|V^J| 2 + | 9 (F)MV^(X + )'| 2 ) 



|/|^fe |Z|<fc 



\I\Ck-l 



Multiplying through by (1 + r 2 ) s+7 and integrating this last line over R 3 , we have 
achieved the bound 13.5011 modulo the estimate: 



V f (l+r 2 y+~<+^\VlE\ 2 dx < V (/ (l+r 2 ) s +^+l J l|V^| 2 cte 

+ || (l + rr+^V^oCO) ||^ f ) . 

Using essentially the same steps which were employed to reduce estimate ()3.43f) 
above to (|3.45|l . this last line is a consequence of the following generalization of 
estimate (I10.2|) in the Appendix: 

(3.53) £ / r 2(, +7+m) | vv / ( l Jo(0) 



—!—) 1 2 dx 

4nr n 



< c, £ lk i+7+|7| viJo(o)|| 2 Ll 

lll<fe 



Estimate l|3.53|l can be reduced to estimate (|3.45f) through a process of induction. 
To see this, we assume that < |/| and integrate by parts a couple of times on the 
left hand side of (|3.53(l . This yields: 



(L.H.S.)£53 = E ({s+~(+\I\+l){s+l+\I\) [r^+i+W-V \vl(lj o (0) + -JL 
/ r3 r 2 <^D (viJo(0)).Vi(ij o (0) + ^-) dx 



1 M 2 
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Applying Cauchy- Schwartz and Holders inequality to the right hand side of the 
above expression in conjunction with the L 2 <^-> L 6 Sobolev embedding and setting: 

A(k) = Y, lk s+7+m vv^lj (o) + £;)lU 2 , 
B{k) = lk s+7+|/| v^Jo(o)|| L | , 

|7|<fc 

we have for 1 ^ k the estimate: 

(3.54) |(L.H.S.) B5^ | ~ A 2 (k) < A 2 (k-l) + B(k)- (A(k) + A(k - 1)) . 

Notice that with this notation estimate (|1U.2|I reads A(0) < B(0). This assumption, 
together with the inductive estimate l|3.54[l shows that A(k) < B(k) for all ^ k. 
This ends the proof of (|3.53|) and therefore the proof of l|3.48ll . □ 



4. Fixed Time and Space- time Energy Estimates for Complex Scalar 
Fields: Conformal Modifications and the Morawetz Theory 

In this section, we prove estimates of the type (|3.2(J|) for the complex scalar field 
(II. 2|) . This is necessarily more involved than in Section [3| because the tensor (|1.13f) 
is not trace- free which is a reflection of the fact that, unlike Maxwell's equations, 
the scalar wave equation is not conformally invariant. However, it is well known 
that there is an analog of the estimate (|3.4() for the wave equation which is also 
known as the Morawetz estimate, and was first used in |13j to prove local energy 
decay for solutions to (|1.13|1 in the case F = 0. We record this estimate here as: 

(4.1) / T\ \L<P\ 2 + T 2 \m 2 + (4 + T 2 _){\?4>\ 2 + |^| 2 ) dx 

J{t=t }nR 3 r 

< f (1 + r 2 )|V0| 2 dx . 

J{t=o}nK 3 

Notice that we have used the usual derivatives above because we are assuming 
that the connection is flat. The usual procedure for proving (|4.1() involves first 
coming up with a certain weighted L 2 type identity, again known as the Morawetz 
identity, and then performing several integration by parts in order to ultimately 
arrive at the estimate (14.11) . In the context we are working in, where we wish to 
prove energy estimates of the type (|3.2UI) involving fractional weights as well as 
characteristic and space-time energy estimates, such a procedure would become 
unduly tedious. Therefore, we provide here a new approach to Morawetz type esti- 
mates which provides l|4.1|) directly in divergence form. Having done this, it will be 
straight forward to modify our divergence identity to include various weights of the 
type H3.15f) - I|3.17[) as well as the fractional Morawetz field (|3.21f) . Furthermore, our 
procedure leads to a more direct geometric insight as to the nature of l|4.1|) which 
is obscured by the usual integration by parts proof. 

Our starting point is to conformally modify the Minkowski metric g a p in such 
a way that the scalar field equation l|1.2fl is preserved. As is well known, if one 
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performs the conformal change of metric: 

for some weight function Q on M, then any solution (b to the inhomogeneous 
equation: 

(4.2) U C (f) = G , 

will transform to ip = flip, where ip is a solution to the inhomogeneous conformal 
scalar field equation: 

(4.3) £&ib - -Rib = fl 3 G . 

6 

Here D c is the covariant wave equation on the space AfxC with connections (V, D) 
and R the corresponding scalar curvature. One readily calculates that: 

(4.4) R = 6ft 3 V tt V a (^) . 

In order that (|4.ri|) match up with 1)4. 2[l . we require to be such that the scalar 
curvature vanishes, R = 0. In light of the calculation 1)4.4(1 . this will be guaranteed 
if one has = where □ is the usual D'Lambertian on Minkowski space. There 

are two interesting and useful choices of £1 which give a singular solution to this 
problem and which are ultimately responsible for the estimate (I4.1|) . These come 
from the fundamental solutions to the Laplace and wave equation respectively: 

'il = r, "n = uu . 

We label the resulting conformal metrics by I g and 11 g respectively. Notice that 
both of these metrics are singular along the varieties r = and |t| = r respectively. 
However, these singularities will not effect what we do here because there will always 
be extra factors involving positive powers of £7 in all the identities we use which 
will cancel the singularities off. A striking property of the metrics I g and 11 g is 
that the vector-field A'o becomes Killing with respect to both of them (away from 
the singular set of course). This is seen simply from the identities: 



(uu) 2 J (uu)' 



Thus, for example one has that: 



(4.5) Cxjg = K (^)g+^C Ko g = 0, 
with a similar calculation showing that: 

(4.6) C Ko n g = . 

Note that the vector-field T is Killing with respect to ! g, but is only conformal 
Killing with respect to 11 g. The identities H4.5[l - i|4.6[) show that if one is to use 
the vector-field Kq to produce energy estimates, it is best done with respect to 
the metrics I 'g and 11 g instead of the usual Minkowski. Accordingly, we define the 
conformal energy-momentum tensors of the first and second kind associated to (14. 2|) 
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to simply be the usual energy-momentum tensors of the equation (|4.3|) with respect 
to the metrics I g and 11 g and the appropriately transformed solutions: 

(4-7) z Q afi [if>] = M(D a (rt)D^)-± I g a0 D^rt)D^j , 

(4.8) U Q a pW\ = ^(D a {uuct>)De{uu<t>)) - ± n g af3 D"<(uu<t>)D y (uu<f>) . 

Here we have used the notation D 1 = g ai D a for g — ! g, 11 g on lines 14.7( 1 and (14.81) 
respectively. As a direct consequence of the equation l|4.3|l and the divergence-fee 
property of the energy-momentum tensor for scalar fields one has the divergence 
laws: 



(4.9) V'Q^] = r 4 (^(G-^D (r^) + F^(^(r^ , 



(4.10) n V an Q a0 [<t>] = {uu) A (^(G- —Dp(uu<f>))+Fp 7 S(<P—D^(uu(f>)) 

\ MM UU 

where 7 V and n V are the Levi-Civita connections of I g and 11 g respectively. Since 
the vector-field Kq is Killing with respect to both of these metrics we may contract 
it with the tensors H4.7(l - I|4.8(l to obtain momentum densities which satisfy the 
divergence laws: 

(4.11) 



'V Q (tQ^MKof) = r 4 [M(G-±D Ko (rct>)) + (K Q fF ^(4,D^yj , 
(4.12) 



ff v" ("QaMiKof) = H 4 — DK^uu^ + iK^Fp^— W(uu<f>)) 

\ ' / V uu uu 



These last two identities can now be integrated over various space-time regions to 
obtain positive quantities for the scalar field 4> a t the cost of estimating the source 
terms on the right hand side of ((4. 11(1 - 1(4. 12(1 . To calculate these, notice that the 
volume forms of the metrics I g and 11 g are: 



1 „, 1 

Juu) 4 



(4.13) dVi = —dV M , dV n = j^cIVm 



while rT and uuT are the respective Lorentzian unit normals to the time slices 
t = const. Furthermore, notice that the vector-fields rL and uuL are the respective 
Lorentzian unit normal to the cones u = const. Therefore, applying the geometric 
Stokes theorem to with respect to these quantities, we arrive at the basic Morawetz 
type energy estimates for the complex scalar field (14.211 : 
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Proposition 4.1 (First and Second Morawetz Estimates for Complex Scalar Fields). 
Let O = r, uu. Then one has the following estimates for solutions to the inhomoge- 
neous equation (|4.2[1 : 

(4.14) / i^ 2 |^L(^)| 2 +« 2 |^L(^)| 2 + (^+« 2 )|^| 2 )^ 

\ (u 2 \^D L (Q0)\ 2 + U 2 |^| 2 ) dV c(u) 
(u)n{Q<t<t } 4 \ J 



<{t=t 
sup 



< / ° / \G ■ ^D Ko (Jty)| + |(A' ) /3 F^3( ( / ) iD7(^))| daaft 

JO JR 3 " " 



J{t=o}nR 3 z i2 



'{*=' 

To derive l|4.1|) from estimate (|4.14|l above is a simple matter. Assume now that 
both G = F = 0. Notice first that it suffices to prove jUJ with the r + , r_ weights 
replaced by u, u and (1 + r 2 ) on the right hand side replaced by r 2 . This follows at 
once from expanding the usual energy identity in a null frame, and using the fixed 
time Poincare inequality: 

(4.15) f |-| 2 dx < f \d r (f>\ 2 dx , 

J{t=const. }nR 3 r J {t=const.} 



to deal with the non-differentiated term. To deal with the homogeneous version of 
estimate (I4.1II . we simply use the identities: 

u 2 \-L(rcp)\ 2 = luL^+^l 2 , u 2 \-L(r4>)\ 2 = 1^0-— | 2 , 

r r r r 

u 2 \—L(uu(f))\ 2 = \uL(f> + 2(j)\ 2 , u 2 \—L(uu(f>)\ 2 = \uLcj) + 20| 2 . 
uu uu 

together with: 

(4.16) ^-y- = (uL^> + 2(f) - (uL<f> + — ) , 

(4.17) =^ = {uL<$> + 20) - (uL(j> - ^) , 
which collectively imply the two pointwise estimates: 

u^L^ + u^m 2 < 2(u 2 \h(r^\ 2 +u 2 \^U^)\ 2 )+^(u 2 +u 2 )\^\ 2 , 

(u 2 +u 2 )\^\ 2 < 2 £ « 2 |lL(^)| 2 + M 2 |iL(^)| 2 . 

Combining these last two, we have completed our demonstration of (|4.1|> . 

We now prove a fractionally weighted and space-time generalization of (|4.14() 
which will be our analog of estimate 13.20|l for the complex scalar field (|4.2() . This 
is: 



Proposition 4.2 (Generalized Morawetz Estimate for Complex Scalar Fields). 
Let 4> be a complex scalar field which satisfies the equation 1)4. 2|) . and let F be the 
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curvature of the corresponding connection. Let < e, 7 be given parameters, and let 
s be chosen so that ^ ^ ,s ^ 1 as well as e ^ s — |. Define the weights To, w 7 ,u> 7i£; 
and wty as on lines (|3.14[1 - l|3.17[) . Also, define the generalized Morawetz type energy 
content of <f> in the time slab {0 ^ t ^ to} l~l R 3 : 

(4.18) 

sup f (rl^-DLir^+r^lDL^+r^flp^ + ^A^w.dx , 

(rl^DUr^+r^lM' + rl^Aw, dV c(u) , 
(u)n{0sct<t } V r W J 

!^LM| 2 +r 1+2e (^•|D^| a +^W| a +7^'|^| 2 ))«/ 7i£ dxdt 

Then one has the following general weighted energy type inequality: 

(4.19) E^ s \0,t M < C?, 7 [|||J'|||i. [0tto](er B< e ^(0,to)M 

r| s r_ |G| 2 w 7e dxdt + / (l+r 2 ) s+7 |£>(<£)| 2 da; 

./{t=c 



to •/{ 

sup 

u 

rto 



where we have set ||| F |||ioo[o,t ](e) egwa/ to £/ie time-slab L°° type norm: 
(4.20) 

III^IIU~[o,t ](£) = sup (T+T ~ 2e |a| + t^tIt " £ |p| + r+r_|a| ) • u; 7 , c (u; 7 ) _1 . 
0^ t ^t v 7 

iG{i}xE 3 

In particular, if we further define the multiindexed L°° norm for electro-magnetic 
fields: 

(4.21) 

III XT' 1 1 12 I 2s+3| |2 1 2 2s+l| |2 , 2s+2 /~2 , 2\\ 

0^ t 

x<£{t}xR 3 

+ \q(F)\ 2 + \\Tl +1 Tl(w% e a\\l HLoo)[0ito] , 

where the quantity q(F) is defined as on line l|3.6[) . then under the more restrictive 
condition that 4e ^ s — ^ we /iawe £/ie bound: 

( 4 - 22 ) III IIU°=[0,to](e) < HI- F lllz°°[Q,to](s,7,e) • 

In particular, the estimate (I4.19|) remains valid with (|4.20|l replaced by (I4.21|) . 

Remark 4.3. Comparing estimate (|3.2(J|) with (|4.19|) we see that there is a close 
analogy, at least with respect to energy estimates, between the various null compo- 
nents of F and the gradient D<f>. Schematically, these are: 

(4.23) a ~ -D L {r<P) , a ~ D L <\> , 

r — 

(4.24) a ~ pj, p ~ ^ . 

r 
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As we shall see in the following two sections, this analogy persists in the discussion 
of L°° type estimates. We will make solid of this in the sequel, where we will use 
the analogy i|4.23[l - (|4.24[l to reduce the field equations l|1.17fl and its commutators 
to a equation for a single abstract vector quantity \& which has weighted energy 
and peeling properties equal to the null decompositions of F and D<f>. 



Remark 4.4. Note that the L°° norm l|4.21[l contains several terms not present in 
the original L°° norm (|4.2U|1 . Since the more involved norm will be used exclusively 
in the sequel, we have stated it here for future reference. Also, notice that on line 
(|4.22|l we have used the formulas Ij3.12|) and (|3.10(l for the subtracted charge p. 
Note that this bound (for the charge) becomes sharp in the far exterior r — t ~ r. 



Proof of estimate Ij4.19|l . The proof is nearly identical to that of (|3.2()(l in the pre- 
vious section. Our first step is to come up with a fractionally weighted momentum 
density associated with the first conformal metric I g. We define this to be: 

(4-25) = 'Q^MKlf-w^, 

where w lt<L is the weight function from line l|3.26fl . This satisfies the divergence law: 



(4.26) 'v a/ Pj>"' e >M = r 4 ^(G-i^(^)) + (l?ofF ;37 3(^)) 

Integrating both sides of the above expression with respect to the volume dVj — 
r~ A dVM over regions of the form: 

n(t Q) u ) = {o^t^t }n{u^u } , 

and noting that rT and rL are the respective Lorentzian unit normals to the time 
slices t = const, and cones u = const., we arrive at the following integral identity: 

(4.27) 

I {t=Q}n{-u ^r} J{t=t }n{t -u ^r} 

+ [ r~ 2 I P s L n ^[(j)\dVc ( u )+ I I r~ 4 (R.H.S.)633dEd*. 

We now calculate each term in (|4.26|1 individually. We will use the same ~ notation 
as in the proof of l|3.20|l . At time t — 0, we compute that: 

(4.28) r" 2 I P (s,7,£) M ~ (l + r 2 ) s+ ''|^(^)| 2 . 

r 

Likewise, at time t = to we have that: 

(4.29) r- 2I P^&] ~ ^|I CL(r0) |2 +T 2 S| I %M| 2 +T 2 S|w|2 y w7j 
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and on the cone uo — const, we have: 

(4.30) r- 2 'P^M ~ ( rl^D L (r^+r^\Pcb\ 



It remains to calculate the expression r~ 4 fR,.H.S.) (|4.26(l . We do this for each 
term separately. Since the first such term does not have a sign, we put absolute 
value signs around it and estimate: 



(4.31) \U(G--D w ( r( j>))\ < IGl-jr^l-^MI+r^l^l + r^l^l) , 

T q \ r r J 

(4.32) \{K s Q fF M ^DH)\ < rf +1 (\a\ ■ \*\ ■ \pcj>\ + \p\ ■ \*\ ■ \-D L (r<f>)\) 

+ T + T^(\a\-\ < t\-\M + \p\-\'t\-\D L 



To calculate the second term on the right hand side of expression (|4.26(l we need 
to compute the deformation tensor ^^Ttap- Using the identities (|2.3|) and (|3.22|) 
this is: 

^^(1) 5+ 1(^, 



2 f v — v 



Contracting this last line with |r 4I Q[4>] ■ w 7 , c yields: 

(4.33) Ir-^Q^M^n^-w^ ~ - (v + *)) • | ■ Wj . 

In particular, we see that this term only contributes a positive addition to the right 
hand side of l|4.26[) . To compute the last two terms of r~ 4 fR.H.S.') (|4.2fci[) . we simply 
use the calculations H3.28(l - I|3.29() and the expansions l|2.26[) to conclude that: 

(4.34) 

-\r- 2 'Q^MKof-Liw^) ~ C^r o 1+e ^|^(r0)| 2 +rf|^| 2 )-^ 
(4.35) 



7,e i 



Now, inserting the terms Q4.28|) - (|4.35ll in the identity (|4.27f) and removing the 
positive contribution (|4.33|) form the right hand side of the resulting estimate, we 
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arrive at our first preliminary estimate to (|4.19|l : 
(4.36) 

{l + r 2 Y+^\-D(r4>)\ 2 dx + f" [ \G\-(Tl a \-D L (r<P)\+Tl s \D L 4>\+Tl s \^\\ Wl dxdt 



{t=o}nR 3 



to 



JS. 3 
to 



r 2s+l 
' + 



|a| • |£| • \pd>\ + \p\. \Z\ .\-D L {rct>)\)- Wl dxdt 



J o ° J r + r 2 _ s (\a\ • |^| • \p4>\ + \p\ ■ |£| • \DjJ>\\ ■ w, dxdt 



la 



sup 

0<t<t o J{t 



-2s I 



\-D L ( r( f>)\ 2 +t 2 _ s \-D l H)\ 2 +4°^ )-w 7 dx 



sup 



( T l s \-D L (r4>)\ 2 + t 2 _ s \P4>\ 2 ) ■ w, dV c(u) 
u)n{o<t^t } \ r / 

?\~D L (r<j>)\ 2 +r 1+e ^|i/?i(r^)| 3 + 7i-|^| 2 ^ dzdi 



By repeating the above steps with the density l|4.25|l replaced by the following 
weighted version of the usual energy density for <f>: 



and using the bounds: 



m{G-D T <j))\ < (L.H.S)g3U) 



(4.37) |Fo 7 9?(0-Dt^)| < (L.H.S)63U 

we also have the estimate: 



sup 

o^t<t J{t}nR 3 



\DL(j)\ 2 ■ Wy dx 



(4.38) (R.H.S.)g3U) > C'l 

\D L cf>\ 2 ■ Wl dV c{u ) + [ f r 1+£ \D L $\ 2 • w' e dxdt 
,} Jo Jw 3 



sup 

« Jc(ti)n{o^t^t } 



By adding together estimates l|4.36|) — 1|4.38)) . we have almost achieved the state- 
ment of estimate <|4.19[1 . What is missing from the right hand side are terms involv- 
ing the weighted quantity u 2s 1 2 1 . As with estimate l|4.14[l , bounds on this quantity 
will come from integrating a divergence identity involving the tensor (|4.8() and then 
combining the resulting estimate with l|4.36|l above. To do this, we form the second 
weighted momentum density: 

(4-39) n Pt n ' e) i4>] = n Q a0 m(Kof-w s ,-y,e, 

where the weight function w s ,-y,e is given by: 

(4.40) w sn , e = (l + u) 2s - 2 .(( 2 - u ) 2 V(-«) + l-X + (-«)) 

+ (l + ur- 2 -^-((2-u) 2 ^ X +(-u) + l- X + (-u)) 
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We now follow through the same steps used to prove (|4.36(l above. Computing the 
divergence of Ij4.39|) we have that: 

(4-41) 

iiyanp( s , y ,e) [(j)] = ( ^ )4 fao . ±-D Ko (uu<j>)) + (K f F Pl ^(cj>D^)\ ■ w s ^, e 



Integrating this expression with respect to the volume dVn — (uu) 4 cIVm and using 
the Stokes theorem we arrive at the integral identity: 



(4.42) 



(uu)- 2D P^' e) [4>]dx = / (uu)- 2n P^ e \<t>]dx 

+ [ {uu)- 2 11 P s L r <>' dVc („ ) + / / (uu)- 4 : li.ll.S. i:fT~m: drill . 

•/c(« o )n{0<t^to} "~ J Jn(t ,u ) 

We now compute in turn each term in (|4.42(l above. When t = 0we have that: 
(4.43) (uu)- 2 a P^' e) [ct>] ~ r 2 ^\^D(r 2 ^)\ 2 . 

Likewise, at time t = to and along the cone u = const, we compute that: 
(4.44) 

(uu)- 2 vp^ty] ~ ( r 2s - 2 u 2 \±-D L (uu<P)\ 2 + Tl°- 2 u 2 \^D k (uu<t>)\ 2 + t 2 + °- 2 u 2 m\ 2 



uu uu 
(4.45) 

(uu)- 2 *P$>~i'%] ~ ^-V|^(u^)| 2 + 4 s -V|^| 2 )-«; 7 . 

We now turn our attention to the terms in the expression (uu) -4 (R.H.S.) H4.41|I . 
First, using the identities (|4. ltipi — (|4. 17|) with L, L replaced by Dl, Dl_ we conclude 
that: 

\uD L ^ + 24>\ ^ u\-D L (rct>)\ +u\*\ , 
r r 

\uD L( j) + 2(j)\ ^ u\D L 4>\ + 2u\-\ . 
— — r 

This then implies that one may estimate: 



3?(G • — D Ko (uu<f>)) -w s , 7 , £ < |G| • (rl s - 1 \uD L (f> + 2<P\ + t 2 + s - 2 u\uD l 4> + 2(t>\) ■ w 7 , 



uu 



(4.40) - C • ( T 2s \^D L (r<f>)\ + t 2 J\D l ^\ + r 2s |^| ) • ,r. 
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To bound the second term on the right hand side of (|4.41|) above, we simply use 
the estimate (|4.32H noting that one may trade tt* < r 2 ?: 

(4.47) 



{K fF^{(j>D^ 



w s 



<if +1 |a|.|?l-IM + IPl-l?l-l -Dl(t4>)\ • 



r + 7d*(|a|.|^|. \f><f>\ + \p\-\f\-\Dt<f>\)'W 1 



To compute the second two terms on the right hand side of (|4.41f) , first notice that 
a simple calculation similar to that used in lines l|3.28|l - H3.29|l above, which we will 
omit, shows that: 



-ii(S 7 , e ) > C-l T 2 + °- 2 T^w> lie , 



^L(w Jt€ ) > 



Notice that we only care that the second term above is non-negative as this will be 
its only use in our proof of (|4.19|l . In particular, the above two lines allow us to 
conclude that: 
(4.48) 



- ^2 "QlM( K o) P -£(&r,c) 
(4.49) 

_ ( UU J 



1 



^ C-^-'r^ u 2 \—D L (uu^\ 2 +u 2 \p<t>\ 
x uu — 



Combining the estimates (|4.43|) - (|4.49|) into the identity (|4.42|) and excluding all the 
terms on the right hand side of the resulting inequality except those which involve 
the Dl derivative and the characteristic term involving the Dl derivative, we arrive 
at our compliment to estimate (|4.36|l : 



(4.50) 



{t=0}nR 3 
to 



r 2( S+7 )| 1 D{r 2^\2 dx 



+ 



[° [ \G\.(rl s \-D L m+r 2 J\D^\ + rl^ 
Jo Jk3 V r w 



w 1 dxdt 



JWL 3 



2s+l 



\a\ ■ |^| • \p<f>\ + \ P \ ■ |^| • \-D L (r4>)\ ) • w 7 dxdt 



+ 









4 


n 


r + rl s ( 






h Js 


S3 ' 




r 



\M + \p\ 



\D L < 



w 1 dxdt 



sup 



sup 



f r^ s 2 u 2 | — Dl l (uu(I))\ 2 w 1 dx 
J{t}nK 3 U M 



2s-2 2 



« JC(u)n{0^t^t o } 
cto 



it | — D L {uu4>)\ w 1 dVc(u) 
uu 



2s— 2 l+2e 2 



1 



— D^(uu4>)\ w' e dxdt 



10 JR3 UU 

We now form estimate (|4.19(l by adding together estimates (|4.36|1 . (|4.38|1 . and H4.50J1 . 
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while also using some further identities to simplify the resulting right and left 
hand sides. First, notice that by again using the expressions (|4.16|) - ()4.17(1 we can 
estimate: 



t 2 _ s \D l 4>\ 2 



^ s |-| 2 < r| s -V|— D L {uu<i>)\ 2 + \D L ^ + rl s \-D L {r q 
r uu ~ ~ r ~ 



rl s \ l -D L (rcp)\ 2 + rl s \^ < r* s " V |^ L ( U ^)| 2 + \D L </>\ 2 + r 2s |^ L (r0)| 2 



Thus, using the bound: 



r ^ S |^|2 «- ---- Inl- 



and the energy notation (|4.18l) , and expanding out the t = energy expressions on 
the left hand side of the resulting estimate we conclude that: 





\G\-( 


JO Jl 


I 3 \ 



t ~ r 



■w 1 dxdt 



2s+l 



o Js. 3 

to 



\a\ • R • \p4\ + \p\ ■ | J| • \-D L {r<f>)\ ) ■ Wl dxdt 



I " I r + r 2 J (\a\ ■ \^\ ■ \p</>\ + \p\ ■ |^| • \D L <f>\) ■ Wl dxdt 
Jo Jr 3 V r r ~ J 



Setting: 



L 2 [0,t ](i 2 )( S ,7) 



(1 + r 2 )^ |^| 2 + |^| 2 ) dx 
{*=0}nR 3 V r 



f ° f t 2s t_|G| 2 w 7 , e dxdt 
Jo Jr 3 



and using the notation i|4.20[) . the condition that e ^ s — i, and again using 
the bound w 2 ^ w~ f ^T-iv' y e along with several different instances of the Cauchy- 
Schwartz inequality we have: 



|G||U 2[0 , to] (L 2 )( s , 7 )-(^) (S ' 7 ' £) M + lll^llk-[0,t ]( £ )-4 S ' 7 ' e) M 

(l + r 2 ) s +^ (j^| 2 + |^| 2 ) dx 



{t=o}nR 3 

Dividing this last expression through by (E^)( s >'>'> c )[0] and squaring, we easily 
achieve the bound: 



E^[<f>] ^ C 2 ^ 



L2[0,t ](L2)(s, 7 ) + lll F llli~[0,t o ](e) ■ E (sn ' e) [(f)} 



(l + r 2 ) s +^ ( \Dcf>\* + \^ ) dx 
{t=o}nR 3 V r 



The estimate Ij4.19|) now follows from this last line and the following gauge covariant 
Poincare type estimate which follows from the same reasoning used to produce ()6.7|) 



STABILITY FOR CSF ON MINKOWSKI SPACE 43 



which we shall proved in the sequel: 



f (l+r 2 ) s+ ^| 2 dx < f (l+r 2 ) s+ T |D0| 2 dx . 

J{t=o}nK 3 r J{t=o}nR 3 



□ 



5. L°° Estimates for Electro-Magnetic Fields 

This section is the sequel to Section above in that it contains the second main 
set of estimates for the field strength F a p. Here we will prove L°° type estimates 
for this quantity without assuming that it necessarily satisfies the field equations 
dl.l7fl . That is, in this section we will only assume that F a p is an arbitrary two- 
form and will prove L°° type estimates for it at the cost of bounds on the energies 
(|3.4ti|l . This is done through the well known procedure of applying global Sobolev 
inequalities in conjunction with geometric identity (|2.1b|) for the component decom- 
position of the Lie derivative of a two-form. This procedure was first explored in 
0, and our treatment adds little except a bound on the L 2 (L°°) term in the norm 
(I4.21|l . This turns out to be a direct consequence of the space-time L 2 estimate 
(I3.48[) and the weighted Sobolev estimates (|1U.14|1(|1(J.15|> and (|1U.19[) . However, 
we choose to prove the entire L°° estimate from scratch here primarily for the sake 
of completeness and because we use the fractional weights t+ s , r 2s . 

Before we proceed with the main estimate of this section, we introduce two 
notational devices which will allow us to reduce things to a more simple form 
where the estimates (|l(J.14[) - (|l(J.15fl and <|10.19[) can be used directly. The first of 
these is a simple algebraic tool which will make dealing with various commutators 
more straight forward. This is the introduction of the "radial boost" field: 

(5.1) flor — Lj l iloi — ——{uL — uL) ■ 

It has the advantage of commuting with the rotation and scaling fields {£lij,S}. 
Furthermore, it turns out that f2or preserves the null decomposition l|2.23[) even 
after one reexpands f^or and passes back to the usual boosts f2,o. This remarkable 
fact 1 as well as the rule that we only need one instance of the boosts in the exterior 
according to 14|) will help to streamline the steps we take in the sequel. 

The second tool we use here is a dyadic decomposition of the distance to the 
standard light-cone u = 0, as well as the distance along its translates u = const. 
We will use this tool to localize all of our L°° estimates to shells adapted to these 
distances. This allows us to add extra powers of t + ,t_ at will because they just 
appear as constants in our localized estimates. The language for this is as follows: 

We first isolate on each fixed time slice t = const, the extended exterior region 
t < 2r. We then chop this into dyadic pieces depending on the distance from the 



It is not a-priori clear that this should happen given the fact that the boosts Qio cause quite 
a bit of permutation in the null decomposition 12.231 . 
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cone u = 0. We decompose the extended exterior region into the set of double 
spherical shells: 

Ji = { x e {t} x M 3 ,t < 2r | 2 l < \t- \x\\ + 1 sC 2 l+1 } , ieN. 

For fixed < i, each of these sets can be further decomposed into exterior and 
inter portions Ji = J^~ U J i ~, where J^ 1 contains only the points where < =pu. 
Note that these sets are connected, and that J~ empty for i sufficiently large. In 
general, we will write J for one of these connected sets. We also form a partition 
of unity adapted to the decomposition Ji which has the homogeneity property: 

\drX3 t \ < tZHJ) , 

where: 

(5.2) T+{Ji) = min t + , T -{Ji) = mm T - ■ 

Furthermore, we note that it is clear that we can also choose the largest J~ so its 
cutoff is supported where r ^ -A. When the context is understood, we will simply 
substitute x f° r XJi an d T +, T - f° r t+ (<J) , T_ {J) . 

Next, we deal with the light-cones u — const. The notation we use here is 
similar to what was outlined in the preceding paragraph, and we will only consider 
the portion of C(u) in the region where t < 2r. For each such cone (omitting the 
dependence on u), we define the dyadic conical spherical shells: 

li = {x e C{u) , < t < 2r | T^t+\x\ + l^ 2 t+1 } , ieN . 

We use the notation r+(2"j),r_(2j) in analogy with line (|5.2|l above. We also intro- 
duce a smooth partition of unity \ii adapted to the Ti which satisfies the bounds: 

In particular, notice that we have not cut off strictly along the truncated cone 
C(u) (~1 {0 ^ t ^ ^o}i but we may assume that each of these x are supported where 
t < 4r. 

Using these notations we now glue together the estimates (|10.14J) - 11().15|I . (|10.19[1 . 
and Hl().2t)(l - 110.21|l of the Appendix into a form which will be used in the sequel: 



Lemma 5.1 (Weighted global Sobolev estimates). Let f be a test function on 
R x R 3 7 and let u> 7 be the weight function defined on line (|3.15l) . Then on each 
time slice t — const, with 1 t one has the following (non-uniform) estimates for 
arbitrary weights i± £ 1: 

(5.3) 

sup \T 1 + +5+ rl +S -f\ 2 w 1 < II r + + ri-w% X*Y J (f) \\h {t<4r) , 

t<9r 

|J|<1,|J|<2 
|Z|+|J|<2 

xe{d r ,s,n r},Yeo 
(5-4) sup |rJ +5+ /| 2 < E ll4 + ^(/)lli 2 (,<3 t) , 

2r<t \I\^2 

XG{s,n i} 
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Furthermore, defining the truncated cones: 

C{u) = C(u) n {1 s$ t ^ t } n {t < 2r} , 

C(u) = C{u) n {1 s$ t ^ t Q } n {t < 4r} , 
we have the characteristic L°° estimate: 

(5-5) sup |t| +5+ /| 2 < £ ^ +XlYJ ^KHc M )- 

(M)ecw |iKi,|J|<2 L A C ^>) 

|J| + |J|<2 

xe{n r-,s},YGO 

In the above estimates, O = denotes the Lie algebra of the rotation group. 

Proof of (|5.3|) - (|5.5[) . These follow almost directly from their local versions ^1(J. 14| — 
(110. 15[) . I|10.19|) . and l|10.20|l - H10.21|) . We'll just deal with (E3J) here as the other 
two estimates follow from similar reasoning and are left to the reader. For a dyadic 
shell J, by combining estimates (|10.14|) - (|10.15|l in tandem (for example with the 
L q set to q = 3) and using the estimate: 

\r-drxf\ < \xf\ + \xS(f)\ + \xn 0r (f)\ + \xd r f\ , 

where x is a cutoff on supp{x}, we arrive at: 

sup rlr^xf? < £ ||x^y J (/)||l, . 

X |J|<1,|J|<2 
|/| + |/|<2 

xe{d r ,s,n 0r },Yeo 

Since both sides of this last estimate are restricted to a doubling of the shell J ', we 
can multiply both sides by the constant ( 1 7)t 2(5 ~ +27 ( 1 7) if t < r, and by the 
constant (J r )r 2 _ (5_ [J) if r < t. Doing this and then summing the result over all 
shells J in the extended exterior t ^ 2r yields the desired result. □ 

We now state and prove the main result of this section: 

Proposition 5.2 (General L°° estimate for electro-magnetic fields). Let F a p be 
an arbitrary two- form on R x I 3 , and let w 1 be the weight function defined on line 
(|3.15|l . Then in terms of null decomposition (12.23(1 and the general energy norm 
(I3.4ti|l . one has the following L°° estimate: 

(5.6) sup (r 2 + s+3 \a\ 2 + T 2 + r 2s+1 \a\ 2 + r 2 + s+2 r^p 2 + a 2 )) ■ w, 

i£{t}xl 3 

+ \\r^rl(w')i e a\\l nLr)[0M < £ (0 , to)^ F] . 

I*l<2 

In particular, using the charge modified energy norm (|3.47|) . and the general weighted 
L°° norm (I4.21|) we have that: 

(5-7) |||^||U~[o,to](w) £ E^' € \o,t )[F] . 
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Proof of estimate Ij5.6|) . In what follows, we shall only consider the estimate (|5.6|l 
in the region 1 $C t. For ^ t < 1 the result follows from a weighted estimate 
on the electric-magnetic decomposition <|1.9ll in conjunction with a simple weighted 
version of the usual Sobolev embedding. We leave the details of this to the reader. 

The proof can now be broken down into three sections. All of the components 
in the first term on the left hand side of (|5.6(l except a can be treated with fixed 
time energy estimates. The pure L°° estimate for a is a consequence of character- 
istic energies, while the mixed norm estimate follows from the space-time energy 
estimate contained in 13.48f) . 



Bounds on {a,p,a}. We first prove this assertion for the a component in both 
the extended exterior region t < 2r and its compliment r < ^i. In fact, it is just as 
easy to demonstrate the full bound: 

(5.8) sup tIt 2 _ s+1 (\E\ 2 + \H\ 2 ) ■ w 7 < 

X 

|/| + |J|<2 

To prove l|5. 8J) . we first combine estimates (|5.3(l and H5.4|l with S + — 0,6" = s and 
5 + = s respectively This yields the following bound for scalar components: 

sup t 2 + t 2s+1 {\E\ 2 + \H\ 2 ) ■ w 7 < 

X 

\\t± V ,*(X i Y j (E),X i Y j (H))\\1 1 . 

|7|<1,|J|<2 
|/| + |JK2 

xe{d r ,s,a 0r },Yeo 

To finish things up, we need only show the following point-wise bounds: 

|7|<1,|J|^2 
|/| + |JK2 

xe{d r ,s,n 0r },Yeo 

| (^(44^)^(44^)) | 2 ■ 

|Z|^l,|J|<2 
|/| + |JK2 
X£{di,S,n i},Y£O 

This, in turn, follows immediately from expanding out the vector-fields d r and Qo r 
in terms of elements of L, and then using the Lie derivative formula l|2.1b|) in con- 
junction with the commutator identity (|2.31a|) . The important thing to note here 
is that in this latter formula one always ends up with constant linear combinations 
of elements of the translation invariant frame {d a }. This allows one to write the 
bracket terms on the right hand side of i|2.1b(l directly in terms of constant coeffi- 
cient linear combinations of (E,H). 

Our next step is to bound the terms p and a in the extended exterior t < 2r. This, 
combined with estimate (|5.8|) above in the deep interior r < ^t, will demonstrate 
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(|5.6(l for these components. We now need to be a bit more careful regarding how 
the operation of Lie differentiation permutes null components. To take this into 
account, we will again make use of (|2.1bfl as well as the expansion identities (|2.33|) 
for Lie derivatives, with some additional help from the frame identities on line 
(|2.21fl . First, notice that to expand vector-fields of the form X r — a/JQ, the 
formula (|2.1bjl implies: 

(5.9) (C Xr F aP ) = J{C Xi F) a(3 + e a {J)X7F^ + e p {u i )X7F cn . 

Therefore, after some simple computations using the first of the special identities 
(the second will be used in a moment): 

(5.10) e A {J)u>i = , e A {Lo l )ujf = -5 A B , 

T 

we arrive at the following equalities: 

(5.11a) d r {p) = u l p{C di F) , n tf (p) = p(Ai/) , 

(5.11b) S{p) = p(C s F)-2p, n 0r (p) = u l p{Cn m F). 

and: 

(5.12a) d r (a) = Ju{L dr F) , %(<r) = a{C^ 3 F) , 

(5.12b) S(a) = p(C s F) - 2a , Q 0r {a) = oj l a{Cn tH F) . 

Notice that the set (I5.11f) is easily proved using the identities (|2.1b(l . I|5.9|l . I|2.21|l . 
and (|2.33|) , while the second set l|5.12|) follows from l|5.11|) and the duality formulas 
(|2~5|) and $TZQ. 

Next, we see that by differentiating both sides of the expressions l|5.11jl - l|5.12(l 
with respect to the rotations Oy and making use of the homogeneity bound: 

(5.13) |tyj(wfc)| < 1 , 

we arrive at the point-wise bound: 

(5.14) _ 

Y, iX'Y^p.a)]" < Y pHC'xC^+a^C^F). 

|/|<1,|J|<2 |I|^1,|J|<2 
|J|+|J|<2 |/| + |J|^2 

xe{d r ,s,n O r},Ye0 X£{di,s,n 0i },Ye<a 
This last estimate, combined with the scalar L°° estimate: 

(5.15) sup rl s+2 T. Wl {p 2 W) < Y, \K'4x I Y J M\\l,( t **)> 

|/|<1,|J|<2 
|/| + |J|<2 
X£{d r ,S,n 0r },Y£O 

which follows from the bound l|5.3|) with S + — s and 6- = 0, implies the estimate 
(|5.6I) for the components p and a in the exterior region t < 2r. We have now 
completed the proof of 1)5. 6|) for the components a, p, and a. 
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The L 2 (L°°) estimate for a. We now prove the mixed norm estimate: 
(5-16) ||< +1 rZK)ia||| 2(iao) < £ E^(0, t )[C x F] . 

|J|<2 

Due to the space-time energy estimate for a contained in the right hand side of 
(|5.16fl . it suffices to prove the two fixed time estimates: 

(5.17) ^ 

sup t 2s+2 t_ H 2 «/ < Yl \\<(™')IM£x£yF) \\h [t< ^ , 

t^Jr 

|/|<1,|J|<2 



|/| + |JK2 

X£{di,s,n i},Yei 



(5.18) sup r 2s+2 - 2 >| 2 < £ || (e{C x F), H (£ x F)) 



r<±t 



2 



\I\^2 
X£{S,Qoi} 



These in turn follow from the scalar estimates 11 (|5.3|) - (|5.4|) and the same procedure 
used above to prove the bounds l|5.8[) and l|5.15|l above. Specifically, notice that by 
the same types of calculations used to produce (|5.11() - (|5.11|l . in particular l|2.1bjl . 
(|2.33f) . I|5.9|l . and the identities (|5.1U|) . we have that: 

(5.19a) d r (a A ) = Ja A (£ di F) , 

(5.19b) fi«(eu) = a A (£n ij F) + [n ij ,e A ] B a B , 

(5.19c) S(a A ) = a A (£ s F) - 2a A , 

(5.19d) n Qr (a A ) = Lo l a A {£n (H F) + a A ■ 

These allow us to pass from a scalar estimate of the form (15.311 (with S + = s 
and 6~ = 0) in (|5.17|1 to the right hand side of this estimate which contains Lie 
derivatives. To do this we simply differentiate both sides of the expressions (|5.19(1 
with respect to rotations, and then express the resulting right hand side in terms 
of the collection {a, a{CxF),a{Cx£vF)} for X 6 L and Y € O. In doing this, 
if the extra angular derivative lands on a a component on the right hand side of 
(|5.19|l , we can use these formulas to re-expand the result using the above remarks to 
bound the resulting coefficients. If, on the other hand, the extra angular derivative 
lands on a coefficient on the right hand side of (|5.19|) we can bound things by using 
(|5.13|l and its companion estimate: 

(5.20) \n jk ([Q lm ,e A } B )\ < 1 , 

which easily follows from the homogeneity property l|2.19|l of the fields {e^}, and 
the fact that Vtfj = rVtfj where flfj is a function of the angular variable only. 



The L°° bound for a. We need only consider the case t < 2r, as the desired 
estimate in the compliment was proved on line l|5.8|l above. That is, we only need 
to prove this bound along the cones C(u) which were introduced just before estimate 



It is clear that one can replace in estimate 15.31 with 
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(|5.5(l above. Using that estimate with 5+ = s , we have the scalar bound: 

sup r^wM 2 £ E Wr+^X'Y^an 2 , v . 
(M)ecw |iKi,|J|<2 

XG{s,n 0r .},'Keo 

Notice that the u> 7 factor can be added in at will because it is a constant on the 
cones C(u). Using this last line, we see that is suffices to prove the point- wise 
bound: 

E i*^»i 2 < E Wx^Ft . 

|/|<1,|J|<2 |/|<1,|J|<2 
|Z|+|J|<2 |/| + |JK2 

xe{s,n 0r },Yeo xe{s,n oi },Yeo 

This in turn follows from the decomposition formulas (|5.19(l and the coefficient 
bounds (|f> . 1 rift and H5.2()(l . The details are left to the reader. This completes our 
proof of estimate (|5.6I) . □ 

We end this section with a discussion of the peeling properties of Lie derivatives 
of the charge portion F of the curvature F. This will be used in the following sec- 
tions wherever the quantities C X F ' a p occur because, while there are nice explicit 
formulas for these objects, they obey no useful space-time estimates. 

Proposition 5.3 (Peeling properties of the pure charge fields F a p). Let F a p be 
an electro-magnetic field, and let F a p be its pure charge component as defined on 
the line H3.10|) . Then for each multiindex I we have the following L°° estimates: 

(5.21) W£ X F)\ < Ci \q(F)\-T Q rfxt<r+i , 

(5.22) \a(C x F)\ , |p(/£F))| , \a(£ x F))\ , < Cj \q{F)\-r+ 2 ~ 
where all IeL 



Proof of the estimates (|5.21l) - (|5.22|) . The proof is a simple inductive procedure 
similar to what was done for J a in section [3] starting on line (|3.51() . We begin 
by making the inductive hypothesis that: 

(5.23) a(C x F) = g . £ . ( r _ t _ 2 ) , 

fe=i 

(5.24) a(£ x F) = q • £ ^ • ( x +£(r - t - 2) , 

fc=0 

(5.25) p(C x F) = q • £ ^ • ( x +)£(r - t - 2) , 

fc=0 

(5.26) a{C x F) = q ■ £ ^ • ( x +)£(r - * - 2) , 

fc=0 

where the are smooth functions of the angular variable whose C°° bounds de- 
pend only on I and the specific component being considered, and the (xt)l( s ) are 
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smooth functions of the single variable s which vanish for s ^ — 1 and satisfy the 
homogeneity bound: 

\di(xt)l\ < Cj s k -> , 

for 1 ^ s. Here, as in previous discussions, the product notation Q ■ x + is symbolic 
for a sum of products of functions with these properties. 

Now, from the formulas ((3.1U|) . it is clear that the assumption ((5.23(1 - 1(5.26(1 holds 
for \I\ = 0. Furthermore, simple explicit calculations using the formulas (|2.32|l . and 
which we leave to the reader, show that differentiating as scalars the quantities on 
the right hand side of the formulas 1)5.23(1 - 115.26(1 yields quantities with these same 
properties. This follows simply from expressing everything in the derivative in 
terms of the variables u,r,u. Therefore, it remains to show that one can induc- 
tively reproduce identities of this form after the bracket terms on the right hand 
side of formula 1(2. lbl) are taken into account. 

First notice that if the Lie derivative is with respect to the fields {ftij,S} then 
the null decomposition 1(2.23(1 is preserved, so the claim follows from the fact that 
the coefficients in the expansions ((2.33(1 can all be put in the form ■ tp(f) for 
some affine function (p. 

Next, notice that while the translation invariant fields {d a } do not preserve 
the decomposition 1(2.23(1 . their coefficients in the formulas ((2.33(1 introduce a factor 
which is bounded by =. Therefore, it remains to deal with the bracket portion when 
dealing boosts {iloi}- Notice that it suffices to limit discussion to the component 
a because the remaining components have the same general form in 1(5. 2411 - 1(5. 26(1 
above and the homogeneity of the coefficients in the formulas 1(2.33(1 guarantee that 
this general form is preserved after multiplication by them. The claim now follows 
from the explicit formulas (|2 . 33±|) and l(2.33h(l , because the only thing that can dis- 
turb the structure of 1(5.23(1 is when one of the other components ^(C^F), p(CxF) 
or ^(CxF) gets moved to the a position. Notice that this in turn can only happen 
when an extra factor of - is introduced, having the effect of causing the sum in 
formulas 1(5.24(1 - 1(5.26(1 to start with k = 1. This completes the demonstration of 



6. L°° ESTIMATE FOR COMPLEX SCALAR FIELDS 

This section is the companion to Section 0] It contains L°° type estimates com- 
parable to 1(5.6(1 for complex scalar fields in terms of bounds on energy norms of the 
kind ((4.18(1 . These will be proved as in the last section through the application of 
global Sobolev estimates and several geometric identities which reduce bounds on 
scalar derivatives to bounds on covariant derivatives. The device which ultimately 
enables us to do this the following simple, although far reaching, estimate of Kato: 

Lemma 6.1 (Kato's "diamagnctic" inequality). Let cf> be a section of V over Ai, 
and denote by \4>\ 2 — (<f>,<j))v- Then one has the point-wise inequality: 

(6.1) \X(\cf>\)\ < \D x <f>\ , 

for any vector-field X on M.. 
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Proof of ()6.1(l . This is a straight forward consequence of the compatibility condi- 
tion p.ljl . We compute: 



X 



= (Dx<t>,<t>)v + (0, 
2|0| 

Applying absolute values to both sides of the above equation and using the Cauchy- 
Schwartz inequality, the desired result follows. □ 

As a first application of estimate (|6.1[l . we prove the following generalization of 
estimates <|5.3[) - (|5.5[) to complex scalar fields: 



Lemma 6.2 (Weighted global Sobolev estimates for complex scalar fields). Let 

and let <f> be a section of V over M, and let w-y be the weight function defined on 
line H3.15f) . Then on each time slice t = const., with 1 ^ t, one has the estimates: 

(6.2) 

sup \rl +&+ Tl +S -^ Wl < Y Wri+ri-wlD^DUWh^r), 

t<2r 

|/|<1,|J|<2 
|/| + |JK2 

xe{d r ,s,n 0r },Yeo 

1+ S + A\2 ^ II J+ /V ± 112 



(6.3) sup |r| +A+ ^| 2 < W T + +D 

xe{s,n (H } 

Furthermore, defining the truncated cones: 

C{u) = C{u) n {1 s$ t ^ t Q } n {t < 2r} , 
C{u) = C{u) n {1 s$ t ^ t Q } n {t < 4r} , 

we have the characteristic L°° estimate: 

(6-4) sup |t| +5+ 0| 2 < \\4 +D xDh 

(t.x)eC(u) |i|^l,|J|<2 

|/| + |JK2 
XG{S-Q 0r },Y£O 



2 



In the above estimates, O = denotes the Lie algebra of the rotation group. 



Remark 6.3. Notice that we have used the vector- field S — fior instead of the 
individuals {S, fior} in the statement of estimate (16.4(1 above. This additional 
structure will be important in the sequel. 
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Proof of estimates (|6.2|> - (|6.4|) . The proof of these are virtually identical to the 
proof of (|5.3(I - H5.5|I . with the added twist that one uses (|6.1|) after each application 
of the scalar estimates l|10.14|l - (|10.15|) . (|1U.19I) . and (|ll).2U[) - (|ll).21[l to expand the 
scalar derivatives to covariant derivatives. The details are left to the interested 
reader. □ 

It will also be useful for us to have a version of the estimate 1|6.2[1 which is tai- 
lored to deal with the initial data (|1.18b|) in the context of the norms which appear 
on the left hand side of (|1.25() . 



Lemma 6.4 (Weighted Sobolev estimates for the initial data). Let and let <fi be a 

section of the hyperplane bundle {0} xl 3 xC and let ^ < Sq be a given parameter. 
Then the following weighted L°° estimates hold assuming that their right hand side 
is finite: 

(6.5) sup (l + r) s °+^| < ]T ||(l + rr- 1+ l J lD^|| i2(K3) , 

x 1<|J|<2 

(6.6) su P (l + rp + tH < ll(l + r-) S0+|J| -D^IU=(R3) ■ 

ve{9,} 



Proof of the estimates (|6.5|I - (|6.6J) . The proof of these is almost identical to the 
proof of estimate (|6.4|l above. We simply apply the estimate (|1U.14|I - H1U.15|I of the 
appendix in order to the functions x4>i X now being the cutoff on a spherical shell 
of dyadic distance from the origin, using the Kato inequality l|6.1|l after each step. 
Then, using the differential bounds: 

\dh\ < (l + r)-' J 'x, 
where x is some cutoff on the support of %, and adding together over all the (finitely 
overlapping) x w e arrive at the set of estimates: 

sup (i + r y° +1 * \d>\ < Yl II 0- + ry-^Wxt || i2(R3) , 

X l J l<2 

sup (i+r) s °+i|0i < J2 wa+rr+^D^u^ . 

X l-f|<2 

ve{3,} 

The second of the above estimates is already of the form (|6.6(l . To conclude Ij6.5|l 
from the first estimate above we simply need to eliminate the zero order derivative 
term. This can be done through the use of the following covariant Poincare type 
estimate which holds for ^ < sq assuming that <f> is both smooth and compactly 
supported (which we may do via a standard density argument): 

(6.7) f {l + r) 2s «- 2 dx < f {l + r) 2sa \D r( j)\ 2 dx . 

JR3 JR 3 

By an application of the Kato estimate l|6.1|l . this last estimate follows from the 
corresponding statement with cj) replaced by a real valued (smooth and compactly 
supported) test function ip, and the covariant derivative D r replaced by d r . This 
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in turn is achieved by integrating both sides of the following identity with respect 
to the measure dr du> and then applying a Cauchy-Schwartz: 

d r ((l + r) 2so -VV) = (2s -l)(l + r) 2s °- 2 rV + 

2(1 + r) 2s "- 1 r(p 2 + 2(1 + r) 2s "- 1 T 2 tp d r (p . 
This completes the proof of l|6.5|l - (|6.6|) . □ 

Before proceeding to the main estimate of this section, we first prove another 
set of preliminary estimates which will be extremely useful in the sequel. These 
are conjugated versions of the estimate . T|> above which also involve the optical 
weight t_ , and which we also call Poincare estimates: 

Lemma 6.5 (Covariant Poincare estimates for time slices). Let <f> be a section to 
V over Ai which has finite H 1 ' E ? R norm (see line <|1.21fl )■ and let D r denote the 
corresponding radial covariant derivative. Then for constants p, q such that both 
\q\ < p + 1 and —l<p one has the following estimate: 

(6.8) / t^tI^ 2 dx < [ T p _ +2 T q + \-D r {rcp)\ 2 dx , 

where the implicit constant depends on both p and q. More specifically, in the 
exterior region one has the estimate: 

(6.9) f r p _T q + \0\ 2 dx < [ T p+2 Tl\-D r (r0)\ 2 dx , 

Jt<r Jt<r r 

whenever < p + 1 + q and — 1 < p, and in the interior one has: 

(6.10) f T p _r q + \(j}\ 2 dx < f T p+2 T q + \-D r (r0)\ 2 dx , 

Jr<t Jr<t r 

for values q < p + 1 and —1 < p. 

Proof of estimate (|6.8|) . Keeping in mind the bound (|6.1|1 . it is clear that we only 
need to consider the incarnation of l|6.8|) with (f> replaced by the real valued test 
function ip = \(f>\ and D r replaced by the usual radial derivative d r . Notice that <p 
is both bounded and has finite H 1 '^ (non-covariant) norm. We now consider the 
regions t < r and r < t separately. 

In the first case, it suffices to prove (|6.8H with t_ replaced by 1 — u and r + 
replaced by 1 + u. Furthermore, since we are assuming that the left hand side of 
(|6.9() is finite we will proceed via a density argument that assumes tp is both smooth 
and compactly supported. This is because in order get things to work below we 
need to use the truncated mollification: 

(6-11) <fie = Xe ■ ft * 

where % e = x(e _1 •) for some 0(1) supported smooth unit function %, and f e — 
e _3 /(e _1 •) for some positive smooth function / with unit mass. Now, a standard 
argument shows this mollification is such that ip e — > ip in the space H 1 ^. There- 
fore, we are reduced to showing that for each fixed value of e we have the bound 
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(16.9(1 for ip e . In particular, we may assume from the start that ip is smooth and 
compact. With this in mind, we integrate the identity: 

d r ((l-u) p+1 (l + u) q (r<p) 2 ) = {p+l){l-u) p {l + u) q (r^) 2 

+ qil-uy^il + uy- 1 ^) 2 + 2(1 - u) p+1 (l + u) q d r (rcj))(r(j)) . 
with respect to the measure drduj over the region t < r to achieve: 

(6.12) / h p ' q (t,r)\ip\ 2 dx + / (l + 2r) q \<p\ 2 dV S 2_ 

Jt<r JS, 2 t=r '~ r 

(1 - u)P +1 (l + uyd r (r(j))(r(j)) drdoj , 

where h p,q {t, r) is the weight function given by: 

h p > q {t,r) = (p + l)(l-u) p (l + u) q + q(l -uf+^l + u)"- 1 , 

= (l-ti^l + tO" ((p + l) + g— - 
V 1 + u 

The conditions —\<p and — g < p + 1 now easily guarantee the existence of a 
constant C p . q such that we have the point-wise bounds: 

Cp,l(l-u) p (l + u) q ^ h p ' q (t,r) C p , q (l~u) p {l + u) q . 

Applying this last bound and a Cauchy-Schwartz to the square of l|6.12|l with the 
angular integral on the left hand side discarded, we arrive at the estimate: 

2 

(l-u) p (l + u) q \cp\ 2 dx' 

t<r 



< 



(l-u) p (l + u) q \(p\ 2 dx ■ I (1 -?i) p+2 (l + u) q \^d r {rtp)\ 2 dx 



>t<r Jt<r 

This proves the assertion (|6.9|l by moving back to the weights t + , t_ . 

We now prove the exterior estimate (|6.10|l . We will again proceed via a density 
argument in which we assume ip is smooth 12 and hence bounded at the origin 
r = 0. The estimate follows essentially the same steps as above. First, we replace 
the weights t + ,t_ by (1 + u) and (1 + u) respectively. The analog of (|6.12|l now 
reads: 

(6.13) / h p > q (t,r)\p\ 2 dx - f (l + 2r) q \p\ 2 dV § 2 

Jr<t JSj =r 

(l + u) p+1 (l + u) q d r (r(f>)(r(f>) drduo , 
where this time h p,q (t,r) is given by the expression: 

h p > q (t,r) = (l + u) p (l + u) q f-( P +l) + q^- 



12 The only thing one needs to check here is that || X<i ' r 1 ( L Pc — <fi) — > as e — > 0, but 
this follows from the vanilla Poincare estimate 14.151 which holds in the spaces H 1 ' . 
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The bounds — 1 < p and q < p + 1 now imply the existence of a constant C Ptq such 
that: 

C-l(l + u) p (l + u) q < -h™(t,r) < C P}q (l + u) p (l + u) q ■ 

The remainder of the proof now proceeds by squaring and using Cauchy-Schwartz 
as above. □ 



We are now ready to state and prove the main estimate of this section: 

Proposition 6.6 (L°° estimates for complex scalar fields). Let <fi be a section to 
V over M. with norm | • | 2 and compatible connection D. Let F a p be the curvature 
tensor of D. Define the k th weighted generalized energy content of <fi in the time 
slab ^ t $C to to be: 

(6.14) 4 W) (O,£ O )[0] = £ E^^t^D 1 ^} , 



\I\^k 

xeL 



where Eq ' 7,e (0, to)[<£] is defined on line (|4.18(l . Now define the time-slab L°° type 
norm: 

(6.15) 

lll0llli°°[O,t o ](s,7,e) = ll T + +lr -K)?,e {-D L {r(j>)Xt<2r+D L (i>Xr<\t) \\h(L<~){0,t o ] 

+ T + +3 \-^ D L(r<P)xt<2r + D L (j)x r< i t \ 2 ■ w l 



xeL 



w 



7 



Then recalling the definition of the energy norm (|3.47|) for the curvature F a p, we 
have the following nonlinear estimate for the field <f> under the additional assumption 
that e < s — i : 

(6-16) ||k|||!oo [0)to](s , 7ie) < E^^(0,t o M-(l + E^- e \0,t o )[F}) 

+ ^ w) (0,to)[F].|||^|||| 00[0 ^ ](s!7i£) . 



Proof of estimate l|6.16f) . To begin with we may assume that 1 ^5 t, as the compli- 
mentary case can easily be dealt with by a straightforward application of radially 
weighted Sobolev estimates and energy bounds similar to what we will use in the 
deep interior region r < \t. We now split cases according to whether t < 2r or 
r<\t. 



Estimate <|6.15[1 for Djjj) in the region t < 2r, and for all components in 
the region r < ^t. As in the proof of Proposition l5. 21 above, we will deal with the 
worst component estimate together with the deep interior estimate. These are both 



5 (HANS LINDBLAD AND JACOB STERBENZ UNIVERSITY OF CALIFORNIA AT SAN DIEGO 



a consequence of the single fixed time bound (where we leave the time dependence 
on the left hand side implicit): 

(6.17) sup rlr^\\D^+ £ |^| 2 ) ■ w, < 

^ s ^(0,* D )M-(l + lll^nii^ro^c,^) +^ S '^ e) (O,*o)[ J P]-||101ilio KS[o , toKs ,^ e) - 

To prove this estimate, we first make a preliminary reduction. Notice that by the 
homogeneity property of any IgL we easily have the bound: 

(6.18) y i^i 2 < i^i 2 + |^| 2 . 

We will now prove the bound 1)6. 17[) for the right hand side of (|6.18l) above. We 
begin by using estimates (|6.2[) and l|6.3[l together with weights S + = 0, S~ = s and 
S + = s (respectively) for the first term on the right hand side of H6.18[l , and weights 
5 + = — 1,5~ = s and S + — s — 1 (respectively) for the second term, which allows 
us to give the full bound: 

(6.19) sup tI^J+Wd^ 2 + \^\ 2 )-w- f < 

x 7 + 

y \\rtw^ D^DHli + E Wt+'tLwI DxtWl* . 

|/|<2 |I|<2 
VGL X6L 

The second term on the right hand side of i|6.19[) is immediately bounded by R.H.S. 
(|6.16fl . so we will concentrate on bounding the second. To do this we will need to 
deal with the commutator [D x , D] . By using the formula i|1.5fl and the basic Lie 
algebra formulas (|2.31(l . we arrive at the point- wise bound: 

(6.20) y \[D' x ,D]<j>\ < t+- y \£ X F\-\DU\- 

|/|<2 |/| + |JK1 

XGL X,Yeh 

Here we have used the notation l-C^-Fl = \E(C X F)\ + \H (C X F)\. Using the decom- 
position l|3.12f> and the L°° estimates l|5.21ll l|5.22|) , as well as the norm definitions 
(|4.21|l and (|6.15|) we can fold estimate l|6.20|l into l|6.19|l above and produce the 
bound: 

y Wr^D'xDHli < ^ W) M-(l + lll^lll^[o,* ]( s , 7 , e )) 

|/K2 

xeL 

+ £ -i S ' 7: ' S) [- F1 ]'lll ( / , ll| 2 ~[0,to](s,7,e) + III F llll~[0,t D ](s,7,6) ' E H T - * lll 2 (K 3 ) ' 

xeL 

The first two terms on the right hand side of this last expression are exactly what 
we see on the right hand side of (|6.17l) . Therefore, we only need to bound the 
last term. Unfortunately, this expression cannot be directly bounded by the energy 
(|6.14|l . This is where we need to make use of the fact that there is the additional 
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room < s — ^ 7 which allows us to replace the r_ 2 weight with the factor t s _ . 
We now prove the following estimate: 

(6.21) \\ T -~ lD xH\h m ;$ E (|I<mI^(^x0)II! 2(R 3) 

+ ||rlH|Z>i£»^||| 2(ffi3) + ||r|( U ;)|iz>^||| 2(K 3 ) ). 

This last line follows from Poincare estimate 13 (|6.8|) with p = 2s — 2 and g = 
along with the simple point-wise bound (note that the extra w 7 factor is just added 
in to match with R.H.S. I|6.17(l and is not needed for things to be sharp): 



(6.22) 




We have now completed our proof of estimate (|6.19|l above. 



Estimate (I6.15|l for 2 and p(j) in the region t < 2r. We begin with the undif- 
ferentiated estimate for ~. Applying estimate (|6.2II with (5+ = s — 1 and 8- = 0, 
and using the point-wise bound: 

(6.23) £ \D X 4>\ < T+ \±D L {rt)>)\+T4DLcl>\+T + (\M + \^\ \ , 

which follows easily from the expansions (|2.32|) . we arrive at the L°° estimate: 
(6.24) 

sup r^r^w, < W (^ s |-I?l(^^)| 2 + r! s |i^(^)| 2 



lsjt<2r r J t<i 



Notice that one can add another derivative to <fi on the right hand side of (|6.24() and 

n bou 
bound for „ . 



still remain bounded by the energy E^' 1 '^ [<fi]. Therefore we also have the desired 



To bound the term p<f> we use the following point-wise estimate which is valid 
in the region t ^ 2r: 



r+\p<t>\ < E i^v 



and then employ estimate l|6.24[) with - replaced by — — 



'We may assume that the H 1 ' 8 1 norm of <f> and Dx4> is finite or there is nothing to prove. 
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Improved L°° estimates for the term £ in the region t < 2r. This estimate 
follows from a simple combination of the global Sobolev and Poincare estimates as 
we now explain. First, we see that by applying l|6.2|) with S + = and J_ = s — 1, 
and expanding out the fields d r and flo r , we have that: 

SUp ^T 2 - 1 ^! 2 ^ < ]T \\Tt- l (w)faxDY<l>\\L»(.t<4r) ■ 

1S ** <2 ' |7|<1,|J|<2 

|/| + |J|<2 
XE{di,S,Q 0i },Y£O 

We are now finished by applying estimate (|6.8I) in the same way as used to prove 
1|6.21|) above, with the help of the point-wise bound l|6.22|l to convert the result into 
energy. 



The L 2 (L°°) estimate for -Di,(r<f)). We'll begin here by estimating things in the 
region 1 < t < 2r. Using estimate (|6.2|l with 6+ = s , 6- = and the weight w 7 
replaced by uu'^ e on each fixed time slice, we directly have the bound: 



(6.25) || T s + +x rl{w')i^D L {r4>) \\mL~)[i,t ]n{t<2r} 

1 



< 



E 



\I\ + \J\^2 

XG{d r ,s,n 0r },Yei 



Ti(vr)l e D 1 x Dt-D L (rq 



|L 2 (L 2 )[l,t ]n{t<4r} 



Keeping in mind now the space-time energy norm contained in E^' 1 '^ [</>], we are 
finished once we have taken into account the commutator [D^Dy, ^Dl{t-)\. Be- 
cause we are dealing with the component with the best decay properties, these 
commutators will need to be dealt with carefully. First, we use the formulas (|1.5|) 
and l)2.31b(l to compute the effect of the rotations (we also use the decomposition 



(6.26a) 



[D Qij ,-D L (r- 
1 



1 [n^aAiCa^F) + ^/fefty) a A 



(6.26b) Q k i([Dn v ,-D L {r-)]) = - 

r 

Similarly, using the formulas l|5.19|l and the fact that 9 r (f2^-) = r -1 ^, we compute 
the effect of the fields {d r , S, f^or} to be: 



(6.27a) 


[D r ,-D L {r-)] = 
r 


- V^ip- 


(6.27b) 


[D s ,-D L {r-)] = 
r 


V— 1m p — 


(6.27c) 


[Da 0r ,^D L (r-)] = 


\p-Yup + 


(6.27d) 


d r ([D nij ,-D L (r-)]) = 
r 




(6.27e) 


S([D Qii ,lD L (v)]) = 
r 


- V=ln3 


(6.27f) 


n 0r ([Dn tJ ^D L (r-)}) = 





u k a A {C dk F) 



(a A {C s F) 



OLA 



-a a 



u 
—( 

r 



STABILITY FOR CSF ON MINKOWSKI SPACE 



-,<) 



Combining the identities Ij6.26|l - (|6.27|) . and expanding out the fields d r and flo r 
using the identity L(u>i) = to move the contractions past the ^-D^(r-) derivative, 
we arrive at the following point-wise bounds for the commutator in the region 
1 t < 2r: 



(6.28) 



|J|<1,|J|<2 |J|<1,|J|<2 
|Z|+|J|<2 |/| + |J|<2 

X£{d r ..S,n r},Y£<0 X£{di,S,floi},YeO> 

+ J2 t+W£x*)\ ■ \ D U\ + E (r-\p\-\Dx<l>\ + ro\^ 
\i\+\j\^i ^ 

X,Y£L XeO 



We now insert each term on the right hand side of l|6.28|) into the right hand side 
of estimate (|6.25|) and bound the result using the norms (|3.47|) . Q6.14p . (|4.21|) . and 
(|6.15[) . Since there are some choices involved and the condition e ^ s — ^ makes its 
appearance here we do this explicitly. Note that the first term on the right hand 
side of (|6.28|) is already of the right form so we move on to the second. For this 
term, we split cases depending on wether the derivatives are falling on a or <j). In 
the first case, we estimate: 

E ll< +1 (^)i-|«(£^F)|-0||! 2(i2)[Oito] , 

X£L 

< E II T U w ')h <£>xF) Wh(Li) ■ II r+4> ||i~ (i oo )[0 , to] , 
|/|<i 

X£L 



In the second case, we use Holders inequality to conclude that: 



E \\r s + +1 {w')\, t .\a\-DU\\l^ mM > 
|/|<i 

< || 4 +1 rIK)t a lli»(L-)[o,to] ' E II T - iD x^\\ 2 L^(L^oM ■ 

X£L 



Notice that the first factor on this last line is contained in the norm (|4.21|) . To 
bound the second factor on the last line above, we again use the extra room ^ < s 
and apply the estimate l|tj.21[l . 

It remains to deal with the last term on the right had side of ({6.28(1 . To handle 
this we expand p — p + p. Using the identity 1)3.10(1 , we see that we may combine 
the term T-pD^cf) with the second term in this expression. Notice that this is 
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automatically included in the energy (|6.14(l . Thus, it suffices to use the bound: 

E ||^r_K)i^x0lli 2 ( i2)[ o, io ] < 

mo 

III F llll«'[0 1 to](s,7,e) ' E W T -( w ')l^ - D X<t>\\h(L2)[0,t ] ■ 

To control the second factor on the right hand side of this last expression, we simply 
use the condition e ^ s — h which guarantees that r_ ^ t^ s Tq +2< \ This concludes 
the t < 2r estimate for ±D L {r4>) in L 2 (L°°). 

To wrap things up for this subsection, we need to deal with the L 2 (L°°) estimate 
in the deep interior region r < \t. In this case, it is just as easy to deal with the 
entire gradient D<p as we did in the first subsection of this proof. Therefore, the 
bound we seek to prove is the following: 

(6.29) \\< +1 - £ DH\l HL ~ )llM{r<kt) < Ei s ^\0,t o M . 

Using the r < range restriction in conjunction with the point-wise bound: 

(6.30) C- 1 ]T \D X 4>\ < t-\D4>\ ^ C \ D ^\ . 

X£L X£L 

as well as the following interior Sobolev estimate on each fixed time slice (from 
ljfi~5J) with S+ = s-l- e): 

SUD J2 <'\ D ^\ < E \\T+*~'Dx<l>\\lZ(r<it) i 

r< \ t xeL 1<|/|<3 
we easily have fixed time bound: 

sup ^-'m < e \\T+*~ eDl x<i>\\ii(r<%t)- 

By integrating the square of this last estimate over the time interval (0,io)j and 
again applying the point-wise bound (|6.30|) . we arrive at the proof of (|6.29|l above. 
This completes our discussion of the L 2 (L°°) estimate in (|6.16ll . 



The pure L°° for iD^(r0) in the region t < 2r. We begin here by applying 
estimate (10.411 to the term ^D^rcf)) with 8+ = s. This yields the bound: 



(6.31) sup rfV^W)! 2 ^ £ 
(t,x)eC(u) r 

E \\rl(w)\ D I x Di l -D L {r^)\ 

|I|<1,|J|<2 
\I\ + \J\^2 

X£{s-(i r-},Yeo 
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To deal with the commutator, we apply the identities (|6.26(I - H6.27(I which provide 
the following refinement of (|6.28|) : 

(6.32) £ < E l-DLirD'xDU)] 

|7|<1,|J|<2 |7|<1,|J|<2 
|/| + |JK2 |/| + |J|<2 

xe{s-n 0r },Yeo xe{s,n i},Yeo 

+ E T + \a(£ x F)\.\DU\ + E (h-H-I^I + 1^ 
|j|+|j|<i |/|^i ^ 

Ye{s-n 0r }uo 

Notice that the weight r_ appearing in (|6.28|) is replaced by |u| in the last term 
on the right hand side of l|6.32[l above. This is because we do not have any term 
involving the radial derivative d r to contend with on the left hand side. It is nec- 
essary to have this if we are to use the characteristic energy bound contained in 
(I6.14|l . We now insert the terms on the right hand side of l|tj.32H into the right hand 
side of (|6.31[) . Notice that the first such expression is automatically bounded by 
the energy (|6.14() . Therefore, it suffices to treat the last two groups of terms which 
we do in reverse order. 

To bound the last group of terms on the right hand side of (|6.32[) , we expand 
p = p + ~p and group the term |u| • |p| • D^cf) with the term r~ 2 |u|_D;^-</> using 
the identity l|3.10|) . The resulting term is contained in l|6.14|l . so we are left with 
bounding: 



E 

AeL 



< \\rl(w)-p\\ 2 L2{d{u)) ■ E K^ S ^l| 2 L~[o, o] , 



A"GL 
L°°[0,to](s, 7 ,e) 



To wrap things up here, we need to bound the middle set of terms on the right 
hand side of (|(j.32|) substituted into the right hand side of <|6.31|) . Dealing first with 
the term where the derivatives fall on a, we estimate: 

E ii<+v)f>(^F)|.^ , 

< E KMfa(£^F)|| 2 L2f v \\r + ^\\l^ mto] , 

<^' 6) (0,io)[^-||^llli-[0 ) to](.,7«)- 



Therefore, it remains to deal with the case where the derivatives fall on <fi. Here 
we have to split cases depending on whether we have Dy4> with Y = S — fio, or 
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Y 6 O. In the first case, we write Y = uL and we expand: 

\uL4>\ T+ \-D L {rcl>)\ + T -±\4>\ . 

r r 

Using the fact that t < 2r this allows us to estimate: 

< || <h| «ll 2 L2(a(u)) ■ (lk^^Mlli» [ o > to] + ll^lli~[o,t ]) . 

<^^(0,to)[F]-||k|||ioc [0ito](5 , 7)e) . 



We are now left with estimating the expression YIypq T ^ +1 ( w )i ' M ' l-^V^I m the 
space L 2 (C(u)). This turns out to require an argument which is a bit more involved 
than what we have been doing so far. Notice that we cannot directly use either 
L°° or L? bounds on the Dq^ cj> term. Furthermore, a Poincare estimate similar to 
(|6.8I) for light-cones will not work in this situation because it introduces boundary 
terms which cannot be controlled. Therefore, being unable to deal with things on 
the light-cone C(uq), we extend the integral to the exterior truncated time slab: 



Tl(u , t ) = ([max{0, u }, t ] x R 3 ) n {u < u } ■ 



Doing this via Stokes theorem yields the estimate: 



(6.33) IK + V)?>|-^|| 2 L2f5 n < 



T 



2s+2 



L (|a| 2 • \D Y 4>\ 2 ) dxdt ■ iu 7 (u ) 



E (,, . 

+ llr^Hl-lal-^lll^^^ga)) . 

In the first term on the right hand side of the above expression, the weight w 1 {uq) 
denotes the constant value of the w 7 weight along the C(uo) cone. In the second 
term on the right hand side above we have switched to the variable iu 7 weight, 
which is possible due to the fact that this is a non-decreasing function for r — > oo. 
To bound this second term we estimate: 



E 

Y60 



rl +1 (w)l-\a\-D Y cj>\\l H{t=Uo}xI 



< 



E 

VGO 



F 



2 



[0,t ](a,7.e) ' 11 T + 



DyV 2 



L 2 ({t=uo}xWL 3 ) 



We are now in the familiar territory where estimate (|6.21|l can be applied to the 
second factor on the right hand side above. 
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Moving on to the first term on the right hand side of JEHU, we expand the L 
derivative to bound: 

(6.34) I / / r + +2 ^(M 2 ' \ d y0\ 2 ) dxdt -w 7 (u ) sC 



Yet 



Tl(u ,t ) 



YGO 

xe{9 r ,9 t ,s,n 0r .} 



n 


t\ s+2 tZ 1 I 


10 Ji 


!3 



Using the formulas i|5.19[) and the fact that d t (oLA) — aA{Cd t F), we expand the 



sum: 



(6.35) J2 \X(\a\ 2 -\D Y ^)\ < 



YGO 
X£{d r ,d t ,S,n 0r } 



£ M 2 -|^| 2 + Yl H-HC Y F)\-\D Y(I 



J|^2 |J|<1 
XeL AeO.YGL 

Substituting the first term on the right hand side of Ij6.35l) in the right hand integral 
in (|6.34|) . and using the mixed norm estimate for a we bound: 



J2 K+V^nl >l-^x<Hli 2([ cM ]x R 3) < 



m<2 



l^lllL»[0,t ]( S , 7 , e ) ■ Y \\ T - 2+e D 



X4> lll»o(l^ ! <|./„| 

|/|<2 



Using now the condition that e ^ s — i, we are in a position to again apply 
(I6.21|l noting that one may add an extra derivative to this estimate and still remain 

bounded in terms of E^'^'^iO, to)[(j>}- 

We have now reduced things to estimating the integral: 

I = Y I ° I t+^tZ 1 \a\ ■ \a{C x F)\ • \D Y <j)\ 2 ■ w 1 dxdt . 

AeL,YGO ^° R3 

Since this cannot be done directly, we use the identities: 

\Dy01 2 = 2Y($t{<i>, D Y <t>) - 25R(0, D 2 Y 4>) , 

2^,D Y cf>) = F(|0| 2 ) . 

to integrate by parts several times with respect to the -Dq ; . derivatives to achieve 
the bound: 

to 



1 ~ Yl 1° I T + +2t - X M • \ a (£xF)\ ■ 101 • \ d y<I>\ ■ w l dxdt 

|Z|,|J|<2 J ° 3 

AeL,YeO 

E r I < S+2 ^ H£xF)\ 2 ■ \<i>\ 2 • «> 7 dxdt , 



\I\^2 

AeL 
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To estimate the first integral on the right hand side above, we use a multiple Holder 
inequality and the simple estimate tZ 1 w 1 t^w' e to conclude: 

h < lk+ +1 ^(w')|,eQf|U=(i~ )[0 , to ] • ||r4.(w')|,ea(£^F) \\ L *(L')[o,t ] 

\I\.\J\^2 

veL,veo 

• II T + T -4>\\L~[<d,t Q ] ' II T - s+£ - D y0lU oo (£ 3 )[o,to] i 

< (lll^lll^ [ o,*o]( W ) + ^ w) ^) • (llkllli-p,*o](., 7 ,.) + ^' ,7,,) M) ■ 

Notice that the passage to the last line above is guaranteed by the condition e ^ s— 5 
and the estimate (|6.21() . To estimate the second integral -Z2 above, we proceed as 
follows: 

h < J2 IK(^)7,ca(^x^lli»(i»)[o ) to]-||T + rl0||l- [D)to ] , 

|/K2 

< £J^ s,7 ' e) [F] • HI ^|||ioo[ 0) t ]( a , 7 , e ) ■ 
This completes the proof of the pure L°° bound for ^D^rcj)) in the region t < 2r. 

Combining all of the estimates proved in the various subsections above, we have 
proved the estimate: 

\\\H\U [0 ,t o] (s,-y,e) < ^ 7 ' £) [O,to][0]-(l + lll^lll|o= [ o^( W ) + ^ 7 ' e) [°' i o]^]) 

+ (ill F llli<»[0,to](a,7,e) + E 2 ^ i° ^ ^[F]^ ■ \ \ \ (j> 1 1 li«.[o,t ](«, 7 , £ ) • 

Using now the L°° estimate (|5.6|l for the curvature, we arrive at estimate (|6.1(i|l . □ 



7. Abstract Weight Notation and the Bilinear Space-Time Estimates 

in General Form 

In this section we will organize and consolidate the L 2 and L°° type estimates 
which have been proved in the preceding sections. This will help to streamline 
notation for the remainder of the paper, and will ultimately reduce a lot of overlap- 
ping which occurs in the commutator estimates needed to deal with the nonlinear 
problem (|1.17|) . This will also give a chance for the reader to review the various 
notations, parameters, weights, etc. which have been introduces thus far. Our first 
order of business here is to set up some generic markers for the quantities which 
arise in the estimates we consider. These are the null components of the curvature 
F given by (|2.23|) . the corresponding null decomposition of the gradient Dcj), and 
the weighted scalar field As we have already mentioned in Remark |4 . 31 above . 
there is a straightforward analogy between the different members of these sets. To 
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make this analogy precise we form the following sets of field quantities: 

(7.1a) = {a, -D L (rcj>) Xt<2r + £>/>X r< | t } > *(-l) = {« > D L<P} , 

(7.1b) * {0) = {p, a, P<j>, Y, T^>' Q(0) = W- 

We will also use the following designation for sets containing the component de- 
compositions of Lie derivatives of the various field quantities listed above: 

(7.2a) /£* (1) = [J {a(4F), V(rD^) X t<2r + ^(0^)XKit}- 



DlD 1 - 



(7.2b) 4* (0) = (J , a(/^F) , p(D* x <f>) , £ _^J^ } 

|/|<fe |,/|^i,yeL + 



(7.2c) = (J , %(^0)} 

|Z|<fe 

(7.2d) 4g (1) - (J {a(£^)} , 



XGL 



(7.2e) 4g (0) = (J , ^(^F)} 



(7.2f) 4Q ( _ 1} - U {a(£^F)} . 



\I\^k 

While it is true that the above objects are properly sets, they can be manipulated 
as if they were functions by making each designator stand for the absolute sum of 
each element in the set, and letting products of the various sets denote the absolute 
sum of products of each element in the corresponding sets. As a first application of 
this principle we use the above notation to recast the norms 1)3. 47[) . 1)6. 14[) . Q4.2ip . 
and iKTTSjl : 

(7.3) 

4*' 7 ' e) (0,to)[*] = 4 s ' 7 ' e) (0,i )[F] + 4*"' e) (0,io)M , 
(7.4) 

lll*llli-IP > to](.,7,e) = \l( F )\ 2 + IK + VV)t,e*(l) lll^-MOA,] 

sup (^Tf-^Ha+^l^j^+^^l^jp+^'+V-l^o)! 2 ).^ . 

Notice that we have included the extra L°° norm for <f> on the right hand side of 
()7.4|) . This extra estimate turns out to be very important for us, but will be used 
in a somewhat auxiliary manner. 

We can now write out the content of line l|7. 3JI with the symbolic estimate (ne- 
glecting the characteristic energy terms which are only used to prove L°° estimates 
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and which don't quite fit into our scheme): 




We can also write the L°° estimates H5.210 - H5.22j l. H5.6j) . and H6.16j) in the following 
consolidated symbolic form: 



Proposition 7.1 (L°° estimates for field quantities in abstract from). Let ^ to 

be a given fixed time, and let < 7, e, s be parameters chosen so that e ^ s — \, 
then one has the following abstract non-linear L°° estimate for the field quantities 
and ClQ: 

(7.6) |||r£*|||i- [0 ,, oK ., 7 , e) < 4+2 e) (o,to)[*]-(i + 4' ,7,e) (o,to)[*]) 

+ ^ s,T ' e) (o,to)[*]-iii^iiiio, [0 , tD](s , 7>£ ) , 

and: 

(7.7a) ^lQ(i) < c k \q\ ■ t t+ 2 ■ xt<r+i , 

(7.7b) 4Q ( o) , < C k \q\ ■ t+ 2 ■ X t<r+i ■ 

We now conclude this section by proving all of the estimates we will encounter in 
the rest of the paper in a simple abstract form. It turns out that everything we 
need can be cast in the language of L 2 bilinear space-time estimates. These in turn 
follow directly from Holders inequality and the L 2 and L°° type estimates contained 
in the right hand sides of (|7~5|l and $TM - \7.7ty . 

Proposition 7.2 (Weighted L 2 bilinear estimates for field quantities in abstract 
form). Let ^ to be a given fixed time, and let < 7, e, s be parameters chosen so 
that e ^ s — o ; and define the auxiliary weight: 



(7.8) w 7 , e = rJxt<r + T_ t \r<t , 
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then one has the following abstract estimates for field quantities ^> and <I>: 

(7.9) \\r 2 + s+1 w^ * W ■ 3(0) Hi a[ o,t ] £ 

4 s ^' e) (0,to)[*(i)]-|||<f(o)lll!o. [0 ,t ](^, e ) , 

(7.10) \\rl s+1 w^ • 3 (0) ||| 2[0>to] < 

i?^ 7 ' e) (0,to)[3 ( o)]-|||3(i)|||io 0[0!to](s , 7!e) , 



(7.11) 



V-"^7,e 3(0) -3(0) ||i 2[0)to] • 



+ 

,T '(°)J ' HI IHi~[0,t ](s,7,e) ' 



4^(0,to)[3 (0 )]-|||3 ( ^ 



(7.12) || <+V^ e *(!)•$(_!) \\h [0M] < 

4'"' r,6) (0.*0)[*( l }]-|||*(-l)llli»[0, t o](.,7,e) 

(7.13) || T S + +1 T S _W 7 ,e ■*(-!) \\h [QM] < 

4'"' r,6) (0.*o)[*(-i)]-lll*(i)llli-[o, I 

(7.14) 



i,t ](s,7,e) ' 



T + T - 3(0) ' 3 ( _i) || L 2 [CMo] < 



4 S ' 7,£) (0,io)[3 (0) ]-|||3 M )|||ioo [0ito](s , 7ie) , 



(7.15) 



lk+ T - W 7,e 3(0) •$(_!) ||| 2[0)to] < 

4* ,7,e) (0,*0)[^(-l)]-|||*(O)llli»[0, t o](.,7,e) ' 

Recalling the definition of the weight function w 7 from line (|3.15|l . we a/so /lave 
the following space-time estimates for the interaction of the charge Q and the field 
quantity ^ : 

(7-16) ||r: +| -Vr 1+e K)^Q (1) -* (0) || i2[0 , to] < \q\ 2 ■E { sn ' e) [<P] , 

(u/ 7 )* (Q ( _i)-3(o)+ Q(o)-3 (0 )+ • *(_i)) ||i=[o,t ] 



»+#-« e 



(7.17) 

< M 2 -4 W) [3], 

(7-18) ||r|- 6 r!+ e K)iQ ( o ) -* ( _ 1 j|| i . [ o, ta] < |g| 2 -4 s,7,e) [3] • 



Remark 7.3. Since the proofs of these are a direct consequence of the estimates (|7.5|) 
and l|7.4[) - l|7.7bjl we simply make a few comments here. First of all, notice that the 
product estimates l|7.10fl and l|7.13|l must use the extra mixed space L 2 {L°°) norm 
on the right hand side of (|7.4|l . For example, if in the latter one were instead to 
use the pure L°° estimate for in conjunction with the space-time L 2 bound for 
3(_i), there would be an additional need for a factor of Tg. As we shall see in Sec- 
tion [51 such an extra convergence factor is not available in our application where 
estimate (|7.13fl comes up (although it is merely a notational convenience where 
(|7.1U|I arises). This is precisely the reason we have included the extra L 2 (L°°) 
norms. Thus, as far as our proof of Theorem 11.11 is concerned we never need to 
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make use of the pure L°° estimate for the quantity ^(i)- 

The second thing we call the readers attention to here is the fact that the esti- 
mates ()7.16|) - (|7.18|l will be used in their precise form for 2t < r. Notice that there 
are no extra factors of tZ which can be put in these estimates as was the case for 
()7. 9)1 - 1)7. 15)1 above. From a technical point of view (aside from convenience), this is 
the reason it is necessary to have the exact space-time energy norms on the right 
hand side of 1)7. 5|) . Without these, estimates of the form 1)7.16)1 - 1)7.181) would lead 
to logarithmic divergences which we see no other way of controlling. 

Before we finish, we take one last look at the estimates 1)7.9)1 - 1)7.18)) . Our pur- 
pose is to recast these in such a way that they conform more closely to how they 
will be applied to estimate the system (|1.17|) . In doing this, a surprisingly simple 
and elegant picture emerges of the underlying structure of the system 1)1.17(1 . This 
reenforces our point of view that the complex scalar field can be treated as if it 
were a tensorial quantity with roughly the same properties as the curvature F a p. 
We will put this structure to use through the so called "parity condition" . This is 
a numerical device used to keep track of weights and components in contractions 
such as the right hand side of 1)2.17(1 , which is typical of the kind of error terms we 
treat in the next two sections. It turns out that all of the estimates we will need 
can be reduced to a single streamlined form, the motivation for which will become 
a bit more clear through its use in the sequel. 

Proposition 7.4 (Abstract parity form of the weighted bilinear I? estimates.). 
Let ^ to be a given fixed time, and let < 7, e, s be parameters chosen so that 
2e ^ s — i. Recall the definition of the auxiliary weight u> 7j£ on line ((3.16)1 . and 
define the parity weights u> 7:C (a) via the formulas: 

(7.19) 

_ I 11 „ _ 1 

lU 7 ,e(l) = T e ■ {W)3,e , W 7 ,e(0) = T Q 2 • {wtf,e Wy, e (-l) = T Q S e • (iu) 7 , e ■ 

Also define the parity optical weights: 

(7.20) T(l) = T_ , T( ) = T(_i) = T+ . 

Then the following abstract bilinear L 2 estimate holds for the field quantities ty, $: 

(7.21) || Wl , e (a + b + c + d) T (a) r (b) * (c) ■ $ (d) \\ 2 L 2 [0M < 

4 s ^ e) (0,t )[vl/].|||$|||^ [oto](s!7ie) . 

whenever the condition —1 ^a + b + c + d^ 1 holds. One also has the analogous 
estimate for the interaction of Q and <J> under the same condition on a, 6, c, d: 

(7.22) || t*^ w 1>e (a + b + c + d) T (a) r (b) Q (c) ■ |||2 [0)to] < 

m 2 -4 s,7,£ \o,t m ■ 

Proof of the estimates ()7.21l) - ()7.22)) . We concentrate on the first estimate 1)7. 21)1 . 
the second being similar and much easier. First of all, we claim that this can be 
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reduced to the more restricted special case: 

(7.23) || r s + +i w 7 , e (b + c + d) r (b) * (c) • $ (d) ||| 2[0ito] < 

4 w) (0,< )[*]-|||$|||^ [0 , fo](s , 7 , e) , 

where now we impose the two conditions: 

(7.24) -1 ^c + d < 1, -1 < b + c + d ^ 1. 

Notice that (|7.23|) implies (|7.21l) with these extra conditions enforced because in 
that case one simply has the bound: 

w 7 , e (a + b + c + d) ■ r+ ■ T( a) < w^^(b + c + d) , 

as long as we are not in the case where both a — 1 and b + c + d — — 1. If this 
happens, we see from the previous conditions that b = —1,0 so we can use the 
bound: 

w lte (a + b + c + d) -T {a) T( b) < w 1 . £ (b + c + d+l) • r (h+1) , 
which holds in this case. Replacing now b' = b + 1 we are back to the case of 

To derive the estimate <|7.21|) in the case where c + d = —1 and b = — 1, notice 
that one must then have a = 1. In this case the bound rjT ^ 1 reduces things 
to H7.23fl with c + d = — 1 and 6 = 0. The case c + d = 1, 6 = 1, a = —1 can be 
treated similarly. 

To derive (|7.21() in the case where c = d = — 1, we use the symbolic inequality: 

t+Vm^-D < *(0) . 
which puts us in a position where we may again apply l|7.23|l under the conditions 

To deal with l|7.21|l in the extreme case where c = d = 1 , we simply use the fact 
that while it is not listed in the estimates (|7.9|) (|7.15|) . the product • is so 
favorable that it satisfies the space-time L 2 estimate (|7.21f) for any of the weights 
io 7)€ (a). 

It now remains to prove l|7.23|l under the conditions (|7.24() . This involves a 
simple case by case analysis split along the value of c + d: 

Case: c + d = 1. In this case we can either have b = — 1 or b = 0. In cither case 
the T( h ) weight is the same, so it suffices to consider the case which maximizes the 
weight w 7i e(6 + c + d). In this case we are dealing with w 7j£ (l). Substituting this 
into estimate (|7.23|) . and using the bound: 

S~\- — S-\- — -\-€. — 

T + 2 W 7 , e (l)T (0) = T + 2 TZ £ (w)l e , 

which holds due to the condition 2e ^ s — i, we see that the desired result follows 
from estimates H7.9f) - (|7.10|l above. 
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Case: c + d = 0. In this case, we can have all choices — 1 ^ b ^ 1. It suffices to 
consider the one which maximizes the product w 7j£ (6) • tvm which is easily seen to 
be b = 0. In this case we compute the total weight in estimate ()7.23(l to be: 



™ 7>e (Q)T (0) = T S + +1 Tl{ W )l 



< T S+1 



T_W) 7 ,e , 



where the last two inequalities follow from the condition 2e ^ s — | . We are now 
in a position to directly apply estimates l|7.11|l (|7.13|l . 



Case: c+d=-l. In this case we can have either b — 1 or b = so it suffices to 
consider the one which maximizes the product u> 7 . e (6 — 1) • t^). This turns out to 
depend on the value of \ < s ^ 1. If it is the case that I + f < s, then the choice 
6=1 maximizes the product in which case we are dealing with the total weight: 

S+- 3 1 

T + 2 W 7 , e (0)T (1) = T+Tl(w)j, e , 

T + T - Wr y,e ■ 

Notice that the inequality on the last line follows simply from the condition \ < s 
and the fact that we may assume e ^ ^. 

If on the other hand it is the case that s ^ | + § , then the choice b — maximizes 
the product u> 7i£ (& — 1) • Tn,) in which case we are dealing with the total weight: 

t + 2 u; 7!£ (-l)r (0) = tI t_ (w)l e , 

where, as before, the inequality is guaranteed by the condition 2e ^ s — i. Sub- 
stituting this last line into estimate i|7.23[l we can then directly apply estimates 
fFTty - QFTty . This completes our proof of (|7~2T|) . □ 



8. Differentiating the Field Equations I: Error Estimates for the 

CURVATURE F a p 

We are now ready to begin in earnest our proof of Theorem 11.11 We assume 
that we have fixed some level of regularity for the problem, say k derivatives with 
2 ^ k. The result will be demonstrated through a bootstrapping argument on the 
energy (|7.3|) . Recalling the definition of the norms for the initial data on the left 
hand side of l|1.25fl . we define: 

(8.1) E[ s ^(0)[F^} = ||£^|| 2 fffc , s+ , (R3) + \\H\\ 2 Hh , B+HR3) + 

|| D(f>0 Hjjfc.s+TpfS) + || 00 |!#fc, S + 7(R3) ■ 

We will now show that: 
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Theorem 8.1 (Bootstrapping form of Theorem 11.1(1 . Let k be a given level of 
regularity, and assume that we are given parameters < s, 7, e with the properties 
that s ^ 1, and Ae ^ s— \, ands + j < |. Let < T be a given time parameter, and 
let \1/ denote the totality of components of the system l|1.17fl as defined by (|7.2() with 
the associated k th level energy content i|7.3[l . Let E^f' 1 ^ (0) [F, (/)] denote the initial 
k th level energy content in the initial data (|1.18|) . Then there exists a constant 
1 *S Ck, s ,-y,e which depend only on k, s, e, 7, but not on T, such that if one first 
assumes that both: 

(8.2) 2#"> (0)[F,d>] , £^< £) (0,T)[*] < C£„, e , 

then the following nonlinear estimate also holds: 



(8.3) 2#"' 6 >(0,T)[*] < C M>7 , e 4^(0)[i^]+E [4''^(0,T)[*] 




In particular, we see from Theorem 18. II that the assumptions: 
together imply that: 

^) (0>T)[$]< l C -i 7ie . 

This, combined with the usual local existence theorem for the system (|1.17f) . and 
keeping in mind the main L°° estimate (|7.6f) (assuming that Cfc )Si7)e is chosen so 
large that the condition l(8.2|l allows one to absorb the extra L°° terms on the right 
hand side of this estimate) implies the claim of Theorem 11.11 Therefore, for the 
remainder of this paper we will assume that the solution exists up to time T and 
we will concentrate on proving the non-linear a-priori estimate 1)8. 3p . 

In this section we concentrate on proving (|8.3[) for the portion of ^Sf which con- 
tains the curvature. By the energy estimate l|3.48fl . this boils down to estimating 
the differentiated current vector C T X J at time t = 0, as well as over the space-time 
slab [0, T]. We do this separately. For the latter it suffices to be able to prove the 
following bounds: 



Proposition 8.2 (Space-time error bounds for the curvature F a p). Let k be a 

given level of regularity, and assume that we are given parameters < s,7, e with 
the properties that s ^ 1, and 2e ^ s — and s + 7 < | . Let < T be a given time 
parameter, and let: 

J = %((f)D4>) , 

be the current vector for the system l|1.17f) • Recall the current vector norm l|3.4Uf) • 
Then there exists a constant 1 ^ Ck, s ,-y,e depending only on k, s, e,7, such that if: 

(8-4) Ei s ^\0,Tm < C£„, e , 

then one has the following weighted space-time estimate: 

(8-5) £ |||^J|||£ 2[ o,T ](s , 7 ,e) < E [4 W) (0,T)[*]]' . 

|Z|<k 1=2 
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Proof of estimate Ij8.5(l . In light of the abstract L°° estimate (|7.6I) and the assump- 
tion H8.4f> with the constant C^l e chosen small enough that the right hand side 
of H7.6fl containing the extra L°° norm can be absorbed into the left hand side, it 
suffices to be able to show that: 

(8.6) 2J III £x J IHi2[0,T](*,7,e) ~ 
\I\^k 

^"'^(O^m • lll^*- 2 « r ]llioo I o, 31( .,- r , c) lll>C*- 2 « r llli«- [ o,r K -,^o) , 

where ||| £. x J |||i2[o,T](s,7,e) is the notation from line l|3.40|l above. Expanding out 
the norm on the left hand side of this estimate, and using the weight notation (|7.19() 
we see that: 

^ III C-xJ llli 2 [0,T]( S ,7,e) = E (ll T + +2 i^X J) L \\l 2 [0,T] 

|J|<fc \I\^k 
XEL XeL 

+ ||< + ^ 7 , e (0) £} x n\^,T\ + K + ^ 7) e(-l) (C x J)l\\l* [0 , T] ) , 

where we are using the notation \C X /| 2 = S AB (C X J)a • {C x J)b- Therefore, using 
the abstract parity estimates (|7.21|l in the form of estimate (|7.23l) as well as estimate 
it suffices to prove the following symbolic bounds: 

(8.7) 

E k £ ^i 

ma 

veL 

< 2J T ( Q ) ( £ L*(f>) + £ L<9(6)) ' *(c) ' (l + III ^l~ 2 ^ \\\L°°[a,T]{s,~i,e)) > 

l+m—k—1 
or Z=0,m=fc 

and a+b+c=i 
(8.8) 

E \^xJ)l\ 

\i\^k 

xgl 

% E T W (4*(6)+4Q(6))^r*(c)-(l + lll4" 2 *IIU»[0 > T]( S ,7, e )) 

l-\-m— k— 1 
Or l=0,m=k 

and a+6+c=-l 

(8.9) 

E 1^/1 

|Z|<fc 

veL 

~ E T ( Q ) ( £ L*(f>) + ^L^W) ' *(c) ' (! + III £l7 2 ^ IIU°°[0,T]( S ,7,e)) ■ 

l+m— k— 1 
or Z=0,m=fc 

and a+b+c=o 

We will only prove the bounds l|8.7|> - (|8.9[) in the extended exterior region t < 2r. 
These bounds in the complimentary region r < \t follow from similar reasoning 
and are in fact much simpler because the weights are all identical there. 
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Our first step is the following simple inductive calculation of the Lie derivative 
C x J which is based on repeated use of the formula (|2.15|1 : 

(8.10) 

= J- E + (d^-l Xi Ki dA , 

1 K 1+ K 2 =I \ " / V / 

= r E ( D xl^ ■ D(D%d>)) + (pzy ■ D(D^<f>)) 

+ r E (d x ^-[c X2 k \d]A + (d k x i^[c Xi k \d^) , 

1 K 1 + K 2 =I \ - / \ / 

= ,4 + 5. 



We now compute each of the A and -B terms on the right hand side of the above 
identities separately. Each of these terms can be seen as real valued two forms, and 
we denote their components by A a and B a respectively. 

To compute the A term, notice that since the sum is symmetric in the K\ and K2 
multiindices, we may replace the Dl derivative by -Dl(t-) in the A^ component. 
Doing this, and putting absolute values around the different components of A while 
using the condition |/| $J k we have the bounds: 



|Ki|+|ifa|<fc-l 



(8.11) \A L \ < J2 T + \—\-\-D L (rD«<P)\ 

\K\ = k + 

]T \DYD%<t>\-\±.DUrD%4>) 

-[JCaKfc— 1 
,X 2 ,YGL 

■12) \Al\ < £ r + |^-| • \D L {D^4>)\ 

\D Y (D^)\-\D L (D^J)\ 

-\K2Kk-1 
,x 2 ,YeL 

■13) 141 < E r + \^\-\P(d^)\ 



\K\ = k 



\K 1 \ + \K 2 Kk-l 
Xi,X 2 ,YeL 



\K\=k 



\ d y(Dx\4>)\ ■ \p(Dxl 



x lt x 2 ,Yeh 



For the first term in each of the above sums, it suffices to merely recall the designa- 
tions l|7.1 fl -l|7.2 |) and the definition of the parity weights T( ) to achieve the bounds 
(|8.7(l - (|8.9(l . To achieve these bounds for the second term on each line l|8.11(l - l|8.13(l 
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above, it suffices to show the bounds: 

Y \dy(d%<i>)\ < Y r w • 

IKil^i a+b=0 

This last line follows from a straight forward application of the bounds: 

(8.14) \X L \ , \X A \ < t + , \X±\ < r_ , l£L, 

which follows from an inspection of the identities (|2.32|) . together with the desig- 
nations H7.1jl - I|7.2jl and the weight definition l|7.20|l . 

We now move on to dealing with the B portion of identity (|8.1()(l . This boils 
down to computing the commutator actions [C x ,D]<p. By a repeated use of the 
formula Ij2.14|l and the fact that L is a Lie algebra (so all its commutators involve 
constant coefficient sums), we easily have the bounds: 

(8-15) \B a \ < Y \(^x[F) a \ ■ \D%4>\ ■ \Dfj\ . 

\K 1 \ + \K 2 \ + \K a \^\I\-l 
Xi,A- 2 ,X 3 ,VeL 

We now use the fact that either \K%\ ^ k — 2 or \K^\ ^ k — 2 which comes from the 
restriction given above on their sum, to employ the L°° estimate for D x cf) contained 
in l|7.6J) to bound the product of the second and third factor in this last sum as 
follows: 

i^i-i^i < £r x{|K2MK3l} *(o)-ii4- 2 *ii^[o,T ](W ). 

Therefore, to achieve the right hand side of l|8.7J) - l|8.9J) for the B term it suffices to 
show the bounds: 



(8-16) Yl \(iY^xF) L \ < £ t (o) (4* (6 )+4Q (6) ) , 

\K\^l a+6=l 
X.Y&u 

(8.17) Y \^c%f) l \< y > 

\K\^l a+b=-l 
X.Y&u 

(8.18) Y Y\^ c xF)a\< Y n«)(4%)+4Q(6)) ■ 

\K\^l A a+b=0 

v.vgl 

Splitting CxF — C X F+C X F we see that it suffices to do this calculation for the F 
portion of things because the same computation for the charge field F is identical. 
Since this is an abstract counting argument for an arbitrary two form we can drop 
the Lie derivatives and the tilde notation. Making now the identifications: 

(8.19) L & (1) , L ^ (-1) , A,B (0) , 

we see that in the groupings (|7.1|I - H7.2|I the component F a p is put in the set that 
corresponds to the sum of a and (3 (according to the numerical designations (|8.19|) 1 
as they range over the null frame {L, L, ca}- Using this observation in conjunction 
with the bound H8.14|l and the weight definitions H7.2()|l . it is seen that the estimates 
(|8.16[1 - H8.18JI are an immediate consequence of translating the contractions X®F a p, 
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computed in a null frame, into the more simple numerical parity sum form. This 
completes the proof of the estimate l|8.5l) . □ 

To complete the proof of the bootstrapping estimate (|8.3I) for the electro-magnetic 
field F a p , it suffices to be able to bound the initial data type norms for the current 
vector J given on the right hand side of (|3.48|) in terms of the energy (|8.1|l . To do 
this it is enough to show: 

Proposition 8.3 (Initial data bounds for the system l|1.17[l )■ Let 2 ^ k be a given 
level of regularity, and assume that we are given parameters < 7, and i < s, 
and s + 7 < I . Let (F, <j)) be a solution to the system (I1.17|) with initial data 
(E, H, 4>q, 4>q). Then there exists a constant Ck, So such that if one first assumes the 
smallness condition: 

(8-20) 4 S ' 7) (O)[F,0] < C£„, e , 

then one has the initial bounds: 

(8.21) 

E II (i + s+7+|/| v£ mo) ||% + E II C 1 + r) s+7+|7|+ X x J(o) || 



2 

LI 



|J|<fc x |/|<fc-l 



< 



4*' 7) (o)[f,0] 



2 



Proof of the estimate H8.21fi . We begin by bounding the first term on the left hand 
side of (|8.21|) . At time t = we directly compute that in terms of the initial data: 

(8.22) E l v ' J ol ^ E \R Kl <h\-\R K 'io\ ■ 

\I\=l \K 1 \ + \K 2 \=l 

This, combined with the two L°° estimates: 

(8.23) (l + r) s +^|0 o | < E ll(l + ^r +7+m ^^o||^(R3) , 

(8.24) (l + r)Wi|0 o | < E II i 1 + r^+Wfa ||l»(r») , 

|Z|<2 

which in turn follow for the estimates H6.5[l - (|6.6|) . achieves the bound l|8.21[) for this 
first term through a use of the L 2 ■ L 3 Holder inequality and directly integrating 
the resulting L 3 estimate using H8.23H - H8.24fl an d the fact that \ < s + 7. 

To estimate the second term on the left hand side of (|8.21|) . we first expand 
the derivatives Vt, x J as we did in line (|8.22() above. Doing this and using the L°° 
estimate: 
(8.25) 

E (l+r) 1+lKl \DU < e ni+rr^^D^D^mm, 

|Jf|<fe-i i^|/i|+|/ 2 |^fc+i 
xe{d a } ^e{0i},xe{9 Q } 

which follows from the estimate 1)6. 5[) applied to functions of the form (l+r)' K ' D^ip 
and the fact that 1 < | + s + 7 as well as | < s + 7, we are reduced to having to 
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demonstrate the general estimate: 

(8.26) H( 1 + s+ ^ 1+|/| ^^lli 2(R 3) < Ei s '^(0)[F], 
i^|/|^fe+i 

xe{d a } 

holds, provided that one first assumes the smallness condition H8.2()(l as well as the 
condition 2 ^ k. 

It is clear that the work in showing l|8.26f) is simply a matter of eliminating the 
D t derivative in favor of Di derivatives. This in turn boils down to controlling the 
commutators of the form [Dt,!) 1 ] because our goal it to move the D t operator as 
far to the right as possible and then to use the field equation i|1.2[l . While this 
overall strategy is simple to describe, the implementation becomes a bit tedious 
because the commutator formula (|1.5(l introduces the field electric Ei which must 
then be differentiated many times. 

The simplest wa to deal with all of this seems to be to fix a value for k and then 
to induct on the number of D t derivatives on the left hand side of this estimate. 
If this number is zero then the claim is obvious, so we may now fix the number at 
some value 1 ^ / ^ k + 1. Our goal is to prove the bound: 

(8.27) \\(l + rr^- 1+ ^D^D t D^ct>f L2m < E^\o)[F} , 

where |_ftTi| + \K 2 \ ^ |/| — 1 and X £ {d a } while Y £ {di}. By induction, we may 
assume that the same estimate holds with the D t derivative replaced by one of Di . 
In this regard it will be useful for us to employ the following notation: If if is a 
multiindex, then we denote by \K°\ the number of time derivatives in the operator 
(note that we are only working with X £ {d a } now). With this convention, we 
have that \K®\ ^ I — 1 in (|8.27|) above. Now, computing some multiple commutators 
we have the point-wise bound: 

(8.28) \D^D t D**4>\ < \D^D«*D t <f>\ 

+ Yl \D^D^[D t ,D z }D^\ . 

\h\ + \l2\ = \K 2 \~l 

xe{d a } 

Vi,y 2 ,Z6{.9i} 

We estimate the second term in the above expression first. Substituting this into 
the left hand side of l|8.27|l and expanding out the single commutator in terms of 
the electric field Ei , distributing the left hand derivatives according to the covariant 
Leibnitz rule, using the field equations (| 1 . 1 Ul) to remove time derivatives from E, 
and splitting result along the Hodge decomposition E — E d f + E c ^ we have the 
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bound: 
(8.29) 

\K 2 \+\Ii\+\h\<k-l 
X£{d a } 
Yi,Y2,ZG{a 4 } 

S E (( II 0- + ry+r+ 1+ W+M(V%E*) ■ (D%4>) h> m 

|J"i|+|i 2 |<fe-i 
\i°\^i-i 
xe{d a } 

+ || (1 + ^^i+l^l+I^I^F) • (D%<t>) || L2(K 3) ) 
+ ||(l + r)-+T +1+ ^W I »l(V^^).(i^0)|U.CR») ) 

+ E 

|ii|+|i 3 |^fc 

= + 5 + C . 

We now bound each of the three terms on the right hand side above separately. For 
the A term we are done by induction after a simple use of the L°° estimate l|8.25|) . 

To bound the B term, we add and subtract off the term V x ■ x + (r — 2) from 
the curl- free component E c f . Then using the estimate H3.53jl in conjunction with 
the bounds we have already established for the first term on the left hand side of 
(|8.21[) . and by making another use of the estimate (|8.25jl this we are reduced to 
showing that: 

\I 2 \^k-l 

\i°Ki-i 

But this last estimate follows from our inductive hypothesis unless I2 = in the 
sum. In this latter case we can use the Poincare estimate (16 .7|) to make our claim. 

Our next task here is to prove the (inductive) bound (|8.21|) for the C term on 
the right hand side of line l|8.29|l above. As before, we use the estimate (|8.25ll to 
deal with the factor involving (1 + r) 1+ ^ l2 ^D-^(f>. Doing this and making use of 
the continuity equation (|1.16() and the smallness condition (|8.20l) . we arrive at the 
bound (as long as < I — 1): 

C< E ll(l + r-) wl+|7l| V t V(0)IU=( R3) ■ 

\hKk-2 

We can now reduced things to the point where we can go back to the paragraph 
containing line l|8.25f) above and repeat the entire argument up to this point which 



II (1 + rr+^+l^+I^V^J) . (D 1 ^) || i2(R3) , 

-2 
-1 
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allows one to bound the right hand side of this last expression via induction. 

To finish here, we only need to deal with the first term on the left hand side of 
(I8.28|l . First of all, if it is the case that \K^\ = then we are automatically done. 
If however 1 ^ \K\\ then we simply repeat the steps outlined above for the second 
term on the right hand side ((8.28(1 , but this time applied to the second term on the 
right hand side of the bound: 

(8.30) \D^D t D^D t (j>\ < \D^D^D 2 t <j>\ 

+ J2 \DfD^[D t ,D z }D^D t( f>\ , 

\Il\ + \l2\ = \K 2 \-l 

xe{d a } 

where this time |-ftfi| = \K\\ — 1 and K\ = 1 — 2. Notice that these numerical 
conditions guarantee that no more than a total of k — 1 derivatives altogether, and 
no more than I — 1 of the D t derivatives in particular, can fall on (j) in this second 
term. Therefore, we are again left with estimating terms which are similar to those 
appearing on the right hand side of 1(8.29(1 above. 

Finally, to deal with the first term on the right hand side of 1(8. 30(1 we simply use 
the wave equation 1(1.2(1 which trades the derivatives D\ for the covariant Laplacian 
D l Di. Having reduced the total number of time derivatives in this last expression 
we are done by induction. This completes the proof of the estimate (18.21(1 . □ 



9. Differentiating the Equations II: Error Estimates for the 

COMPLEX SCALAR FIELD (j) 

To wrap things up, we need to prove the boot-strapping estimate ((8.3(1 for the 4> 
portion of things. Using the energy estimate 1(4.19(1 for complex scalar fields, and 
using the assumption that C^l e in 1(8.4(1 is chosen so small that the L°° norm term 
on the right hand side of 1(4.19(1 can be absorbed into the right hand side of that 
inequality, and using the estimates 1(6.7(1 and 1(8.26(1 above to deal with the terms 
involving the initial data we see that it is enough to prove the following L 2 estimate 
for the commutator of covariant derivatives and the complex D'Lambertian: 

Proposition 9.1 (Error bounds for the complex scalar field 4>). Let k be a given 
level of regularity, and assume that we are given parameters < s,7, e with the 
properties that s ^ 1, and 2e ^ s — \, and s + 7 < |. Let < T be a given time 
value. Then there exists a constant 1 ^ C\ S!7ie depending only on k, s, e, 7, such 
that if: 

(9.1) E^\0,Tm < C£„, e , 

then one has the following weighted space-time estimate for commutators: 

(9.2) £ ||<r^Hi [□ C ,Z^]0|| L2[O , T] < J2 [^ 7 ' £) (0,r)[*]]' , 
|j|<fc 1=2 

VGL 
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where u> 7i£ is the weight function from line 1(3.16(1 above. 

Proof of estimate 1(9. 2|) . As in the previous section, we may assume that C^ s e 
is chosen small enough that the right hand side of 1)7.6(1 containing the extra L°° 
norm can be absorbed into the left hand side. Taking this into account, it suffices 
to be able to show: 

(9.3) Yl K T -(™)i [□ C ,^x]0II^[o,t] < 

J|^fc 

Our first step in this process is to expand the commutator [□ , D^-] . This is done 
through the following multiindex identity: 

(9.4) P C ,Dx]<l> = E D%P°,D%]D%<t>. 

K 1 + K 2 +K 3 =I 
\K 2 \=1 

x 1 ,x 2 ,x 3 eh 

Now, using the single derivative commutator formula (|2.17|) in conjunction with a 
repeated application the two identities: 

D x (Y a F a(3 D l3 c/ ) ) = [X,Y) a F af3 DP<i> 

+ Y a £ x F aP D^ + Y a F a0 D (D X (j ) ) 

- Y a F aP ( X V^ZV + \/~^\Y a FapF^ 3 X 1 ■ <j> , 

and: 

D x (V ff (Y a F a/3 ) ■ <f) = VP{[X,Y} a F a p) ■ </> + V f3 (Y a £ x F a p) ■ <j> 

+ V (Y a F a/3 ) ■ D x <t> - W7r^V 7 (y Q F Q/3 ) • , 

which in turn follow from the geometric formulas 1(2. 4JI . 1(2.11(1 . 12.16(1 . 1(2.14(1 . and 
(12.9(1 (respectively), we arrive at the following point-wise bound for the left hand 
side of 1(9.4(1 above: 

(9.5) J2 l[ nC >^l> 

\I\<k 

< £ (1^(^)1-1^1 

\Ki\+\Kz\4,k-l 

Xi,x 2 ,Yeh 

+ \2Y«{C«\F a p) ■ D^Dfj) V a (Yl) ■ {C K x \F a(j ) ■ 

+ J2 \m^xlF^)\ ■ \Y 2 \C^)\ ■ \D*14>\ , 

|Jfi| + |K 2 | + |if s |<fc-2 
X 1 ,X 2 ,X 3 ,Yi,Y 2 eh 

= A + B + C . 

Our goal in now to prove the bound 1(9.3(1 for each of these three terms. We do this 
separately and in order. As with the proof of 1(8.7(1 - 1(8.9(1 above, we will only carry 
out our calculations in the region t < 2r as the calculations necessary to handle 
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the region r < \t are essentially trivial because all components of (F, D<p) can be 
treated on an equal footing there. 

Our first step is to deal with the A term above. Taking in to account the abstract 
parity estimates H7.21fl - H7.22fl . our task is reduced to showing the symbolic bound: 

(9.6) A < 

E T + lr (a) T (b) ( £ L*(c) + £l<2(c))-£l *(<*)•(! + III C-l7 2 ^ IIU°°[0,T]( s , 7 ,e)) • 

l+m— k— 1 
a+b+c+d=0 

Expanding out the contraction with respect to (3 in the first factor of the A term, 
and using the identifications (|8.19fl we have the straight forward bound: 

(9-7) \Y P C%\M < E T M \£xA>>)\ ■ 

a+b=0 

To bound the C^J^ term itself, we can simply use lines (|8.7fl - (|8.9fl above. How- 
ever, it will be necessary for us to use the following refinement of those bounds, 
which can easily be checked by recalling lines H8.llfl - H8.15fl : 

E \ C xl J (b)\ ~ 

E ^ 1t (o) (4:%+4 1 Qm)' E i^i ■(! + hi 4~ 2 * iiiL»[o,n ( ., 7 , e )) ■ 

c+d=b |Ki|<ia 
h+h=l ViGL 

Expanding this into H9.7fl above, and tacking on the extra factor of D^fi while 
using the bound: 

E k^i • \d^j\ < £r x{Wn So) ■ ni4- 2 *iiiL^o,T ](s , 7 , e ) , 

l2-\-m^k — l 

x u x 2 eh 

noting that by design we have l± + max{i 2 ,m} ^ k — 1, we see that we have 
achieved H9.6fl . Notice that here we are making (crucial) use of the special bound 
|Z>^0| < tTJ , when |/| ^ k — 1, contained in the norm H7.4fi . 

We now move on to bounding the B term in line H9.5fi above. Here we will prove 
the symbolic bound: 

(9.8) B< E T (a) (4*W+4Q( b ))-^*(c) • 

a+b+c=0 
Z+m^fc-1 

Through an application of estimate (|7.21|> . this will prove the bound H9.3fl for the 
B portion of things. To simplify matters, it suffices to prove the above bound 
with I = m = 0. The same estimate containing derivatives is simply a matter of 
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notation. With this in mind, it suffices to show that both: 

(9-9) \Y a {F aP )-^{rcj>)\ < £ r (a) (Q {b) + • * (c) , 

a+b+c=0 

(9.10) 

|(V(r).^-V Q (^)).^|.|0| < £ r (a) (Q (fe) + *«,))• * (c) . 

a+b+c=0 

Notice that the proof (|9.9() is a simple matter of applying the definitions l|7.1[) . 
(I7.2()|l , and the weight bounds (|8.14|) . The key observation here is that since this is 
a full contraction, its parity weight must be zero. Also, note that the term ^D^^cj)) 
can be treated as a ^(-i) term on account of the symbolic bounds (i.e. any bound 
for the right hand side is satisfied by the left hand side): 

*(o) < ■ 

Therefore, to prove l|9. 8fl it now remains to show the bounds (|9.1U|) . For the most 
part this is a simple matter of noticing that from the formulas (|2.34|) (or by ho- 
mogeneity!), each component of the covariant derivatives on the left hand side of 
(|9.10|l satisfies the bounds: 

|-V Q (r) -X ] + |V Q A^| < 1 . 
r 

Thus, if the parity weight of F a p is (0) or (1) we can pass to the right hand side of 
(I9.10|l after multiplication by the bounded factor r + • r -1 . The only case where this 
general procedure does not work is when F a p has weight (—1). In this case, to pass 
to the right hand side of (|9.1Q(I we must pick up an extra factor of tq. This indeed 
turns out to be the case, and is perhaps the most striking structural property of 
the commutator (|2.17|) as was pointed out in (see p. 226 of that work) in a 
somewhat different form. What we need to show is that: 

(9.11) l-V^fr) -X A - V L {X A ) + V A {X L )\ < to, 

for each X £ L in the region t < 2r. This follows from a direct use of the identities 
(I2.34|l . Notice that for the case X £ {d a ,S}, the above sum either vanishes or 
is of the order 0(r~ r ). Thus, the main thing to check is that l|9.11|l holds when 
X £ {Qij, f^io}- This in turn follows easily from the formulas: 

To finish our proof of estimate l|9.3|) , we only need to show this bound for the 
C term on line 1)9. 5|) above. Writing F = F + F and expanding out this term into 



-— n n 

2r 



2 Wi 



v A (^) 

V A (ttf ) 



2r ^ 



2^ 
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parity notation, we have the symbolic bound: 
C < 

^ r + 1 T (a)T(b) (^L*(c) + £LQ(c)) ■ £™ *(<J) ' III -^l" 2 * IIU~[0,T]( s , 7 ,e) 
Z+m=fc-2 

+ X! T + T ( a ) T w £lQ(c) • £Z l Q(d) ■ £iT 2 *(o) > 

J+m=fc-2 
a+b+c+d=0 

= Ci + C 2 . 

To prove the bound l|9.3l) for the term C\ is a simple matter of referring to the 
abstract parity estimate (|7.21(l . Thus, we are left with bounding the term C^- This 
reduces directly to a special case of the estimate l|7.22|) after an application of the 
L°° bound (|7.7|) which easily implies that their exists a cutoff function \ + such 
that: 

T + T (a) T (b) £lQ(c) ■ £hQ(d) < t+ ■ |g| ^ ^ J — — ■ \q\ . 

a+b+c+d=0 

Notice that the right hand side of this last line behaves the same in estimate (|7.22|) 
as the expression: 

T + 1 T (0) r (0) £jT 2 Q(0) ■ |||£]L~ 2 *IIU o °[0 ) T]( s , 7 ,e) ■ 

Multiplying this by the factor £^ _2v I / (o) an d then substituting the result into left 
hand side of l|9.3[) . we are done by an application of the parity estimate l|7.21fl . This 
completes the proof of (|9.2I) . □ 



10. Appendix 

This appendix contains the proofs of several more or less standard weighted 
Sobolev type estimates which are used at various places in the main paper. In 
particular in Sections |2I and EHEl We make no claim of originality here, rather our 
purpose is to have things stated in the precise form in which we find them conve- 
nient to use in our work. The first such estimate is a simple weighted version of 
the usual L 2 L 6 embedding: 

Lemma 10.1. Let p be a real valued test function on K 3 , and set: 
(10.1) q = / p dx , 



the average of ip. Let 5 be given such that i < 6 < 1. Then the following weighted 
integral inequality holds: 

(10.2) / R3 ,-|V (I, + £)f^ < U^lllf 

In the above estimate, the implicit constants depend on S. 
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Remark 10.2. In effect, what estimate ()10.2(l shows is that is possible to commute 
the weight r s past the Riesz operator VA _1 so long as one first subtracts off the 
leading order term in the asymptotic expansion of A -1 <^>. This is: 

(10.3) V,-^ 



'A^ Airr 2 ' 

That such a correction is necessary is clear from the range of weights we are consid- 
ering because for ^ < S a function with this asymptotic behavior at infinity cannot 
be in L? weighted by r s . The idea behind p0.2(l is, of course, that once the right 
hand side of 1)10. ?>\ is subtracted off from VA^'ip the resulting function will decay 
enough to be in r s weighted 1? as long as the appropriately weighted Ls norm of ip 
is bounded. However, there is a limit to how much weight one can expect to apply 
this way. This is because one can further expand (|10.3|) to include higher order 
terms: 



(io.4) V ^~- V >(^+E^r + {M9her} 

where: 



Ik = / Vk ¥>{y) dy , 
Jm. 3 

are the higher moments of ip. Without the additional vanishing of these other 
quantities, it is clear that one cannot put r^VA^ 1 ^ in L 2 for | ^ 5. (|10.4|) 
is sometimes referred to as the multipole expansion of V^y. Since we are only 
interested in the decay of the initial data (|1.18|l modulo charge of the order no 
greater than 0(r -3 ), the first term in I|1U.4|I will be the only one which concerns 
us in this work. 



Proof of (|10.2I) . Our first step here is to simply integrate by parts several times on 
the left hand side of I|1U.2|I . This yields the identity: 

(10.5) (L.H.S.)ECU = 5(26+1) ( r 2S ~ 2 (4v> + dx 

To bound the integral in the first term on the right hand side of (|10.5|) above, we 
use the definition UlU.lfl to compute: 

/ ^ {^+-T-) 2 dx < [ r 2 «- 2 ([ |_J— --i-l |^(y)| dy] dx. 
Jm? V A 4?rry J R3 \J R3 \x - y\ \y\ 1 J 

We now split cases depending on whether |<5^1orl<5<|. In the first case 
we use the bound: 
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This reduces our work to estimating the two integrals: 

i / r 1 \ 2 1 

y\ 5 My)\ dy) dx , - < S < 1 



A \Jw> \x-y\ J 2 

B =f r^ii r^-is \y\ 5 \^y)\dy) dx, \ < s < 1 

In the second case, we simply use the bound: 

i 1 1 i / \y\ 

5> 



\x-y\ \yy \x\\x-y\ 
This reduces us to bounding the integral: 

T^=2s(l 1 \tts^t \v\ S W{y)\ dy] dx, 1<S<^. 

R \Jr3 \x-y\\y\ d J 2 

In all cases, the bound we wish to prove is: 

(10.6) A,B,C, < \\r B <p\\ Li . 

This last estimate follows from several different instances of the generalized frac- 
tional integration Lemma Tl U . 31 below . Notice that in the case of integral A we have 
(notation of Lemma [10.31) a = j3 = 1, and 7 = 0. In case of integral B we have 
a = 2 — 5, j3 = S, and 7 = 0. In the case of integral C above we have a = 2 — 6, 
j3 = 1, and 7 = 5—1. In all cases we have that p = p' = 2 and q = | and 
a + /3 + 7 = 2so the scaling condition (|10.10f> is satisfied. Also, note that in each 
case we have that a < § and 7 < i so the "gap" condition I|1U.11[1 is satisfied. This 
completes the proof of the bound: 



(10.7) 



r ^-2 (}_ v+JLY dx 
V A Attt 



It remains to bound the second term on the right hand side expression ()10.5f) 
above. By Holders inequality we have that: 



.,2(5 



d.r 



\\r S <P\\ L % 



A iirr J 



By the usual L 2 L 6 Sobolev embedding, the Leibnitz rule, and the triangle 
inequality, the second factor on the right hand side above can be bounded by: 



|| r s 



Airr J 



L 6 



< 



1 



-<P- 



47rr J 



\L 2 



r s V 



i 2 



Therefore, using the bound l|10.7fl above, we can estimate: 



(10.8) 



r 2 V 



-tp- 



< 



Aitr I 



dx 



IkVll 



L~> 



\\r 5 v\\ 



Anr 



)H 



L 2 



Adding estimates (|lf).7fl and l|10.8|l into the identity 1)10. 5|l . dividing through by 
the quantity || r s V(^<yJ + -£-) ||x,2 and resquaring, we have achieved the desired 
estimate (tTU^|) . □ 
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To complete the proof of l |10.2|l , we need to show the bound H10.6|) for the integrals 
A,B, and C above. This will be a consequence of the following generalization of 
the classical Hardy-Littlewood-Sobolev fractional integration theorem (see |llp: 

Lemma 10.3 (General fractional integration lemma). Let 1 < q < p < oo and 

^ a, [3, 7 be given parameters. Then the following integral estimate holds for 
(positive) test functions F and G on M 3 : 

(10 - 9) /,./,. GW WV^WW n ' )ixiy 5 '' G ''-'-'' f l'"- 

where the following conditions hold on the various indices: 

(10.10) 3-(a + /3 + 7 ) = 3(---), (scaling), 

q p 

(10.11) a < 3(1--) , 7 < 3(1--) , ("gap") . 

p q 

The implicit constant in the above estimate depends on p,q,a, /?, 7. 



Proof of (|10.9H . By the restricted weak type version of the Marcinkiewicz interpo- 
lation theorem, it suffices to prove that: 

1 J- 1 

\x\ a \x - y\P \yp 

for measurable sets E and F where 1 < q < p < 00 and the conditions (|1C). 1 Of) 
(II L) . hold. By the Riesz rearrangement inequality (see ^T]) this is further reduced 
to showing that: 

(10.12) f f ( X e(x)-\x\-<*)* —t—g (xF(y)-\yn)* dxdy < \E\7.\F\ k «, 

Jr3 Jm* f - y\ p 

where /* denotes the symmetric decreasing rearangement 14 of the function /. We 
now compute that: 



{xe(x)-\x\ a )*{r) = J X{ XE .\ x \-° >s y(r) ds , 

X mE nB { s-i)li) ir) dS 



= (4 7 r)- 1 ( X£ )*(r).r-i , 

and similarly for (xf(x) • 1^1 ~ 7 )*- We now plug these formulas into the right hand 
side of 1|1U.12H and apply the usual fractional integration theorem to yield: 

(10.13) f [ ( X e(x)T ■ \x\~ a ( XF (y)r ■ \yr dxdy 

Jr3 Je? I* - y\P 

< \\(XE(x)r-\x\- a \\ LSl -UxF(x)r-\x\^\\ L e 2 , 



14 Here we use a definition of /* which does not include the usual normalization factor of 4ff 
to avoid additional typesetting. This is xj;( r ) = ^ £ • Clearly this convention does not effect 

the use of rearrangements in the Riesz inequality. 
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where we choose 0i and 02 to be given by the expressions ^- = y + f and = i + ^ 
respectively. Notice that from the scaling condition (|10.10|) and some quick algebra 
we have that : 6-/3 — 3(^- + ^-), which is the scaling condition needed for I|1U.13|) 
to hold. Furthermore, notice that by the gap condition we have that both 

1 < 0i,02 < 00 which is also needed for (|10.13[1 to hold. The proof of (|10.12|) is 
now accomplished by simply computing: 

II (Xe(x))* ■ \x\~ a ll^'ej = (J^ r^dr) k , 

< |£?|W (3_Qei) 
= \E\ ■ 

Notice that the inequality follows because we automatically have that a0i < 3. An 
identical computation shows that: 

UxF(x)y-\xr\\i„ 2 < \f\ . 

This completes the proof of I|1U.12|) and therefore the proof of estimate (|1U.9|I . □ 

The second main set of estimates we prove here are a localized version of the so 
called global Sobolev inequalities. These were first utilized in [5] to prove global 
existence and asymptotic behavior for general non- linear wave equations. The ver- 
sions which we state here are sufficiently "atomic" to provide all of the L°° we 
will need in this paper. To discuss these, we shall employ the same notation as 
introduced in the beginning of Sectional In particular, the notion of a dyadic cone 
distance (CD) shell and the associated t + (J),t^(J) notation. 

Lemma 10.4 (Exterior global Sobolev estimate). In the exterior region 1 ^ t < 2r 
let J be a given dyadic CD shell, and let q be given such that 2 < q < 00. Then 
for test functions f the following estimates hold: 

(10.14) sup \ X f(t,r,Lu)\ < t~«(J) V || ^- x/(i,ry)|U«(s=) , 

where \ * s the cutoff on the dyadic CD shell J . Furthermore, in the case where 

2 ^ q < 4 we also have the estimate: 

(10.15) ||x/IU-(£«(s»)) < 

T~ + 2(h ~ k \j)Tzhj)(\\xf\\u> + \\r-d rX f\\ L , + ]T W^xfW^) ■ 

i<j 

Proof of estimates (|10.14fl - HlU.15fl . The proof of both of these is entirely standard 
as they are essentially just rescaled versions of the usual translation invariant 
Sobolev estimates. For example, (|ll).14f> is just the Sobolev theorem on spheres 
§ r . For the sake of completeness we give here a proof of <|10.15|) • 

The first step is to introduce a new set of variables which we denote by (r, to) and 
are supported in the unit dyadic annular region A — [|, 8] x S 2 (where S 2 without a 
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subscript denotes the unit sphere), and which satisfy the change of variable formula: 

(r-(J)-r + t,r- u) = (r,w r )elx§J. 

That is, we rescale the standard radial variable around the point r = t by the factor 
T-J and we rescale the angular variable by the factor of r. This procedure gives 
the change of volume formula: 

(10.16) drdujr = v{r)-(Fduj, v{r) ~ r+{J) - T-(J) , 

where the comparison on the right hand side above holds in the region x £ J . Using 
the analog of l|10.16|) for integration on a fixed sphere we have the equivalence: 

||x/IU~(z«(sa)) ~ t+(J) ■ II xfh^(LH^)) ■ 
Furthermore, a direct calculation using mm gives: 

T + (J)T?(J)-(\\ X f\\L* + ^2\\n ijX f\\^ + \\T-drXfWv) ~ 

i<j 

\\X.f\\L^(A) + \\drXf\\L*(A) + II fxf \\l*(A) ■ 

Using these last two expressions in tandem we see that the proof of (|10.15() reduces 
to being able to show that: 

(10-17) ||x/IU r( L« ( s=)) < WxfWm(A) ■ 

Notice that by design, the rescaled X f does not intersect the boundary dA. Using 
two charts on the interior of A, estimate (|10.17t follows from the general mixed 
Sobolev embedding on R 3 for test functions ip: 

(lo.is) IMUsm) ~ II ^ll^ 1 . 

where (x,y) elxl 2 and | — i — | < 1. This last estimate in turn follows from 
running the Sobolev lemma in the x and y variables separately to achieve: 

IMUczd ^ II W a ~' )+ W 2(i ~«Vll^ • 

where the "+" notation denotes an arbitrarily small positive constant which gets 
us around the case p — oo. (|10.18|l now follows from this last estimate and the 
symbol bounds: 

(i + |t| a )4(*-*)+(i + |^| a )<*-i) < (i + |&,| a + |&| a )i, 

whenever | — i — | < 1. We have now shown (| 1 () . 1 5|) . □ 

The second global Sobolev estimate we prove here is the analog of Lemma fl 0.41 
for the interior region r < ht. This is: 

Lemma 10.5 (Interior global Sobolev estimate). Fix 1 t, then in the interior 
region r < jt one has the following weighted estimate for compactly supported (in 
that region) test functions f: 

(10-19) \\f\\ L « < t- 3 ^-^(\\f\\ L * + Yl H*(/)IU")' 

xe{s,n 0i } 

whenever - — - < %. 

p q 3 
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Proof of estimate 1(10. 19j) • Note that this estimate is taking place in a ball of radius 
t. By simply rescaling the usual Sobolev estimate we have that: 

ll/IU. < t- 3 ^-^{\\f\\ LP + t \\vf\\ LP ) . 

Therefore, the desired result follows directly from the bound: 

* iiv/iu, < E n*(/)iu* . 

xe{s,n 0z } 

which in turn follows from the identities: 

tdi = to i0 -Xid t , 
d t = (t 2 - r 2 )- 1 (tS - x l n iQ ) . 

a 

Finally, we end this appendix with a characteristic version of the exterior es- 
timates (|10.14fl - (|10.15|) . This will be used to prove the peeling properties of the 
best decaying components of F and <fi in the main part of the paper. Although 
strictly speaking it is not necessary for our proof, we would like to know that these 
good decay estimates hold a-priori at each fixed time, without having to assume 
that the solution to p. 17(1 satisfies appropriate global decay estimates. Because of 
this we will need to cut the characteristic estimate off sharply inside the time slabs 
1 t < to. 

Lemma 10.6 (Characteristic (truncated) global Sobolev estimate). Let C(u) de- 
note the cones u — const, and let to be a fixed parameter. Now define the truncated 
cone: 

C{u) = C{u) n {1 < t < t } n {t < 2r} , 
LetT be a dyadic shell along the extended exterior region t < 2r of the cone C{u), 
and x Us smooth cutoff. Then for any test function f one has the following estimate 
for 2 < q < oo : 

(10.20) sup \xf{t,r,u)\ < r~«(l) II%X/(*>?V)IU*(S ? ) • 
Furthermore, in the case where 2 q < 4 we also have the estimate: 

(10.21) \\xf\\L?(L*(SZ)(U(u)) ~ 

r; f " f (J) (\\ X f\\ LH c(u)) + \\uL(xf)\\ L2( c iu)) + E ■ 

Proof of estimates (|10.20|l - l(10.21|l . The first estimate (|10.20(l is of course just a 
restatement of (|10.14H . Therefore, we concentrate on (|10.21fl . This estimate is 
proved using the same steps for l|10.15|) above, that is through a rescaling argument 
followed by the mixed norm Sobolev estimate i|10.18|) . However, we need to take a 
little care to deal with the sharp cutoff inside the slab {1 ^ t ^ to}. This is done 
by utilizing a extension argument for H 1 functions (see 
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The details are as follows. We first rescale all variables along the cone xC(u) by 
the factor r+(Z). Doing this, we are again reduced to showing the annular estimate 
(llO.f 7|l . However, this time it is possible for xf to extend past the boundary dA 
even though we are only allowed to estimate its behavior inside the region A. We get 
around this problem by using the fact that A is a C°° submanifold with boundary 
of R 3 , which allows us to extend xa xf to a function / with the properties: 

fU = X/U > ||/||i?i(K 3 ) ~ II xf\\w(A) ■ 

Applying estimate (|10.18|l to the function / and using this last bound yields the 
desired result. □ 
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